Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,707 C→T 6.4% intergenic (+9/‑257) mtnA → / → PP_mr25 methylthioribose‑1‑phosphate isomerase/gyrA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,975,7070CT 6.4% 1163.0 / 40.4 371intergenic (+9/‑257)mtnA/PP_mr25methylthioribose‑1‑phosphate isomerase/gyrA regulatory region
Reads supporting (aligned to +/- strand):  ref base C (170/177);  new base T (15/9);  total (185/186)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 376 total aligned reads displayed.

GCCGGGTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAAC  >  NC_002947/1975557‑1975856
                                                                                                                                                      |                                                                                                                                                      
 ccGGGTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATc                                                                                                                                                        <  2:4280941/150‑1 (MQ=255)
 ccGGGTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATc                                                                                                                                                        <  2:7505266/150‑1 (MQ=255)
   gggTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGc                                                                                                                                                      <  1:1311152/150‑1 (MQ=255)
   gggTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGc                                                                                                                                                      <  1:4622849/150‑1 (MQ=255)
        tCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCAC                                                                                                                                                <  1:1151836‑M1/150‑1 (MQ=255)
          gAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCG                                                                                                                                              <  1:7185660‑M1/150‑1 (MQ=255)
            ggTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGC                                                                                                                                            <  2:3890676‑M1/150‑1 (MQ=255)
            ggTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGC                                                                                                                                            <  1:7081940‑M1/150‑1 (MQ=255)
               ttgtgccatggatctgcgcaatgaaggccgccgtgatgggcggtcttggaaactgccggatgttatcacaaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATCGCCTCACCGGCctc                                                                                                                                         <  2:2245683‑M1/23‑1 (MQ=255)
                ttCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCT                                                                                                                                        >  1:6771126‑M1/1‑150 (MQ=255)
                ttCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCT                                                                                                                                        >  2:7007944‑M1/1‑150 (MQ=255)
                ttCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCAGCACTGATTTTGTATCGCCTCACCGGCCTCg                                                                                                                                        >  1:7019521‑M1/1‑149 (MQ=255)
                 tCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTg                                                                                                                                       >  1:2447178‑M1/1‑149 (MQ=255)
                   aaCCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCG                                                                                                                                     >  1:7387051‑M1/1‑150 (MQ=255)
                    aCCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGC                                                                                                                                    >  1:6570430‑M1/1‑150 (MQ=255)
                     cccGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCG                                                                                                                                   >  2:1197875‑M1/1‑150 (MQ=255)
                     cccGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCG                                                                                                                                   <  1:3106258‑M1/150‑1 (MQ=255)
                      ccGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGG                                                                                                                                  <  1:5616083‑M1/150‑1 (MQ=255)
                       cGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGG                                                                                                                                 >  2:4985957‑M1/1‑150 (MQ=255)
                        ggTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGC                                                                                                                                <  1:8043108‑M1/150‑1 (MQ=255)
                        ggTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGC                                                                                                                                <  2:3254915‑M1/150‑1 (MQ=255)
                           gtTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTg                                                                                                                             <  1:8325255/150‑1 (MQ=255)
                                     aCGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCa                                                                                                                   <  2:6810179/150‑1 (MQ=255)
                                       gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                 >  1:8647459/1‑150 (MQ=255)
                                       gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                 <  2:4274600/150‑1 (MQ=255)
                                         cggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTgg                                                                                                               <  1:269676/150‑1 (MQ=255)
                                         cggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTgg                                                                                                               <  2:2582767/150‑1 (MQ=255)
                                                ccTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCGCCCACAGGGATATCg                                                                                                        >  1:4065078/1‑150 (MQ=255)
                                                 cTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCAGTGGGCTGTCGc                                                                                                       >  2:8274964/1‑150 (MQ=255)
                                                 cTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGc                                                                                                       <  1:2985707/150‑1 (MQ=255)
                                                 cTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGc                                                                                                       >  2:4859380/1‑150 (MQ=255)
                                                  tGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTtgcccgcccccacaggttcaccacg                                                                                                      >  1:8785059‑M1/1‑125 (MQ=255)
                                                  tGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCa                                                                                                      >  1:6497887/1‑150 (MQ=255)
                                                   gATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCAt                                                                                                     <  1:5134648/150‑1 (MQ=255)
                                                        ggtcttggaaactgccggatgttatcacaaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGcc                                                                                                 >  2:5512317‑M1/87‑150 (MQ=255)
                                                          gtgATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTa                                                                                              <  1:7495018/150‑1 (MQ=255)
                                                            gATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAgg                                                                                            >  2:6168209/1‑150 (MQ=255)
                                                            gATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAgg                                                                                            >  2:604773/1‑150 (MQ=255)
                                                             aTCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAggg                                                                                           <  1:4498733/150‑1 (MQ=255)
                                                              tCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCGTAGGGc                                                                                          >  2:7216996/1‑150 (MQ=255)
                                                               cGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGcc                                                                                         <  1:5777928/150‑1 (MQ=255)
                                                                  gACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgtg                                                                                      <  1:2016631/150‑1 (MQ=255)
                                                                      gAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCTCCCACGGGGATATAGCATGGCCTTAGGGCCggggtgt                                                                                  >  1:2366345/1‑149 (MQ=255)
                                                                        aCGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAAc                                                                                <  1:1556116/150‑1 (MQ=255)
                                                                             gCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCtgtg                                                                           <  1:6936357/150‑1 (MQ=255)
                                                                             gCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCACGGCCCTAAGGCCatgcga                                                                                   >  1:7722133/1‑139 (MQ=39)
                                                                             gCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCACGGCCCTAAGGCCatgcga                                                                                   <  2:7722133/142‑4 (MQ=39)
                                                                               aTCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTggg                                                                         <  2:7268236/150‑1 (MQ=255)
                                                                               aTCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTACCCGCTCCCACTGGGATATCGCAAGGCCTTAGGGCCGTGGTGAACCGGTggg                                                                         >  1:5526395/1‑150 (MQ=255)
                                                                                tcgtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTtgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtggga                                                                        <  2:8785908‑M1/150‑56 (MQ=255)
                                                                                tcgtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTATGGCCGTGGTGAACCTGTgggg                                                                        <  2:5376187/150‑1 (MQ=255)
                                                                                tcgtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTgggg                                                                        <  2:5902282/150‑1 (MQ=255)
                                                                                  gtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggc                                                                      >  1:5040128/1‑150 (MQ=255)
                                                                                  gtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggc                                                                      >  2:5631830/1‑150 (MQ=255)
                                                                                  gtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTAGCGGGCTTGCCCGCGCCCACAGGGATATCGCATGGCCGTAGGGCCGTGGTGAACCggggggggg                                                                      >  1:5184291/1‑149 (MQ=255)
                                                                                   tcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcg                                                                     >  1:6303492/1‑150 (MQ=255)
                                                                                   tcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCACCCACAGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcg                                                                     >  1:3378473/1‑150 (MQ=255)
                                                                                   caaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggc                                                                      >  1:3755046‑M1/60‑150 (MQ=255)
                                                                                   caaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTCGCCCGCCCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggc                                                                      >  1:3537307‑M1/60‑150 (MQ=255)
                                                                                    cgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcgg                                                                    <  2:3565576/150‑1 (MQ=255)
                                                                                    aaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcg                                                                     >  2:4370938‑M1/59‑150 (MQ=255)
                                                                                    aaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcgg                                                                    >  1:5914745‑M1/59‑150 (MQ=255)
                                                                                     gAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcggg                                                                   >  2:944800/1‑150 (MQ=255)
                                                                                     aacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcgg                                                                    >  2:4640934‑M1/58‑150 (MQ=255)
                                                                                        gcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtTTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGa                                                                >  2:2488800‑M1/55‑150 (MQ=255)
                                                                                        cgcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTcgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaa                                                                <  2:1504566‑M1/150‑64 (MQ=255)
                                                                                        cgcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGa                                                                >  2:4889798‑M2/1‑149 (MQ=255)
                                                                                         gcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTcgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaa                                                                <  1:401508‑M1/150‑64 (MQ=255)
                                                                                          cccGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTcgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaagc                                                              <  2:1725376‑M1/150‑66 (MQ=255)
                                                                                          cccGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagc                                                              >  1:432003‑M2/1‑147 (MQ=255)
                                                                                            cGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTtgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaagccc                                                            <  1:6137441‑M1/150‑68 (MQ=255)
                                                                                            cGGATGCTGCAAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGAGCGGGCAagccc                                                            >  1:174517‑M2/1‑145 (MQ=255)
                                                                                              gATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagcccgc                                                          >  1:4815988‑M2/1‑143 (MQ=255)
                                                                                               aTGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagcccgcg                                                         <  2:2706766‑M2/150‑9 (MQ=255)
                                                                                                   tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTtgcccgcccccacaggttcaccacggccctaaggccatgcgatatcccagtgggag                                                                       <  2:2834595‑M1/132‑57 (MQ=255)
                                                                                                   tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagcccgcgaaga                                                     >  2:8426395‑M2/1‑138 (MQ=255)
                                                                                                   tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagcccgcgaaga                                                     <  2:5734451‑M2/150‑13 (MQ=255)
                                                                                                   tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCACTTAGCGGGCTTGCCCGCTCCCAAAGGGATATCGCATGGCCGTAGGGGCGAGGTGAACATGTGGGGGCGGGCGagcccgcgaata                                                     <  1:2836223‑M2/150‑13 (MQ=255)
                                                                                                    gcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCGCCCACAGGGAAATCGCAAGGCCGTAGGGCCGGGGTGAACCGGTGGGGGCGGGCGagcccgcgcagag                                                    <  1:445230‑M2/150‑14 (MQ=255)
                                                                                                      cAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGagcccgcgaagaggc                                                  >  2:6876121‑M2/1‑135 (MQ=255)
                                                                                                           tGGCTCAGTTGGTTTGCCGCATGCGCCTGAACTTATTTTTTAACGCCTCACCGGCCTGTTCGGGGGGTTGCCGGCTCCCCCGGGTATATCGCAGGGCCTTCGGGCCGGGGTGAACCTGTGGGGGCGGGCGCGCCCGCGACGAGGCCGGTg                                             <  2:3739170/150‑1 (MQ=255)
                                                                                                             gCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCa                                           >  2:2885624/1‑150 (MQ=255)
                                                                                                              cTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATGTTGTATCGCCTCACCGGCCTCTTCGCGGGGTTGCCCGCTCCCACTTGGATATCGCATGGCCTTTGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGATGAGGCCGGTGCAg                                          <  2:6410156/150‑1 (MQ=255)
                                                                                                               tCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCGTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCCGTGCAgg                                         >  1:4263176/1‑150 (MQ=255)
                                                                                                                 aGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCa                                       >  2:2778385/1‑150 (MQ=255)
                                                                                                                       tttGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACaa                                 >  2:7950669/1‑150 (MQ=255)
                                                                                                                           ccGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTAGCGGGCGTTGCCGCGCCCACGGGGCTATCGCAAGGGCGTAGGGGCGTGGTGTACCGGTGGGGGCGGGCGGGCCCCCGGACAGGCCGGTGCAGGCAAGACAAAAAc                             >  1:3170721/1‑150 (MQ=255)
                                                                                                                            cGCAAGCGCCTGCACTGTTTTTGTATTGCCTGCACCGGCCTCTTCGCGGGCTcgcccgcccccaccggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaagcccgcgaagaggccggtgaggcgatacaaaatca                             <  2:7668109‑M1/150‑99 (MQ=255)
                                                                                                                             gCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCGCCCCCGGGGATATCGCATGGCCGTAGGGCCGGGGGGAACCTGTGGGGGCGGGCGGGCCCGCGGAGAGGCCGGTGCCGGCCAGACAAAAAACa                           >  1:3406644/1‑150 (MQ=255)
                                                                                                                              cacgcgcCTGAACTTATTTGGTAGCCCCTCACCGGCCTCTTCGCGGGCTTGCCCGCACCCACTGGTTTACCGCATGGCCTTCGGGCCGTGTTGAACCTGGGGGGGCGGGCGCGCCCGCCAAGAGGCCGGTGCAGGCAATACAAAAACCAg                          <  2:4003812/147‑1 (MQ=255)
                                                                                                                              cAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAg                          >  2:4171147/1‑150 (MQ=255)
                                                                                                                                aGCGCCTGCACTGATTTTGTAGCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCAGGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGt                        <  2:8194758/150‑1 (MQ=255)
                                                                                                                                     cTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGGGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCACTACAAAAACCAGGCTTCaa                   >  2:8069774/1‑150 (MQ=255)
                                                                                                                                        cACTGATTTTGCATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAgcg                <  2:6657298/150‑1 (MQ=255)
                                                                                                                                         aCTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAgcgc               >  2:2109196/1‑150 (MQ=255)
                                                                                                                                         aCTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAgcgc               >  2:5213547/1‑150 (MQ=255)
                                                                                                                                              TTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTTTGGCCGTGGTGAACCTGTTGGGGCGGGCGAGCCCGCGAATAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTg          <  2:447845‑M1/150‑1 (MQ=255)
                                                                                                                                               aTTGTCTCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTga         >  2:1806036‑M1/2‑150 (MQ=255)
                                                                                                                                               TTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTga         >  2:189899‑M1/1‑150 (MQ=255)
                                                                                                                                                     TCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCaa   <  2:8131284‑M1/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCCGGGTGTCGAGGTGTTCAACCCGGTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAAC  >  NC_002947/1975557‑1975856

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: