Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,430 A→C 8.8% intergenic (+26/+48) sbcB → / ← PP_1366 exodeoxyribonuclease I/transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,555,4300AC 8.8% 443.5 / 24.8 152intergenic (+26/+48)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base A (69/68);  new base C (6/8);  total (75/77)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 153 total aligned reads displayed.

CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1555281‑1555576
                                                                                                                                                     |                                                                                                                                                  
cTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGCCaa                                                                                                                                                    >  2:3365508/1‑150 (MQ=255)
 ttGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                   >  2:4205662/1‑150 (MQ=255)
    ggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgtg                                                                                                                                                >  1:4803320/1‑150 (MQ=255)
       gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                             <  1:1197709/150‑1 (MQ=255)
        aCTTCGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCCGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGGTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCgg                                                                                                                                            <  2:3280476/150‑1 (MQ=255)
         cTTCGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGc                                                                                                                                           <  2:1284805/150‑1 (MQ=255)
          ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCCGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCCGGAACGCTATTCAATCGGCTGATGGCAGATCAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:1549463/1‑150 (MQ=255)
          ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:2074233/1‑150 (MQ=255)
          ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:2146142/1‑150 (MQ=255)
          ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:4048732/1‑150 (MQ=255)
           tctaaatGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCTGGCAGGCGCTTTTTGGAGGCCTGGCAACTGAATGCTCGACAATTGCAGGAACGCTATAAAATCGGCTTATGGCATATAAAAAAAAACGCCGGCAATGTGCCGGCGt                                                                                                                                         <  1:999314/147‑1 (MQ=255)
             gAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGttt                                                                                                                                       <  2:1662184/150‑1 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  2:3272077/1‑150 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       <  1:2625365/150‑1 (MQ=255)
              aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtttt                                                                                                                                      >  1:628041/1‑150 (MQ=255)
               aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCTGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:3666284/150‑1 (MQ=255)
               aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACTAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     >  2:5381052/1‑150 (MQ=255)
               aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:4024863/150‑1 (MQ=255)
                aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGCGGCCGGCCGGCGCGTGTTGGAGGCCTGGAAACGGAATGCCCGAAAATTGCAGCAACGCTTTCCAATCGGCTGTTGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTATttg                                                                                                                                    <  1:4272530/150‑1 (MQ=255)
                 aGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgt                                                                                                                                   >  2:4455961/1‑150 (MQ=255)
                  ggCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgtt                                                                                                                                  >  1:3547350/1‑150 (MQ=255)
                   gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgttt                                                                                                                                 <  1:3114319/150‑2 (MQ=255)
                     gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                               >  1:858265/1‑150 (MQ=255)
                      cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     >  1:4999663/1‑143 (MQ=255)
                      cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:4999663/143‑1 (MQ=255)
                       gcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCg                                                                                                                             >  2:1613210/1‑150 (MQ=255)
                            cAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCa                                                                                                                        <  1:5418833/150‑1 (MQ=255)
                            cAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTGTTGTGCGAAGCaa                                                                                                                       >  1:4389823/1‑150 (MQ=255)
                             aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGAGCCGGCGTTGTTGTAGTGCGAAGCaa                                                                                                                       >  2:1164641/1‑150 (MQ=255)
                                 cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                   <  2:1564093/150‑1 (MQ=255)
                                 cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                   >  2:3462713/1‑150 (MQ=255)
                                 cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAgg                                                                                                                   <  2:3753593/150‑1 (MQ=255)
                                      tgggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCACCTGCATGCTCGAAAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTa                                                                                                              <  1:1683108/149‑1 (MQ=255)
                                        ggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATg                                                                                                            >  2:2623010/1‑150 (MQ=255)
                                           gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGAATGCTCGACAATTGCAGTAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                         <  1:1965008/150‑1 (MQ=255)
                                              gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                      >  1:966907/1‑150 (MQ=255)
                                               acgaGGCTGGCTGGCGCGTGTTGGAGGCCTGGCAACTGCTTGCTCTAAATTTGCAGGAACGCTATTCAATCGGCTTATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                     <  1:3004803/150‑1 (MQ=255)
                                                cgaGGCCGGCAGGCGCGTGTTGGAGGCCTTGCAACTGCATGCTCGACAATTGCAGGATCGCTATTCAATCGGCTGATGGCAGATTTCCAAAAACGCCGGCATTGTGCCTGCGTTTTTGTTGTTCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                    <  2:3621567/150‑1 (MQ=255)
                                                 gaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATc                                                                                                   >  2:2444595/1‑150 (MQ=255)
                                                  aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                  <  2:4994200/150‑1 (MQ=255)
                                                        gcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcag                                                                                            >  1:339584/1‑150 (MQ=255)
                                                            gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCTGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        <  2:3028088/150‑1 (MQ=255)
                                                             cgcgTGTTGGAGGCCTGGCAACTGCATGCTCGACACTTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                       <  2:1498815/150‑1 (MQ=255)
                                                               cgTGATGGAGGCCTGGCAACTGCATGCACGACAATTGCACGAACGCTATACAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGtt                                                                                     <  1:5019409/150‑1 (MQ=255)
                                                                  gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc                                                                                  <  2:4314184/150‑1 (MQ=255)
                                                                    tGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCa                                                                                >  1:5124803/1‑150 (MQ=255)
                                                                     ggcggCCTGGCAACTCCAGGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAg                                                                               <  1:2276417/147‑1 (MQ=255)
                                                                       aGGCCTGGCATCTGCATGCTCGACATTTGCAGGAAGGCTATTCAATCGGCTGATGGCATATAAACAAAACCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                             <  2:3343436/150‑1 (MQ=255)
                                                                       aGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCCGTCCAGCAGGGTTGCCCAGCt                                                                             <  2:3121044/150‑1 (MQ=255)
                                                                         gCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCttt                                                                           <  2:5186304/150‑1 (MQ=255)
                                                                           cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                    >  2:4395889/1‑122 (MQ=255)
                                                                           cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                    <  1:4395889/122‑1 (MQ=255)
                                                                            tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTTCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCCGCTTTCaa                                                                        >  1:3346278/1‑150 (MQ=255)
                                                                             ggCAACTGCATGCGCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAc                                                                       <  1:315208/150‑1 (MQ=255)
                                                                                          tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                           >  1:5286278/1‑101 (MQ=255)
                                                                                            gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGGTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCa                                                        >  1:2613923/1‑150 (MQ=255)
                                                                                             aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCac                                                       >  2:4434975/1‑149 (MQ=255)
                                                                                              cAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAg                                                      >  2:2190973/1‑150 (MQ=255)
                                                                                                aTTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg                                                    <  2:2518139/150‑1 (MQ=255)
                                                                                                     aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTGTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  2:1516796/1‑150 (MQ=255)
                                                                                                       gAACGCTATTCAATCGGCTGATGGCATATAAACAAACACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACTTCGCCGCCCCAAGTGGCTTTcc                                             <  1:4727904/150‑1 (MQ=255)
                                                                                                             tATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                       >  2:3616630/1‑150 (MQ=255)
                                                                                                             tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCAAGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                       >  1:1174212/1‑150 (MQ=255)
                                                                                                                tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtttaa                                    >  1:2361963/1‑148 (MQ=255)
                                                                                                                tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    >  1:4453641/1‑150 (MQ=255)
                                                                                                                tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    >  1:4453635/1‑150 (MQ=255)
                                                                                                                tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    <  1:3789500/150‑1 (MQ=255)
                                                                                                                 cAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                   >  2:4763522/1‑150 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               <  1:4685453/137‑1 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  2:4685453/1‑137 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:4816433/150‑1 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:4841720/150‑1 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:244532/150‑1 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  >  2:4299209/1‑150 (MQ=255)
                                                                                                                       gCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCtt                             >  2:2740533/1‑150 (MQ=255)
                                                                                                                        cTGTTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                            <  1:5190196/150‑1 (MQ=255)
                                                                                                                        cTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                            <  1:3230308/150‑1 (MQ=255)
                                                                                                                            tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTAAGGGTCTTGTGGt                        >  1:2389604/1‑150 (MQ=255)
                                                                                                                            tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACTCTTTCAGGGTGTTGTGGt                        >  1:2389755/1‑150 (MQ=255)
                                                                                                                                 gcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        >  2:5239978/5‑81 (MQ=21)
                                                                                                                                 gcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        <  1:5239978/77‑1 (MQ=21)
                                                                                                                                 gcacaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccg                   >  1:4026770/5‑150 (MQ=255)
                                                                                                                                     aaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt               <  1:963987/150‑1 (MQ=255)
                                                                                                                                      aCAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt              <  1:4240732/150‑1 (MQ=255)
                                                                                                                                       cAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg             <  1:5124205/150‑1 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  1:4789947/1‑150 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTGTgg            >  2:1156082/1‑150 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGGATGCGGATAATAAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  1:196073/1‑150 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  2:1208935/1‑150 (MQ=255)
                                                                                                                                         aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGAATAATAAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCGTGTGGTTGCCGCCGTTTGt           >  1:4564892/1‑150 (MQ=255)
                                                                                                                                          aaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    <  1:3966766/124‑1 (MQ=255)
                                                                                                                                          aaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    >  2:3966766/1‑124 (MQ=255)
                                                                                                                                           aaCGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttttcgaa         <  2:1549463/150‑6 (MQ=255)
                                                                                                                                            ccGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc        <  2:4753239/149‑1 (MQ=255)
                                                                                                                                            aCGCCGGCACATTGCCGGCGTTTTTTTTGTGCGAAGCAATATGGCGTATTCGGATAATCATTACAGCAGGGTTTCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACTATTTAAGGGTCTTGTTGTTTCCGCCTTTTGtttg        >  1:1481825/1‑149 (MQ=255)
                                                                                                                                              gCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa      <  2:2146142/150‑1 (MQ=255)
                                                                                                                                               ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     <  2:2074233/150‑1 (MQ=255)
                                                                                                                                               ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     <  2:176822/150‑1 (MQ=255)
                                                                                                                                               ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     <  2:3547350/150‑1 (MQ=255)
                                                                                                                                                  gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc  <  1:1760659/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1555281‑1555576

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: