Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,431 T→A 7.9% intergenic (+27/+47) 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,4310TA 7.9% 463.0 / 21.4 149intergenic (+27/+47)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (67/70);  new base A (6/6);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 152 total aligned reads displayed.

TTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1555282‑1555576
                                                                                                                                                     |                                                                                                                                                 
ttGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCCGGCAGGCGCGTTTTGGAGGCCTGGCAACTGCATGCTCTACATTTGCAGGAACGCTATTCAATCGGCTTAGGGCAGATAAACAAACACGCCGGCAAt                                                                                                                                                   <  2:3642383/150‑1 (MQ=255)
ttGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                   >  2:4205662/1‑150 (MQ=255)
   ggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgtg                                                                                                                                                >  1:4803320/1‑150 (MQ=255)
      gACTTCGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                             <  2:117682/150‑1 (MQ=255)
      gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                             <  1:1197709/150‑1 (MQ=255)
        cTTCGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGc                                                                                                                                           <  2:1284805/150‑1 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCCGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCCGGAACGCTATTCAATCGGCTGATGGCAGATCAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:1549463/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  2:3230308/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:2074233/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:2146142/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:4994200/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                                                                                                                                       <  1:2625365/150‑1 (MQ=255)
            gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  2:3272077/1‑150 (MQ=255)
             aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtttt                                                                                                                                      >  2:1571833/1‑150 (MQ=255)
              aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCTGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:3666284/150‑1 (MQ=255)
              aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACTAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     >  2:5381052/1‑150 (MQ=255)
               aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGCGGCCGGCCGGCGCGTGTTGGAGGCCTGGAAACGGAATGCCCGAAAATTGCAGCAACGCTTTCCAATCGGCTGTTGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTATttg                                                                                                                                    <  1:4272530/150‑1 (MQ=255)
                aGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgt                                                                                                                                   >  2:4455961/1‑150 (MQ=255)
                 ggCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgtt                                                                                                                                  >  1:5138460/1‑150 (MQ=255)
                 ggCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgtt                                                                                                                                  >  1:3547350/1‑150 (MQ=255)
                  gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgttt                                                                                                                                 <  1:3114277/150‑2 (MQ=255)
                  gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgttg                                                                                                                                 >  1:4753239/1‑150 (MQ=255)
                    gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                               >  1:858265/1‑150 (MQ=255)
                     cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:4999663/143‑1 (MQ=255)
                     cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGc                                                                                                                              <  2:1373401/150‑1 (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)
                                  gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                 >  1:3205385/1‑150 (MQ=255)
                                    gAGGGCGCCGGCGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGt                                                                                                               >  2:1047906/1‑150 (MQ=255)
                                     tgggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCACCTGCATGCTCGAAAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTa                                                                                                              <  1:1683108/149‑1 (MQ=255)
                                       ggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATg                                                                                                            >  2:2623010/1‑150 (MQ=255)
                                        gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                           <  2:231460/150‑1 (MQ=255)
                                        gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                           <  1:5247658/150‑1 (MQ=255)
                                          gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGAATGCTCGACAATTGCAGTAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                         <  1:1965008/150‑1 (MQ=255)
                                             gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                      >  1:966907/1‑150 (MQ=255)
                                              acgtGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGCAACGCTATTCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                     <  1:3786253/150‑1 (MQ=255)
                                              acgaGGCTGGCTGGCGCGTGTTGGAGGCCTGGCAACTGCTTGCTCTAAATTTGCAGGAACGCTATTCAATCGGCTTATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                     <  1:3004803/150‑1 (MQ=255)
                                                 aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                  <  2:4994200/150‑1 (MQ=255)
                                                      ggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGTTGGAATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                             <  1:5060897/150‑1 (MQ=255)
                                                           gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCTGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        <  2:3028088/150‑1 (MQ=255)
                                                              cgTGATGGAGGCCTGGCAACTGCATGCACGACAATTGCACGAACGCTATACAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGtt                                                                                     <  1:5019409/150‑1 (MQ=255)
                                                                 gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc                                                                                  <  2:4314184/150‑1 (MQ=255)
                                                                 gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGATATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc                                                                                  <  2:34556/150‑1 (MQ=255)
                                                                   tGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCa                                                                                >  1:5124803/1‑150 (MQ=255)
                                                                    ggcggCCTGGCAACTCCAGGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAg                                                                               <  1:2276417/147‑1 (MQ=255)
                                                                        gCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCttt                                                                           <  2:5186304/150‑1 (MQ=255)
                                                                          cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCa                                                                         >  1:3371243/1‑150 (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)
                                                                                         tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  1:61940/1‑150 (MQ=255)
                                                                                           gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGGTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCa                                                        >  1:2613923/1‑150 (MQ=255)
                                                                                            aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCac                                                       >  2:4434975/1‑149 (MQ=255)
                                                                                            aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa                                                       <  2:2904586/150‑1 (MQ=255)
                                                                                             cAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAg                                                      >  2:2190973/1‑150 (MQ=255)
                                                                                             cAATTGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAg                                                      <  2:2091698/150‑1 (MQ=255)
                                                                                                    aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  2:1516797/1‑150 (MQ=255)
                                                                                                    aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTGTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  2:1516796/1‑150 (MQ=255)
                                                                                                      gAACGCTATTCAATCGGCTGATGGCATATAAACAAACACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACTTCGCCGCCCCAAGTGGCTTTcc                                             <  1:4727904/150‑1 (MQ=255)
                                                                                                       aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa                                            <  2:1693425/150‑1 (MQ=255)
                                                                                                          gCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActc                                         <  2:2592915/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)
                                                                                                                cAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                   >  2:4763522/1‑150 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  2:4685453/1‑137 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:244532/150‑1 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  >  2:4299209/1‑150 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:4816433/150‑1 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:4841720/150‑1 (MQ=255)
                                                                                                                           tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGt                        >  2:382874/1‑150 (MQ=255)
                                                                                                                           tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTAAGGGTCTTGTGGt                        >  1:2389604/1‑150 (MQ=255)
                                                                                                                           tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACTCTTTCAGGGTGTTGTGGt                        >  1:2389755/1‑150 (MQ=255)
                                                                                                                            ggCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTGTTGGTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCCGGGTTGCCCAGCTTTCAACCCCGGCGCCGCCCCAAGGGGCTTTCCACTCTTTCCGGGTCGTGTGGtt                       >  2:3329850/1‑150 (MQ=255)
                                                                                                                                gcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        <  1:5239978/77‑1 (MQ=21)
                                                                                                                                gcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        >  2:5239978/5‑81 (MQ=21)
                                                                                                                                gcacaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccg                   >  1:4026770/5‑150 (MQ=255)
                                                                                                                                   aaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCt                <  2:3609657/150‑1 (MQ=255)
                                                                                                                                    aaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt               <  1:963976/150‑1 (MQ=255)
                                                                                                                                     aCAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt              <  1:4240732/150‑1 (MQ=255)
                                                                                                                                     aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAACCACGGCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCgtt              >  1:4670673/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGGATGCGGATAATAAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  1:196073/1‑150 (MQ=255)
                                                                                                                                        aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGAATAATAAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCGTGTGGTTGCCGCCGTTTGt           >  1:4564892/1‑150 (MQ=255)
                                                                                                                                         aaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCCAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGGTGCCGCCTGTGGtt          >  2:2451882/1‑150 (MQ=255)
                                                                                                                                         aaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    >  2:3966766/1‑124 (MQ=255)
                                                                                                                                         aaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                    <  1:3966766/124‑1 (MQ=255)
                                                                                                                                           ccGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc        <  2:4753239/149‑1 (MQ=255)
                                                                                                                                           aCGCCGGCACATTGCCGGCGTTTTTTTTGTGCGAAGCAATATGGCGTATTCGGATAATCATTACAGCAGGGTTTCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACTATTTAAGGGTCTTGTTGTTTCCGCCTTTTGtttg        >  1:1481825/1‑149 (MQ=255)
                                                                                                                                              ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     <  2:2074233/150‑1 (MQ=255)
                                                                                                                                                 gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc  <  1:1760659/150‑1 (MQ=255)
                                                                                                                                                 gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc  <  2:4188714/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1555282‑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: