Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,513 T→C 5.2% 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 5.2% 1054.6 / 29.2 307M1A (GTG→GCG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base T (158/133);  new base C (7/9);  total (165/142)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 317 total aligned reads displayed.

AGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGG  >  NC_000913/4291367‑4291662
                                                                                                                                                  |                                                                                                                                                     
aGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                    <  1:2602607/150‑1 (MQ=255)
aGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                    <  1:1917298/150‑1 (MQ=255)
  tCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                                                                  >  1:3571497/1‑150 (MQ=255)
     ggTCGGCGTGTTGATGCTGCTCTTTTTCATTCTCTACGCCGTACAGTTGTCGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGTCGTATTCCTTCCTGAATTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGTccc                                                                                                                                               <  2:158322/150‑1 (MQ=255)
     ggTCGGCGTGATGATGCTGCGCTTTTACAGTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAccc                                                                                                                                               <  2:3449746/150‑1 (MQ=255)
      gTCGGCGTGTTGATGCTGCGCTTTTTCCTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                                                                              <  2:2043237/150‑1 (MQ=255)
        cGGCGTGTTGTTGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                            <  2:2898284/150‑2 (MQ=255)
        cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGt                                                                                                                                            <  1:3502141/150‑1 (MQ=255)
        cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                            <  2:2191772/150‑2 (MQ=255)
         ggCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                                                                           <  1:1441495/150‑1 (MQ=255)
         ggCGTGTTGATGCGGCGCTTTGTCCGTCTCGACGCCGGACAGTGGCCGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                                                                           <  2:2985545/150‑1 (MQ=255)
                       gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                             >  2:4099268/1‑150 (MQ=255)
                       gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                             >  1:4150846/1‑150 (MQ=255)
                         gcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGa                                                                                                                           >  1:4010448/1‑150 (MQ=255)
                           tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                                                                         >  2:975402/1‑150 (MQ=255)
                           tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                                                                         >  1:3031908/1‑150 (MQ=255)
                           tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                                                                         >  1:795405/1‑150 (MQ=255)
                           tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                                                                         >  1:966679/1‑150 (MQ=255)
                            ttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTg                                                                                                                        <  2:1340584/150‑1 (MQ=255)
                             tttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                                                                       >  2:1159025/1‑150 (MQ=255)
                               tCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                                                                     <  1:4040544/150‑1 (MQ=255)
                               tCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                                                                     >  1:3936279/1‑150 (MQ=255)
                               tCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                                                                     >  1:72296/1‑150 (MQ=255)
                                 aTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                                                                   >  2:1386523/1‑150 (MQ=255)
                                  ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                  <  2:866762/150‑1 (MQ=255)
                                  ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                  >  2:1669211/1‑150 (MQ=255)
                                  ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                  >  2:1305670/1‑150 (MQ=255)
                                  ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                  >  2:1790033/1‑150 (MQ=255)
                                  ttCTCTACGACGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                  >  1:1764731/1‑150 (MQ=255)
                                     tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               <  1:864733/150‑1 (MQ=255)
                                     tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               >  1:3078508/1‑150 (MQ=255)
                                     tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               <  2:2944662/150‑1 (MQ=255)
                                     tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               >  1:1570616/1‑150 (MQ=255)
                                     tctACGACGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGGTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                               >  1:551622/1‑150 (MQ=255)
                                       tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                                                                             >  1:1019323/1‑150 (MQ=255)
                                        ccGCCGGACAGTTGACGGTGGCGTCAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                                                                            <  1:2651496/149‑1 (MQ=255)
                                        aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                                                                            <  1:2059711/150‑1 (MQ=255)
                                         cGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                                                                           >  1:335858/1‑150 (MQ=255)
                                           ccGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCt                                                                                                         <  1:2141319/150‑1 (MQ=255)
                                              gACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                                                                      >  1:2086683/1‑150 (MQ=255)
                                               aCAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCa                                                                                                     <  1:291223/150‑1 (MQ=255)
                                                cAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCGTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCat                                                                                                    <  1:3410738/150‑1 (MQ=255)
                                                    tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                <  1:1234522/150‑1 (MQ=255)
                                                           ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTAGTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                                                                         <  2:3537905/150‑1 (MQ=255)
                                                             cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACt                                                                                       >  1:3844647/1‑150 (MQ=255)
                                                              gTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                      >  1:453521/1‑150 (MQ=255)
                                                                    cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                                                                                >  1:2118089/1‑150 (MQ=255)
                                                                     aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                               >  2:786217/1‑150 (MQ=255)
                                                                            agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTg                                                                        >  2:3522314/1‑150 (MQ=255)
                                                                             ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                       <  1:1386523/150‑1 (MQ=255)
                                                                              gtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                                                                      >  1:2292529/1‑150 (MQ=255)
                                                                                ggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                    >  2:1910263/1‑150 (MQ=255)
                                                                                    tCTGGACGTATTCCTTCCTGAAGTCGGTTTTGTGGGGTTGTTGTTAAGTATAGGGGTCAAAGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATTGGACTCACTTGAATAGGCATTCTGGCGCAGtt                                                                >  1:25133/1‑150 (MQ=255)
                                                                                       ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                             <  2:2386591/150‑1 (MQ=255)
                                                                                       ggACGTAGTCCTTCCTGAAGTCGGTGTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                             <  2:3171464/150‑1 (MQ=255)
                                                                                            tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                        <  1:2869522/150‑1 (MQ=255)
                                                                                              ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTccag                                                      >  1:2166303/1‑148 (MQ=255)
                                                                                              ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg                                                      <  2:1040554/150‑1 (MQ=255)
                                                                                                  ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCtgctgc                                                  <  1:477468/150‑1 (MQ=255)
                                                                                                   tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCTTTCTGGCGCAGTTCGCGCTCCTgcttct                                                 <  1:2681135/150‑1 (MQ=255)
                                                                                                    cctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                                >  2:3609386/1‑150 (MQ=255)
                                                                                                       gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGcc                                             >  2:1345485/1‑150 (MQ=255)
                                                                                                       gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTAAACGTGTTGACCCCCATGACGGCCATAGCGGGAGTGCGGTTGCGCTGGCCAGCCAAGATTCGACTAAAAATCAAAGGCATTCTGGCGCAGTTCGCGCTCCAGCAGCACGcc                                             >  1:936426/1‑150 (MQ=255)
                                                                                                        aaGTCGGTTTTCTGGCGTTGTTGTGAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCt                                            <  2:3756247/150‑1 (MQ=255)
                                                                                                        aaGTCGGTTTTCAGGCGTTGTTGTTAAGTCTCGGGGTCAAAGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCAGGCCAGCCAAGATGCGACTCACATGCATCGGCAATCAGGCGCAGTTAGCGCTCCAGCTGCTCGcca                                            >  1:2772056/1‑149 (MQ=255)
                                                                                                           tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCCTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCGTTg                                         >  1:2848628/1‑150 (MQ=255)
                                                                                                             ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                       <  1:1524677/150‑1 (MQ=255)
                                                                                                             ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                       <  1:1649467/150‑1 (MQ=255)
                                                                                                              gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATAGGCATTATGGCGCAGTTAGCGCTCCAGCTGCTCGCCGTTGGCg                                      >  1:1841122/1‑150 (MQ=255)
                                                                                                              gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCAAGATGCGACTCACTTGCATCGGCAATCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCg                                      >  1:1559748/1‑150 (MQ=255)
                                                                                                                tttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTa                                    <  2:3913160/150‑1 (MQ=255)
                                                                                                                    tGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGa                                <  2:619026/150‑1 (MQ=255)
                                                                                                                    tGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTAATGa                                >  1:598863/1‑150 (MQ=255)
                                                                                                                     ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGAc                               <  2:1080352/150‑1 (MQ=255)
                                                                                                                     ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGAc                               >  2:4051778/1‑150 (MQ=255)
                                                                                                                      gCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACg                              <  1:3723445/150‑1 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGGTGGCCTGCCATGATGCGACTCTCTTGCATCGGGTTTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                         >  1:3497467/1‑150 (MQ=255)
                                                                                                                              gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                      >  1:2389468/1‑150 (MQ=255)
                                                                                                                               ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt                     >  2:1524938/1‑150 (MQ=255)
                                                                                                                               ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt                     >  2:1198504/1‑150 (MQ=255)
                                                                                                                               ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTGTGGCGTACTGACGGAGTGtttt                     >  2:3303513/1‑150 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                    <  1:335065/150‑1 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                    <  1:1112295/150‑1 (MQ=255)
                                                                                                                                 aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtct                   <  2:2741160/150‑1 (MQ=255)
                                                                                                                                 aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtct                   <  2:70020/150‑1 (MQ=255)
                                                                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca                 <  2:3070523/150‑1 (MQ=255)
                                                                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  1:4063292/1‑150 (MQ=255)
                                                                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  2:2518348/1‑150 (MQ=255)
                                                                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              <  1:3363813/150‑1 (MQ=255)
                                                                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  1:1776028/1‑150 (MQ=255)
                                                                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca              >  2:596008/1‑150 (MQ=255)
                                                                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCGTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCGGCGCGCCTGGGGCGTACTGCCGGATTGTGGTCTCATCAg             >  2:1920576/1‑150 (MQ=255)
                                                                                                                                        ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGc            >  1:3553527/1‑150 (MQ=255)
                                                                                                                                        ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGc            >  2:1724612/1‑150 (MQ=255)
                                                                                                                                         gggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCg           >  2:1465343/1‑150 (MQ=255)
                                                                                                                                              aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc      >  2:2607150/1‑150 (MQ=255)
                                                                                                                                              aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc      >  1:3531178/1‑150 (MQ=255)
                                                                                                                                               aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct     >  2:4183209/1‑150 (MQ=255)
                                                                                                                                               aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct     <  2:95088/150‑1 (MQ=255)
                                                                                                                                                 gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg   <  2:796139/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGG  >  NC_000913/4291367‑4291662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: