Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,513 T→C 6.6% M1A (GTG→GCG) †‡ nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,5130TC 6.6% 345.7 / 11.0 122M1A (GTG→GCG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base T (59/53);  new base C (6/2);  total (65/57)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGG  >  NC_000913/4291368‑4291662
                                                                                                                                                  |                                                                                                                                                     
gTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCCTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                    <  1:107805/150‑1 (MQ=255)
gTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTg                                                                                                                                                   <  2:558672/150‑1 (MQ=255)
   tGGTCGTCGTGTTGATGATGCGCTTTTTCATTCTCTACGACGGACAGTTGACGGTGGCGTAAGCCAGAATAGAGGTGGTTTCTGGACGTATTCCTTCCTTAAGTCGG‑GTTTCTGGCGGTTTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                                <  2:834224/150‑1 (MQ=255)
   tGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCATAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                                <  2:1193658/150‑1 (MQ=255)
     gTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                                                                              <  2:1662930/150‑1 (MQ=255)
      tCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCg                                                                                                                                             >  1:1333874/1‑150 (MQ=255)
       cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGt                                                                                                                                            >  2:1541577/1‑150 (MQ=255)
       cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                            <  2:1688485/150‑2 (MQ=255)
         gCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGGGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                                                                          <  2:1520133/150‑1 (MQ=255)
            tgtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGGTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                                                                       <  1:16599/150‑1 (MQ=255)
             gtTGATGCGGCGCTTTCTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTGGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                                                                      <  2:962436/150‑1 (MQ=255)
              ttGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGc                                                                                                                                     <  2:180477/150‑1 (MQ=255)
                 aTGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                                                                                  >  1:1704823/1‑150 (MQ=255)
                   gctgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                                                                                <  2:43232/150‑1 (MQ=255)
                      gcgcTTTTTCCTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCCGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                             <  1:488582/150‑1 (MQ=255)
                      gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                             >  2:588346/1‑150 (MQ=255)
                      gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                             >  1:1576079/1‑150 (MQ=255)
                       cgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                            <  2:92207/150‑1 (MQ=255)
                       cgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGGTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                            <  1:1126640/150‑1 (MQ=255)
                       cgcTTTTTCATGCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTGTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                            <  1:1398808/150‑1 (MQ=255)
                        tcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGa                                                                                                                           <  2:1002544/149‑1 (MQ=255)
                            tttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                                                                       >  2:903882/1‑150 (MQ=255)
                            tttCATTCTCTACGCCGGACAGTTGACGGGGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                                                                       <  1:104119/150‑1 (MQ=255)
                               cATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                    >  2:338335/1‑150 (MQ=255)
                               cATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                    >  2:144320/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTTACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                >  1:1401685/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                >  1:828222/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcc                                                                                                                >  2:539643/1‑149 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                <  2:325452/150‑1 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TGTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGCCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                >  1:319463/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCATTCCTGAAGACGT‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGAAGTGCGGTTgcg                                                                                                                >  1:729386/1‑150 (MQ=255)
                                    tctACGCCGGACAGTTGACGGTGGCGTAAGCCGGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               >  1:1064387/1‑150 (MQ=255)
                                     ctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                              >  2:429022/1‑150 (MQ=255)
                                         gCCGGACAGGTGACGGGGGCGTAAGCCAGAAAAGAGGGGGGTGCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                          <  1:116684/150‑1 (MQ=255)
                                         gCCGAACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                          <  1:1691200/150‑1 (MQ=255)
                                                 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACTTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                  >  1:860275/1‑150 (MQ=255)
                                                 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                  >  2:26962/1‑150 (MQ=255)
                                                 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                  >  1:683519/1‑150 (MQ=255)
                                                 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                  >  1:8620/1‑150 (MQ=255)
                                                   tGACGGTGGCGTAAGCCAGAAAAGTGGGGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                <  2:1201608/150‑1 (MQ=255)
                                                   tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                >  2:1536696/1‑150 (MQ=255)
                                                   tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                >  1:996246/1‑150 (MQ=255)
                                                   tGACGGGGTCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                <  2:297815/150‑1 (MQ=255)
                                                    gccGGGGGCGTAAGCCAGAAAAGAGGTGGGTTCTGGACGGATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                                                                               <  2:834922/148‑1 (MQ=255)
                                                      cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCgc                                                                                             >  2:946378/1‑149 (MQ=255)
                                                       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                            <  2:241146/150‑1 (MQ=255)
                                                       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                            <  1:156375/150‑1 (MQ=255)
                                                       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTACCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGGCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                            <  1:236120/150‑1 (MQ=255)
                                                         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTCGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTc                                                                                          >  1:146185/1‑150 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                  >  2:964811/1‑150 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                  <  1:223051/150‑1 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                  <  1:801075/150‑1 (MQ=255)
                                                                  ccAGAAAAGAGGGGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTCTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATc                                                                                 <  1:75825/150‑1 (MQ=255)
                                                                    aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGCTGCGACTCACTTGCATCgg                                                                               >  1:507996/1‑150 (MQ=255)
                                                                    aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                               <  1:964811/150‑1 (MQ=255)
                                                                     gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGc                                                                              >  1:997793/1‑150 (MQ=255)
                                                                            ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                       >  2:212938/1‑150 (MQ=255)
                                                                            ggtggtTTATGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGAGTTGTTGTTAAGTCTCGGGGTCAACGTGTTTACCCCGATTACCGACTTCGCGGGAGAGCGGTTGCGCTGCCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                       >  1:1280589/1‑150 (MQ=255)
                                                                              gggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                                                                     <  2:884011/149‑1 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCGACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGAAGCGACTCACTTGCATCGGCATTCTGgcgc                                                                    >  1:934702/1‑150 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                    <  2:539210/150‑1 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                    >  2:1005665/1‑150 (MQ=255)
                                                                                    cTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc                                                               >  1:751654/1‑150 (MQ=255)
                                                                                     tGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg                                                              >  2:1463635/1‑150 (MQ=255)
                                                                                      ggACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                             <  1:213601/150‑1 (MQ=255)
                                                                                            aTTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                       >  1:1250470/1‑150 (MQ=255)
                                                                                            aTTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAACCCGTTGACGGCCTTCGCGAGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                       >  1:1028097/1‑150 (MQ=255)
                                                                                             ttccttccTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg                                                      <  2:524193/150‑1 (MQ=255)
                                                                                                cttcctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgctg                                                   >  2:1203107/1‑150 (MQ=255)
                                                                                                   cctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGTGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                                <  2:1514385/150‑1 (MQ=255)
                                                                                                   cctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                                >  1:563301/1‑150 (MQ=255)
                                                                                                      gAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGGGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGcc                                             <  1:1360736/150‑1 (MQ=255)
                                                                                                      caaGTCGG‑TTTTCTGGCGTGGTTGTTCAGGCTCGGGGTCACCGTGTTGACCCCGGTGACGGCCTTCGCGGGCGTGCGGTTGCGCTGGCCGGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGcc                                             <  1:424831/149‑1 (MQ=255)
                                                                                                         gTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCttt                                          >  1:1333397/1‑150 (MQ=255)
                                                                                                            cgtttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAGCCCGTTGACGGCCTTCGCGAGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                        <  2:1028097/147‑1 (MQ=255)
                                                                                                             g‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCACTCCTGCTGCTCGCCTTTGGCg                                      <  2:1704823/150‑1 (MQ=255)
                                                                                                                  tCTGGCGTTGTTGTTTAGTCTCGGGGTCAACGGGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACt                                  <  2:1064387/150‑1 (MQ=255)
                                                                                                                   cTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTg                                 >  2:331786/1‑150 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGTCTACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGGTGCGCTGGCCTGCCATGATGCGACTCTCTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                         <  1:239431/150‑1 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                         <  2:374795/150‑1 (MQ=255)
                                                                                                                             tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtt                       >  2:1394535/1‑150 (MQ=255)
                                                                                                                              gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                      >  1:1022907/1‑150 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttttg                    >  1:1138610/1‑149 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                    >  2:1216256/1‑150 (MQ=255)
                                                                                                                                  gtcaacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtctc                  >  1:1285820/6‑150 (MQ=255)
                                                                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCGGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca                 <  2:1062166/150‑1 (MQ=255)
                                                                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGACTGCGGTAGAGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca                 <  2:1067664/150‑1 (MQ=255)
                                                                                                                                    caacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcat                >  1:760118/4‑150 (MQ=255)
                                                                                                                                     ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc               >  2:1488551/1‑150 (MQ=255)
                                                                                                                                     ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc               <  1:1583569/150‑1 (MQ=255)
                                                                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  2:1004012/2‑150 (MQ=255)
                                                                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGAGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  1:1179462/2‑150 (MQ=255)
                                                                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAg             >  1:1256812/1‑150 (MQ=255)
                                                                                                                                        ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCGTTGGCGTACTGACGTATTGTTTTCTCATCAGc            >  1:1347154/1‑150 (MQ=255)
                                                                                                                                          ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGa          >  1:188051/1‑150 (MQ=255)
                                                                                                                                           gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAt         <  1:429022/150‑1 (MQ=255)
                                                                                                                                           gTCAACGTCTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCTCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTACGCGCTCCTGCTGCACGCCATTGGCCTACTGACGTATTGTTTTCTCATCAGCGAt         >  1:338549/1‑150 (MQ=255)
                                                                                                                                              aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCTTCAGCGGTTtc      >  2:25045/1‑150 (MQ=255)
                                                                                                                                                cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc    <  1:1461090/150‑1 (MQ=255)
                                                                                                                                                  tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCgg  <  2:39713/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGG  >  NC_000913/4291368‑4291662

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