Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,380 T→A 6.4% 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,560,3800TA 6.4% 935.8 / 25.7 284intergenic (+27/+47)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (143/120);  new base A (14/6);  total (158/126)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

TTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCA  >  NC_002947/1560231‑1560529
                                                                                                                                                     |                                                                                                                                                      
  gggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                                      <  2:1692496/150‑1 (MQ=255)
     cGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGcc                                                                                                                                                   <  2:7035501/150‑1 (MQ=255)
      gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGAGGCCGGGCACCTGCATGCTCGACAATGGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                                  <  2:3961107/150‑1 (MQ=255)
      gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                                  >  1:5823139/1‑150 (MQ=255)
         ttCGAAAAGGCGCGGCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  2:6266952/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  2:1987235/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  1:2520181/1‑150 (MQ=255)
         ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGGTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCg                                                                                                                                               >  2:3736842/1‑150 (MQ=255)
          tCGAAACGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGTGGCAGGCTGGCGCGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATGGCAGGAACGCTTTTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGt                                                                                                                                              <  1:5650022/150‑1 (MQ=255)
          acGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCCCGTGTTGGAGGCCTGGCAACTGCATGCTCGAAAATTGAACCAACGCTATTCAATCGGATGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGt                                                                                                                                              <  2:6121758/149‑1 (MQ=255)
           cGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                             <  1:992866/150‑1 (MQ=255)
           cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                             <  2:1785354/150‑1 (MQ=255)
            gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                            >  1:5111423/1‑150 (MQ=255)
            gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                            >  1:2571247/1‑150 (MQ=255)
            gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                            >  1:1022036/1‑150 (MQ=255)
             aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtttt                                                                                                                                           >  2:4977650/1‑150 (MQ=255)
              aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGttttt                                                                                                                                          <  2:2420237/150‑1 (MQ=255)
              aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                          >  2:3164817/1‑150 (MQ=255)
                    gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                                    >  1:4645390/1‑150 (MQ=255)
                    gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                                    >  1:4116556/1‑150 (MQ=255)
                    gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATGGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTTTttgtttat                                                                                                                                    <  2:3589174/150‑4 (MQ=39)
                    gcgcCAGCAAGCCGGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTtgtg                                                                                                                                    >  2:6524132/1‑150 (MQ=255)
                      gcCAGCACGCCGGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTGGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCg                                                                                                                                  <  2:954554/150‑1 (MQ=255)
                            aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCac                                                                                                                            >  2:5152662/1‑149 (MQ=255)
                            aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCaa                                                                                                                            >  1:776963/1‑150 (MQ=255)
                             aGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCCCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAag                                                                                                                           >  1:1712405/1‑150 (MQ=255)
                               ccTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAagag                                                                                                                         >  1:3677662/1‑150 (MQ=255)
                                cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                        <  1:5499611/150‑1 (MQ=255)
                                cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                        <  2:1745035/150‑1 (MQ=255)
                                  gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                      <  1:432148/150‑1 (MQ=255)
                                      gggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTAt                                                                                                                  <  2:4019693/150‑1 (MQ=255)
                                        gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                                <  1:7404644/150‑1 (MQ=255)
                                           ccgccgCGGCCGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa                                                                                                             <  2:445451/150‑1 (MQ=255)
                                           cccacgaGGCAGGCAGGCGCGTGTTGGAGGCCGGCCAACTGCATGCTCGACAATTGCAGGAACGCTATGCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa                                                                                                             <  1:3880456/147‑1 (MQ=255)
                                           ccGACGATGCAGGCAGGCGCGTGTTGGAGGCCTCGCAACTGCATGCTCGACAATTCCAGGAACGCTATTCAATCGGCTGATGGAAGAAAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa                                                                                                             <  2:6899451/150‑1 (MQ=255)
                                            cgacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGAt                                                                                                            >  2:6764729/1‑150 (MQ=255)
                                             gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                           >  1:3606448/1‑150 (MQ=255)
                                                gaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATc                                                                                                        >  2:2667847/1‑150 (MQ=255)
                                                 aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  1:2717676/150‑1 (MQ=255)
                                                 aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  2:3828337/150‑1 (MQ=255)
                                                 aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGATATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  1:2284370/150‑1 (MQ=255)
                                                   gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt                                                                                                     >  2:6740781/1‑150 (MQ=255)
                                                    ctggctggcGGGTTTTGGAGGCCTGGCAACTGCTTGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                                    <  1:6376913/148‑1 (MQ=255)
                                                       gcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcag                                                                                                 >  2:25352/1‑150 (MQ=255)
                                                         aggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagca                                                                                               >  1:4349410/1‑150 (MQ=255)
                                                         acgcgcGTGTTCGAGGCCTGCCAACTGCATGCTCCACAATAGCAGGAACGCTATACAATCGGCTGATGGCACATAAACAAAAACGCCGGCAATGTGCCGGCGATATTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagca                                                                                               <  2:5686319/148‑1 (MQ=255)
                                                           gcgcGTGTTGGAGGCCTGGCTAGTGCATGCTCGACAATGGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                             <  2:6022332/150‑1 (MQ=255)
                                                            cgcgTTTTGGAGGCCTGGCAACTGCATGCCCGACAATTGCAGGAACGCTATACAATCGGCTGATGGCAGAAAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                            <  1:133501/150‑1 (MQ=255)
                                                             gcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagctcagc                                                                                           >  1:2209051/1‑145 (MQ=255)
                                                               gtgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAATTCGCCGGCAATGTGCCGGCGTTTTTGTTGTGGGAAGGAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTg                                                                                         >  1:5884287/1‑150 (MQ=255)
                                                                    ggAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCTCAg                                                                                    <  1:573707/150‑1 (MQ=255)
                                                                     gAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGc                                                                                   >  1:1257719/1‑150 (MQ=255)
                                                                       ggCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCtt                                                                                 <  1:6899766/150‑1 (MQ=255)
                                                                         ccTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTc                                                                               >  1:1461229/1‑150 (MQ=255)
                                                                          cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGTGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCTGCTTTct                                                                              >  2:7122226/1‑149 (MQ=255)
                                                                           tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                             >  2:1248424/1‑150 (MQ=255)
                                                                                aCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACg                                                                        >  2:4207334/1‑150 (MQ=255)
                                                                                 cTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                       >  2:1885858/1‑150 (MQ=255)
                                                                                  tGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc                                                                      <  1:1773794/150‑1 (MQ=255)
                                                                                       gCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgccgc                                                                >  1:198154/1‑150 (MQ=255)
                                                                                        cTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCgccgcc                                                               <  2:2310680/150‑1 (MQ=255)
                                                                                           gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCa                                                            >  1:3645745/1‑150 (MQ=255)
                                                                                            aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGTTTTCAACCACGTCGCCGCCCCaa                                                           >  1:505738/1‑150 (MQ=255)
                                                                                            aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa                                                           <  2:730151/150‑1 (MQ=255)
                                                                                               aTTGCAGGAACGATATTCAATCGGCTGTTGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg                                                        <  1:3394800/150‑1 (MQ=255)
                                                                                                ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAg                                                                                                                               <  2:7172345/79‑1 (MQ=255)
                                                                                                    aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCTCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                   >  2:4436918/1‑150 (MQ=255)
                                                                                                    aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCGTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                   >  2:7097575/1‑150 (MQ=255)
                                                                                                    aGGAACGCGATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                   <  1:7042634/150‑1 (MQ=255)
                                                                                                      gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTTTGCGAAGCAATAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCAAGATTTCAACCACGTCGCCGCCCCAAGTGGCtttac                                                 >  2:5006381/1‑148 (MQ=255)
                                                                                                       aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa                                                <  2:684963/150‑1 (MQ=255)
                                                                                                          gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCActc                                             >  1:1633997/1‑150 (MQ=255)
                                                                                                            tATTCAATCGGCTGATGGCAGATAAAGAAAAACGCCGGCAATGTGCCGGCGTTTGGGTTGTGCGAAGCAAGTGGGCGTATGGGGATAATCAGTCCAGCAGGGTGGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCtt                                           >  1:2960119/1‑150 (MQ=255)
                                                                                                              ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc                                         >  2:1269176/1‑150 (MQ=255)
                                                                                                               tCATTCGTCTGAGGGCATATAAACAAAAACGCCGGCAATGTGCCGGCTTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                        <  1:5006381/150‑1 (MQ=255)
                                                                                                               tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                        <  1:5006072/150‑1 (MQ=255)
                                                                                                               tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCa                                        <  2:2755104/150‑1 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAgg                                      >  1:2177156/1‑150 (MQ=255)
                                                                                                                 aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTGTCAgg                                      >  2:5448024/1‑150 (MQ=255)
                                                                                                                      gCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCAAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCtt                                 <  2:1423739/150‑1 (MQ=255)
                                                                                                                       cTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                                <  2:5823139/150‑1 (MQ=255)
                                                                                                                        tGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCGACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt                               >  1:1131945/1‑150 (MQ=255)
                                                                                                                         taTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                              <  1:5581299/149‑1 (MQ=255)
                                                                                                                         gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTtgtg                              <  1:4478240/150‑1 (MQ=255)
                                                                                                                          aTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCTAAGTGGCTTTCCACTCTTTCAGGGTCTTGTgg                             >  2:2252794/1‑150 (MQ=255)
                                                                                                                             ttcgcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCTCAAGTGGCTTTCCACTCTGTCAGGGTCTTGTGGTTg                          >  2:4082167/8‑150 (MQ=38)
                                                                                                                              cAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATACTCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCTCTCTTTCAGGGTCTGGTGGTTgc                         <  2:6408537/150‑1 (MQ=255)
                                                                                                                               aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGCTAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCGTTCAGGGTCTTGTGGTTgcc                        >  1:3284194/1‑150 (MQ=255)
                                                                                                                               aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                        <  2:3813228/150‑1 (MQ=255)
                                                                                                                               aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                        >  1:5520843/1‑150 (MQ=255)
                                                                                                                                     aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGGGGTTGCCGCCttt                  >  1:1577798/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGGCTTGTGGTTGCCGCCTTTgg                >  1:481604/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                                  <  1:3971824/69‑1 (MQ=255)
                                                                                                                                           aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGTGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTc            >  1:7174995/1‑150 (MQ=255)
                                                                                                                                           aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt             >  2:3236014/1‑150 (MQ=255)
                                                                                                                                           aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCCCTATTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTc            >  2:6311579/1‑150 (MQ=255)
                                                                                                                                           aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCTCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            >  2:4686075/1‑150 (MQ=255)
                                                                                                                                           aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTc            >  1:7134309/1‑150 (MQ=255)
                                                                                                                                                ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTAGATCa       >  1:535255/1‑150 (MQ=255)
                                                                                                                                                 gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc      >  2:6669324/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCA  >  NC_002947/1560231‑1560529

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: