Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,379 A→C 9.7% 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,560,3790AC 9.7% 615.4 / 37.9 223intergenic (+26/+48)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base A (107/92);  new base C (14/9);  total (122/101)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1560230‑1560525
                                                                                                                                                     |                                                                                                                                                  
cTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCaa                                                                                                                                                    >  1:1013463/1‑150 (MQ=255)
 ttGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                   >  2:401168/1‑150 (MQ=255)
 ttGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                   >  2:4516118/1‑150 (MQ=255)
   gggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                                 >  2:4993954/1‑150 (MQ=255)
        aCTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCgg                                                                                                                                            <  2:2108723/150‑1 (MQ=255)
         cTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGc                                                                                                                                           >  1:4402818/1‑150 (MQ=255)
          ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:1736219/1‑150 (MQ=255)
           tctaaaaGGCGCGCCAGCACGCCTGGGTGGGCGCCGACGCGGCCGGCCGGGGGGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATTGCAGCAACGCTTTTCAAGCGGCTGTTGGCTTATAAACAAACACGCCGGCAATGTGCCGGCGt                                                                                                                                         <  1:1427578/147‑1 (MQ=255)
           tCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGt                                                                                                                                         >  2:2709769/1‑150 (MQ=255)
           tCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGAAACTGCATGATCGACATTTGCAGTAACGCTATTCAATCGGCTGTTGGAAGATAAACAAAAACGCCGGCACATTGCCGGCGt                                                                                                                                         <  2:1716410/150‑1 (MQ=255)
            cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGtt                                                                                                                                        <  2:2245878/150‑1 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCTGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       <  1:3648529/150‑1 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  1:4256182/1‑150 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  2:240859/1‑150 (MQ=255)
             gAAAAGGCGCGCCAGCAAGCCGGGGCGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCACCTGCCTGCTCGACACTTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCACATTGCCGGCGttt                                                                                                                                       <  2:3294687/150‑1 (MQ=255)
               acaGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCCGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAAGCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  1:5399706/148‑1 (MQ=255)
               aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:5236962/150‑1 (MQ=255)
                 aGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgt                                                                                                                                   <  1:2709769/150‑1 (MQ=255)
                  ggCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttgtt                                                                                                                                  <  1:2778957/150‑1 (MQ=255)
                   gcgcgcCAGCACGCCTGGGAGGGCGCCGACGAGGCTGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATCCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTTTttgttt                                                                                                                                 <  1:4105731/150‑2 (MQ=39)
                    cgcgcCCGCACGCCTGGGAGGGCGCCGACGAGGCAGGCTGGCGCGTGTGGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTtgttgt                                                                                                                                <  1:3938527/150‑1 (MQ=255)
                    cgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTtggtgt                                                                                                                                >  1:5091601/1‑150 (MQ=255)
                     gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTCTTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCACTGTGCCGGCGTTTTTGTtgtg                                                                                                                               >  1:5321762/1‑150 (MQ=255)
                     gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                               >  2:4644470/1‑150 (MQ=255)
                         cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                           <  2:4214695/150‑1 (MQ=255)
                          aGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAg                                                                                                                          >  1:2518349/1‑150 (MQ=255)
                           gcccgCCTGGGAGGGCGCCGCCGCGGCAGGCCGGCGCGGGTTGGAGGCCGGGCAACGGCATGCGCGACAATGGCCGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGc                                                                                                                         <  1:5052473/146‑1 (MQ=255)
                           gCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGAAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGc                                                                                                                         <  2:712249/150‑1 (MQ=255)
                              aGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCTACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTGTTTGTTGTGAGAAGCAag                                                                                                                      >  2:4433017/1‑150 (MQ=255)
                               gCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAaga                                                                                                                     >  1:969155/1‑150 (MQ=255)
                                ccTGGGTGGGCGCCTACGCGGCAGGCTGGCGCTTTTTGTTGGCCTGGCAAGGGCTGGCTCTACAATGGCTGTAACGCTATTCAATCGGCTGATGGCAGATAAAAACAAACGCCGTCAATGTGCCGGCGATTTTGTTGTGCGAAGCattag                                                                                                                    <  2:1560206/150‑5 (MQ=255)
                                 cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                   >  2:3941147/1‑150 (MQ=255)
                                 cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATGGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                   <  2:2527556/150‑1 (MQ=255)
                                  tGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGTAACGCTATACAATCGGCTGTTGGCAGATAAAAAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAggg                                                                                                                  <  1:1135849/150‑1 (MQ=255)
                                   gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGTTAAACAATATCGCCGGCAATGTGGGGGCGTTTTTGTTGTGGGAAGCAAGAgggg                                                                                                                 >  1:5075808/1‑149 (MQ=255)
                                   gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACACAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                 <  2:3237901/150‑1 (MQ=255)
                                   gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                 <  1:3224749/150‑1 (MQ=255)
                                      aGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTa                                                                                                              <  2:3524185/150‑1 (MQ=255)
                                         gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                           >  1:2165573/1‑150 (MQ=255)
                                            ccGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGAAAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa                                                                                                        <  1:3945941/150‑1 (MQ=255)
                                             cgacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGAt                                                                                                       <  2:832724/150‑1 (MQ=255)
                                                 gcggctggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCCGGAACGCTATTCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATc                                                                                                   <  1:2417869/148‑1 (MQ=255)
                                                 gaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATc                                                                                                   >  1:4154143/1‑150 (MQ=255)
                                                   ggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg                                                                                                 <  1:1303230/150‑1 (MQ=255)
                                                    gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt                                                                                                >  1:3135725/1‑150 (MQ=255)
                                                      aggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTcc                                                                                              <  2:3643878/150‑1 (MQ=255)
                                                        gcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTGTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcag                                                                                            >  1:4267188/1‑150 (MQ=255)
                                                           ggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagcag                                                                                         >  1:1059947/1‑150 (MQ=255)
                                                                   ttGGAGGCCTGGCAAGTGCATGCTCGACAATTGCTGGAACGCTATTCAATCGGCTGATGGGAGTTAAACAATTTCGCCGGCAATGTGCCGGGGTTTTTGTTGTGGGTAGGAAGTGGGCGTTTGCGGTTTATGAGTCCAGGAGGGTTgggg                                                                                 >  2:4822771/1‑147 (MQ=255)
                                                                   ttGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGccc                                                                                 >  2:102577/1‑150 (MQ=255)
                                                                    tGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCa                                                                                <  1:3554145/150‑1 (MQ=255)
                                                                           cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAATTCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCa                                                                         >  2:493680/1‑150 (MQ=255)
                                                                           cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCa                                                                         >  1:1524107/1‑150 (MQ=255)
                                                                            tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGATAATCAGTACAGCAGGGTTGCCCACCTTTCaa                                                                        >  1:563095/1‑150 (MQ=255)
                                                                             ggCAACTGCATGCTCGACAATTGCAGGAACGCTATTCACTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAc                                                                       <  2:2062829/150‑1 (MQ=255)
                                                                             ggCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAc                                                                       <  2:601038/150‑1 (MQ=255)
                                                                              gCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAcc                                                                      >  1:4510218/1‑150 (MQ=255)
                                                                                  cTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                  <  1:872958/150‑1 (MQ=255)
                                                                                   tGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc                                                                 <  2:5427464/150‑1 (MQ=255)
                                                                                    gCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcg                                                                >  2:5129882/1‑150 (MQ=255)
                                                                                     cATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgc                                                               <  1:4452087/150‑1 (MQ=255)
                                                                                         cTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCgccgcc                                                           <  2:5057054/150‑1 (MQ=255)
                                                                                         cTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGGATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCgccgcc                                                           >  1:779344/1‑150 (MQ=255)
                                                                                          tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  2:3553835/1‑150 (MQ=255)
                                                                                          tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  2:2653428/1‑150 (MQ=255)
                                                                                               aaTTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGt                                                     <  1:4993954/150‑1 (MQ=255)
                                                                                                 ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTgg                                                   <  2:2970459/150‑1 (MQ=255)
                                                                                                     aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                               >  1:2375146/1‑150 (MQ=255)
                                                                                                       gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                             <  1:5468701/150‑1 (MQ=255)
                                                                                                          cGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCAct                                          >  2:2789264/1‑150 (MQ=255)
                                                                                                           gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActc                                         >  2:2489172/1‑150 (MQ=255)
                                                                                                           gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActc                                         >  2:3191508/1‑150 (MQ=255)
                                                                                                             tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                       >  1:1921117/1‑150 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:3827255/150‑1 (MQ=255)
                                                                                                                  aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  <  2:1751975/150‑1 (MQ=255)
                                                                                                                     cGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTc                               <  1:2170320/150‑1 (MQ=255)
                                                                                                                       gCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCtt                             >  2:2366279/1‑150 (MQ=255)
                                                                                                                          gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGGCTtgtg                          >  1:1629390/1‑150 (MQ=255)
                                                                                                                           ttGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTgg                         <  1:2589509/149‑1 (MQ=255)
                                                                                                                              gCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTg                      <  1:4951636/150‑1 (MQ=255)
                                                                                                                                aGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                    >  1:2696376/1‑150 (MQ=255)
                                                                                                                                   tAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCCGTCCAGCAGGGTTGCCCAGCATTCAACCACGTCGCCGCCCCAAGTGGCGGTCCACTCGTTCAGGGTCTTGTGGTTgccgcc                 >  2:3873288/1‑150 (MQ=255)
                                                                                                                                     aaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGTCtt               <  2:1633722/150‑1 (MQ=255)
                                                                                                                                      aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGACGCAAGAGGGCGTATGCGGATAATCAGTCCCGCCGGGTTGCCCCGCTTTCAACCACGGCGCCGCCCCACGTGGCTTGCCACTCGTGCCGGGGCTGGGGGGTgccgccgtg              >  2:1419629/1‑147 (MQ=255)
                                                                                                                                      aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                             <  2:194762/119‑1 (MQ=255)
                                                                                                                                      aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                             >  1:194762/1‑119 (MQ=255)
                                                                                                                                       ctcataCGGCGGCTATGTGCCGGGGTTTTTGTTGTGCGAAGCTTGTGGGCGTATGCGGATTTTCTGTCGAGCAGGGTTGCCCAGCTTTCATCCTCGTCGCCGCCCCAAGTGGGTTTCGACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg             <  1:175446/145‑1 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCAATGTGCCGGTGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  1:2953298/1‑150 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  2:4905673/1‑150 (MQ=255)
                                                                                                                                        aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCGTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  2:3103211/1‑150 (MQ=255)
                                                                                                                                         aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt           >  1:2246854/1‑150 (MQ=255)
                                                                                                                                         aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCCCCCCAACTGGCTTTCCACTCTTTCAGGGGCGTGTGGTTGCCGCCGTTGGt           >  1:3780691/1‑150 (MQ=255)
                                                                                                                                            aCGCCGGCACATTGCCGGCGTTTTTGTTTTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCACCAGGGTTTCCCAGCGTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACAATTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTc        >  1:2059320/1‑150 (MQ=255)
                                                                                                                                            aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAAACACGTCGCCGCCCCAAGTGGCTTTCCAATATTTAAGGGTCTTGTGGTTGACGCCATTTGttta        >  1:5227871/1‑149 (MQ=255)
                                                                                                                                            aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc        <  1:2406886/150‑1 (MQ=255)
                                                                                                                                             cGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCg       >  2:460351/1‑150 (MQ=255)
                                                                                                                                             cGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTCg       >  2:889069/1‑150 (MQ=255)
                                                                                                                                                 ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   >  2:1141265/1‑150 (MQ=255)
                                                                                                                                                 ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   >  1:1349718/1‑150 (MQ=255)
                                                                                                                                                 ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   >  1:4567836/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1560230‑1560525

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: