Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,046,787 A→C 6.2% intergenic (‑185/‑179) yihG ← / → polA inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,046,7870AC 6.2% 334.8 / 8.5 113intergenic (‑185/‑179)yihG/polAinner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base A (53/52);  new base C (2/5);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 117 total aligned reads displayed.

TAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTAC  >  NC_000913/4046638‑4046936
                                                                                                                                                      |                                                                                                                                                      
tAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATTAAAAACAAAGACAAACA‑TCATGCTGTaaaa                                                                                                                                                        <  2:192631/150‑1 (MQ=255)
  tcAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAg                                                                                                                                                      <  1:886344/149‑1 (MQ=255)
  gCAGAAAACGACCCAAATAACGGATGATCCTTAAGGTGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAATAACAAAGACAAACA‑TCATGCTGTATAAAg                                                                                                                                                      <  2:1090537/150‑1 (MQ=255)
    aGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCa                                                                                                                                                    >  2:334645/1‑150 (MQ=255)
     taaaaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCAt                                                                                                                                                   <  1:744643/149‑1 (MQ=255)
        aaCGACCCAAATAACGTATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTTTTTACAGCATGat                                                                                                                                                <  1:907765/150‑1 (MQ=255)
        aaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGat                                                                                                                                                <  2:721123/150‑1 (MQ=255)
         aCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTTTTTACAGCATGatg                                                                                                                                               <  2:1004979/150‑2 (MQ=255)
         aCGAACCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGata                                                                                                                                               <  2:348156/150‑1 (MQ=255)
           gACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGataat                                                                                                                                             >  1:195555/1‑150 (MQ=255)
           caCCCAAATAACGTATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTAAAAGGGATGGATTTATTTTTAACTGCATGAAAAACACAGACAAACA‑TCATGCTTTTTACAGCATGatgtt                                                                                                                                             <  1:185947/149‑4 (MQ=255)
             cccAAATAACTGATTATCCTTAATGAGAAAAATAATTCATATCTATACACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTTTTTACAGCATGatgtttg                                                                                                                                           <  2:622664/150‑6 (MQ=255)
             cccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTTTTTACAGCATGatgtttg                                                                                                                                           <  2:1314939/150‑6 (MQ=255)
              ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTATAAAAAATCGCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑aaa                                                                                                                                         <  2:826168/150‑1 (MQ=255)
                aaaTAACTGATGATCCTTAATGAGAATAATAATTCATCTATAACCACATTAGAAAAAATCACATTATCTCAATTATTAGGGATGGATTTATTTATAACTGCAGGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAAtt                                                                                                                                       <  1:201097/150‑1 (MQ=255)
                aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAAtt                                                                                                                                       >  2:1042058/1‑150 (MQ=255)
                 aaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTa                                                                                                                                      <  2:113725/150‑1 (MQ=255)
                    aaCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTaaaa                                                                                                                                   >  2:744475/1‑150 (MQ=255)
                    aaCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTA‑TTTTAAAAGGAAGTAAAACAAATGCAAACATTAATGC‑GTAAAAAGGATGATAATAAAATTaaaa                                                                                                                                   >  1:467796/1‑150 (MQ=255)
                      cGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGc                                                                                                                                 <  2:265747/150‑1 (MQ=255)
                       ggATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCg                                                                                                                                >  1:1042339/1‑150 (MQ=255)
                       ggATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCg                                                                                                                                >  2:649791/1‑150 (MQ=255)
                           gatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGt                                                                                                                            >  2:882500/1‑150 (MQ=255)
                            atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAATAGCATGATAAT‑AAATTAAAAGCGATGTa                                                                                                                           >  1:54432/1‑150 (MQ=255)
                             tCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTaa                                                                                                                          >  2:44841/1‑150 (MQ=255)
                              ccTTAAGTATAAAAATAATTATTATCTATCAACATTAGAAAAAATACCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTaaa                                                                                                                         <  1:828053/150‑1 (MQ=255)
                               cTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTTTTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATATT‑AAATTAAAAGCGATGTAaat                                                                                                                        >  2:482029/1‑150 (MQ=255)
                               cTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAATAAGCATGATAAT‑AAATTAAAAGCGATGTAaat                                                                                                                        >  2:122343/1‑150 (MQ=255)
                               cTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAaat                                                                                                                        <  1:451651/150‑1 (MQ=255)
                               attAAGGAGAATAATAATTCATATCTATCCAAATTAGAAAAAATCCCATTATCTCAATTATTAGGTATTGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAaat                                                                                                                        <  1:15005/149‑1 (MQ=255)
                                 tAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATAATTAGGGATTGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCTTGCTGTAAAAAGCATGATTAT‑AAATTAAAAGCGATGTAaataa                                                                                                                      >  2:1075356/1‑150 (MQ=255)
                                      agaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTAt                                                                                                                 <  2:961013/150‑1 (MQ=255)
                                       gaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATg                                                                                                                >  1:273702/1‑150 (MQ=255)
                                       gaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATg                                                                                                                >  2:479157/1‑150 (MQ=255)
                                                    tatctatcCACATTAGAAAAACTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATATTTTATGCACAAAGTTATcc                                                                                                   <  2:1116429/150‑1 (MQ=255)
                                                    tatctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATcc                                                                                                   <  1:406116/150‑1 (MQ=255)
                                                    tatctatcAACATTAGAAAAAATCACATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATTCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATcc                                                                                                   <  2:1121614/150‑1 (MQ=255)
                                                     atctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCca                                                                                                  >  1:205218/1‑150 (MQ=255)
                                                      tctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCcac                                                                                                 <  2:195555/150‑1 (MQ=255)
                                                        tatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACAt                                                                                               >  2:963078/1‑150 (MQ=255)
                                                            cacaTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACg                                                                                           >  2:1251905/1‑150 (MQ=255)
                                                            cacaTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACg                                                                                           >  2:1186572/1‑150 (MQ=255)
                                                            cacaTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACg                                                                                           >  2:702801/1‑150 (MQ=255)
                                                              caTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAt                                                                                         >  2:591846/1‑150 (MQ=255)
                                                                 tAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTa                                                                                                      <  1:340391/134‑1 (MQ=255)
                                                                 tAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTa                                                                                                      >  2:340391/1‑134 (MQ=255)
                                                                     aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                                  <  1:633831/150‑1 (MQ=255)
                                                                     aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                                  <  2:520173/150‑1 (MQ=255)
                                                                      aaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgagc                                                                                 <  1:1084007/150‑1 (MQ=255)
                                                                         aTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TAATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCCCATGACGATTTGAGAGCGAt                                                                              >  2:352796/1‑150 (MQ=255)
                                                                           cccATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATcc                                                                            >  1:751640/1‑150 (MQ=255)
                                                                                 aTCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCagaaga                                                                      >  2:280379/1‑150 (MQ=255)
                                                                                   ctcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATc                                                                    >  2:887747/1‑150 (MQ=255)
                                                                                    tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                                   <  1:629518/150‑1 (MQ=255)
                                                                                    tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGAGAGCGATCCAGAAGATCt                                                                   >  1:349195/1‑150 (MQ=255)
                                                                                     cAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTa                                                                  >  2:1221055/1‑150 (MQ=255)
                                                                                      aaTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                                 <  1:284761/150‑1 (MQ=255)
                                                                                       attattAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACa                                                                <  1:891160/150‑1 (MQ=255)
                                                                                                   tgatTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCAc                                                    <  1:649791/149‑1 (MQ=255)
                                                                                                   ggatTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCAc                                                    >  1:345200/1‑150 (MQ=255)
                                                                                                    gatTTATTTTTAACTGCATGAAAAACAAAGCCAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATATATGCACAAAGTGATCCACATGACGATTTGCGAGCGCTCCAGAAGATCTACAAAAGATTTTCACg                                                   >  2:205364/1‑150 (MQ=255)
                                                                                                        tatttTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGaaaa                                               >  1:397748/1‑150 (MQ=255)
                                                                                                          tttttAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                             <  2:482937/150‑1 (MQ=255)
                                                                                                          tttttAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                             <  2:751640/150‑1 (MQ=255)
                                                                                                          tttttAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                             <  1:711944/150‑1 (MQ=255)
                                                                                                               aaCTGCATGAAAAACAAAGACAAAAA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTATAAGCGATGTAAATATTTTATGCAAAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAATAGATTTTCACGAAAAGCGGTGa                                        <  1:787335/150‑1 (MQ=255)
                                                                                                                 cTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGAGATGTAAATAATTTATGCACAAAGTTATCAACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaa                                      >  1:571038/1‑150 (MQ=255)
                                                                                                                  tGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaa                                     >  2:1133248/1‑150 (MQ=255)
                                                                                                                  ggCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAAATAAAAGCGATGTAAATAATATATGCACAACGTTATCCACATGACGAATTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaa                                     >  1:492551/2‑150 (MQ=255)
                                                                                                                   gCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaaa                                    >  2:663929/1‑150 (MQ=255)
                                                                                                                       taaaaaCAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCa                                <  1:629425/149‑1 (MQ=255)
                                                                                                                       gAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCa                                >  1:1096525/1‑150 (MQ=255)
                                                                                                                       gAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCAACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAAATCa                                >  2:37810/1‑150 (MQ=255)
                                                                                                                        aaaaaCAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                               >  2:1218057/1‑150 (MQ=255)
                                                                                                                        aaaaaCAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                               >  2:789806/1‑150 (MQ=255)
                                                                                                                           aaCAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGATGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGtt                            >  2:857926/1‑150 (MQ=255)
                                                                                                                           aaCAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGtt                            >  2:1135976/1‑150 (MQ=255)
                                                                                                                                aGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGAAGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCAt                       <  2:104368/150‑1 (MQ=255)
                                                                                                                                  acacacA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATGTTTGCACAACGTTCTCCCCATGACGATGTGCGCGCGATCCAGCAGCTCTACAAAAGAGTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATcc                     <  1:657312/150‑1 (MQ=255)
                                                                                                                                  aCAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAATATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATcc                     <  2:225613/150‑1 (MQ=255)
                                                                                                                                  aCAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATcc                     <  2:105839/150‑1 (MQ=255)
                                                                                                                                  aCAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATcc                     <  1:455483/150‑1 (MQ=255)
                                                                                                                                  aCAAACA‑TCAAGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATATATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGACGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATcc                     <  2:951587/150‑1 (MQ=255)
                                                                                                                                   cAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCct                    >  1:10911/1‑150 (MQ=255)
                                                                                                                                     atta‑tcatGCTTTTTACAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgt                  >  1:422096/4‑150 (MQ=255)
                                                                                                                                      aCA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtc                 <  1:871624/150‑1 (MQ=255)
                                                                                                                                       ca‑tcatGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtct                <  1:397515/150‑1 (MQ=255)
                                                                                                                                       ca‑tcatGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtct                <  1:924280/150‑1 (MQ=255)
                                                                                                                                        a‑tcatGCTTTTTACAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtctg               >  1:401606/1‑150 (MQ=255)
                                                                                                                                          tcatGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtctgt              >  1:398690/1‑150 (MQ=255)
                                                                                                                                            atGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATTACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTgg            <  2:354595/150‑1 (MQ=255)
                                                                                                                                            atGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTgg            <  1:334645/150‑1 (MQ=255)
                                                                                                                                            atGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTgg            <  1:1113772/150‑1 (MQ=255)
                                                                                                                                             tGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGc           <  1:702801/150‑1 (MQ=255)
                                                                                                                                               cTGTAAAAAGTATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAATTTATCCAGTTGACTATTTGCTAGCGATCCAGAATATCTACTAAATATTTTCACTAAAAGCGGGGTAAAACTCATGTTTTCATCCTGTCTGTGTCAt         <  1:704635/150‑1 (MQ=255)
                                                                                                                                                 gTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATcc       <  2:516413/150‑1 (MQ=255)
                                                                                                                                                 gTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGAACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATcc       >  1:1182500/1‑150 (MQ=255)
                                                                                                                                                    aaaaGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCttt    <  2:611140/150‑1 (MQ=255)
                                                                                                                                                    aaaaGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCttt    >  1:458043/1‑150 (MQ=255)
                                                                                                                                                      aaGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGAGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTAc  <  1:146099/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACA‑TCATGCTGTAAAAAGCATGATAAT‑AAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTAC  >  NC_000913/4046638‑4046936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑