Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,970,782 T→G 5.4% intergenic (+33/‑233) 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,970,7820TG 5.4% 876.2 / 11.0 293intergenic (+33/‑233)mtnA/PP_mr25methylthioribose‑1‑phosphate isomerase/gyrA regulatory region
Reads supporting (aligned to +/- strand):  ref base T (109/157);  new base G (17/4);  total (127/166)
Fisher's exact test for biased strand distribution p-value = 4.39e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

GTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTCCCCACCGCCCCCG  >  NC_002947/1970633‑1970931
                                                                                                                                                     |                                                                                                                                                     
 tctgcgcaatgaaggccgccgtgatgggcggtcttggaaactgccggatgttatcacaaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtttttgtattgcctgCACCGGCCTCTTCGCGGGCtt                                                                                                                                                      >  2:3157058‑M1/130‑150 (MQ=255)
     cGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  <  2:7183373/150‑1 (MQ=255)
            aCGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCAAGCCCGcgaaga                                                                                                                                           <  1:6169768/150‑6 (MQ=255)
             cGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTCGCCCGCCCCCAC                                                                                                                                          >  1:2943200‑M1/1‑150 (MQ=255)
              gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                                         <  1:1628154/150‑1 (MQ=255)
              gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                                         >  2:7026405/1‑150 (MQ=255)
              gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                                         >  2:7947073/1‑150 (MQ=255)
              gccggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACt                                                                                                                                         >  2:4820883/1‑150 (MQ=255)
                cggccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTgg                                                                                                                                       <  2:915780/150‑1 (MQ=255)
                  gccgACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGa                                                                                                                                     <  1:5380155/150‑1 (MQ=255)
                         tGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCa                                                                                                                              >  1:7467727/1‑150 (MQ=255)
                           aTCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATg                                                                                                                            <  1:505195/150‑1 (MQ=255)
                               aTGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCt                                                                                                                        >  2:1169825/1‑150 (MQ=255)
                                   gATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAgg                                                                                                                    <  2:7601140/150‑1 (MQ=255)
                                    aTCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAggg                                                                                                                   <  1:7329578/150‑1 (MQ=255)
                                        tGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgt                                                                                                               <  2:1949748/150‑1 (MQ=255)
                                         gACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCTCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgtg                                                                                                              <  2:2258400/150‑1 (MQ=255)
                                         gACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgtg                                                                                                              <  2:9130813/150‑1 (MQ=255)
                                         gACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGATTGCCCCCTCCCACAGGGATAACGCATGGACTTAGGGCCgtg                                                                                                              >  1:2211791/1‑150 (MQ=255)
                                           gccggatgttatcacaaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtttttgtattgcctgCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgtggt                                                                                                            >  1:660376‑M1/88‑150 (MQ=255)
                                           ccGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCgtggt                                                                                                            <  1:4937315/150‑1 (MQ=255)
                                            cGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCGCATGGCCTTAGGGCCgtggtg                                                                                                           >  1:1433075/1‑150 (MQ=255)
                                            cGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCTCCCCCAGGGATATCGCATGGCCTTAGGGCCggggtg                                                                                                           >  1:4737063/1‑150 (MQ=255)
                                               aCGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAAc                                                                                                        <  1:9021648/150‑1 (MQ=255)
                                                 gcgGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCt                                                                                                      <  1:8277895/150‑1 (MQ=255)
                                                   ggCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgt                                                                                                    <  1:1374203‑M2/150‑3 (MQ=255)
                                                    gCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACaggttcaccacggccctaaggccatgcgatatcccagtg                                                                                                   <  1:7235218‑M1/150‑40 (MQ=255)
                                                    gCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtg                                                                                                   >  2:2359357‑M2/1‑147 (MQ=255)
                                                     cATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgg                                                                                                  >  1:1039298‑M2/1‑146 (MQ=255)
                                                     cATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgg                                                                                                  >  1:4183055‑M2/1‑146 (MQ=255)
                                                     cATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgg                                                                                                  >  1:4056562‑M2/1‑146 (MQ=255)
                                                         gtcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCGTCGCGGGCGTGCCCGCTCCCACTGGGATATCGCATGGCCGTAGGGCCGTGGTGAACCGGGGGgggc                                                                                              >  2:5763640/1‑150 (MQ=255)
                                                         caaacgcgggggcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtttttgtattgcctgCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggc                                                                                              >  1:6821225‑M1/74‑150 (MQ=255)
                                                          tcgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcg                                                                                             <  2:6469915/150‑1 (MQ=255)
                                                           cgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcgg                                                                                            <  2:7867075/150‑1 (MQ=255)
                                                           cgAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTCGCCCGCCCCCACaggttcaccacggccctaaggccatgcgatatcccagtgggagcgg                                                                                            <  2:2943200‑M1/150‑47 (MQ=255)
                                                            gAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCGTAGGGCCGTGGTGTACCTGTGGgggcggg                                                                                           >  1:8875982/1‑150 (MQ=255)
                                                             aGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGgggcgggc                                                                                          <  1:38336/150‑1 (MQ=255)
                                                               cgcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGa                                                                                        <  2:1403773/150‑1 (MQ=255)
                                                               cgcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTCGCCCGCCCCCACaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaa                                                                                        <  2:8233284‑M1/150‑51 (MQ=255)
                                                               cgcCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGGATCGCCTCACCGGCCTCTTCGCGGGCTCGCCCGCCCCCACaggttcaccacggccctaaggccatgcgatatcccagtgggagcgggcaa                                                                                        <  2:8970373‑M1/150‑51 (MQ=255)
                                                                  ccGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGcc                                                                                     <  2:383417/150‑1 (MQ=255)
                                                                    gcggtggggagaggtggcagttggtctcagccgcgcttgaaacctggtttttgtattgcctgCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCcg                                                                                   >  2:4356675‑M1/63‑150 (MQ=255)
                                                                      aTGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCcgcg                                                                                 <  2:585418/150‑1 (MQ=255)
                                                                        gctgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTTTTTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTGGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGACCCTGTGGGGGCGGGCGAGCCCGCGaa                                                                               <  1:3209379/150‑1 (MQ=255)
                                                                        gctgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGaa                                                                               <  1:4324702/150‑1 (MQ=255)
                                                                          tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAaga                                                                             <  1:3513330/150‑1 (MQ=255)
                                                                          tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCCGGCATGACCGCAACCACAGGGATATCGCATGGCCATAGGGCCGGGGTGAACCTGTGGGGGCGGGCGGGCCCCCGaaaa                                                                             >  1:3563430/1‑148 (MQ=255)
                                                                          tgcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCGTCGCGGGCGTGCCCGCGCCCCCTGGGAGATCGCATGGCCGTAGGGCCGGGGTGAACCGGGGGGGGCGGGCGAGCCCGCGAagc                                                                             >  1:1414470/1‑149 (MQ=255)
                                                                           gcCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGGCCCCGCCCACTGGGGTATCGCCTGGCCGTGGGGGCGGGGTGAACCGGTGGGGGCGGGCGGGCCCGCGCagag                                                                            >  1:273098/1‑150 (MQ=255)
                                                                            ccAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCGTGGCGGGCGTGCCCGCGCCCACTGGGGTATCGCATGGCCGTAGGGGCGTGGGGAACCAGTGGGGGCGGGCGCGCCCGCGGAGCgg                                                                           >  1:1111128/1‑150 (MQ=255)
                                                                             cAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGc                                                                          <  1:3732449/150‑1 (MQ=255)
                                                                             cAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCCCTCCCACGGGGATATCGCATGGCCGTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGCGCCCGCGCAGAGGc                                                                          >  1:7070624/1‑150 (MQ=255)
                                                                               aGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCg                                                                        <  2:6097135/150‑1 (MQ=255)
                                                                               aGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCg                                                                        >  2:7564318/1‑150 (MQ=255)
                                                                                 cTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGccc                                                                                    >  1:5130719/1‑136 (MQ=255)
                                                                                 cTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGccc                                                                                    <  2:5130719/136‑1 (MQ=255)
                                                                                      tCAGCTGGTTTGCCGCTAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  <  1:2517149/69‑1 (MQ=255)
                                                                                      tCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  >  2:2517149/1‑69 (MQ=255)
                                                                                        aGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCa                                                               >  2:6776192/1‑150 (MQ=255)
                                                                                           tGGTTTGCCGCACGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATa                                                            >  2:4854683/1‑150 (MQ=255)
                                                                                              tttGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACaa                                                         >  2:3671453/1‑150 (MQ=255)
                                                                                              attGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCAGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCCGGCAATACaa                                                         >  2:8448652/2‑150 (MQ=255)
                                                                                                tGCCGCAAGCGCCTGCACTGATTTGGTCTCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATCTCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACaaaa                                                       <  2:561669/150‑1 (MQ=255)
                                                                                                 gCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACaaaaa                                                      <  2:7423194/150‑1 (MQ=255)
                                                                                                  ccGCACGCGCCTGCACTTATTTTGTACCGCCCCACCGGCCTGTGCGGGGGCTTGCCCGCCCCCCCGGGTATCTCGCAGGGCGTTCGGGCCGTGGTGAACCTGTGGGGGCGGGCGCGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAAc                                                     <  2:35569/150‑1 (MQ=255)
                                                                                                  ccGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCACCCACTGGGAAATCGCATGGCCTTTGGGCCGTGGTGTACCTGTGGGGGCGGGCGAGCCCCCGAAGAGGCCGGTGCAGGCACAACaaaaaa                                                     >  1:5877042/1‑149 (MQ=38)
                                                                                                    gCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCa                                                   <  2:8106952/150‑1 (MQ=255)
                                                                                                    gCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCTCCCCCTGGGATATCGCATGGCCGTAGGGCCGTGGGGTACCGGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACCAAAAACa                                                   >  1:3616747/1‑150 (MQ=38)
                                                                                                    gCAAGCGCCTGCACTGATTTTGTAGCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGGGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCa                                                   <  2:2720863/150‑1 (MQ=255)
                                                                                                        gcgcCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGGGTGCCCGCTCCCACCGGGGTATCGCATGGCCGTAGGGCCGTGGGGCACCGGGGGGGGCGGGCGAGCCCGCGAAGAGGCCGGGGCCGGCAATACACAAACCAGGtt                                               >  1:3371074/1‑150 (MQ=255)
                                                                                                          gcCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGGGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTc                                             >  2:4176051/1‑150 (MQ=255)
                                                                                                            cTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTCGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCaa                                           <  2:8713519/150‑1 (MQ=255)
                                                                                                            cTGCACTGATTTTGTATCGCCTCACCGGCCTCTCCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCaa                                           <  1:8151958/150‑1 (MQ=255)
                                                                                                              gCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACtggttcaccacggccctaaggccatgcgatatc                                                                                                         <  2:8307814‑M1/86‑34 (MQ=255)
                                                                                                                aCTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAgcgc                                       >  2:1610525/1‑150 (MQ=255)
                                                                                                                aCTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCGTGCCCGCGCCCACTGGGATATCGCATGGCCGTAGGGCCGTGGGGTACCGGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAgcgc                                       >  1:4473890/1‑150 (MQ=38)
                                                                                                                     ttttGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATTGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGGAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTg                                  >  2:67879/1‑150 (MQ=255)
                                                                                                                     ttttGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTg                                  >  2:4405472/1‑150 (MQ=255)
                                                                                                                        tgtcgcGCCTCACCGGCCTCTTCGGGGGCTGGCCCGCTCCCACGGGGATATCGCCGGGCGTTAGGGCCGGGGTGAACCTGGGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTgaga                               <  2:1753653/145‑1 (MQ=255)
                                                                                                                        tGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGTATATCGCTTGGCCTTCGGGCCGTGTTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTgaga                               <  2:325447/150‑1 (MQ=255)
                                                                                                                        tGTACCGCCTCCCCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTgaga                               <  2:944327/150‑1 (MQ=255)
                                                                                                                            tCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCaa                           >  2:5281226/1‑150 (MQ=255)
                                                                                                                             cGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAAc                          <  2:4737063/150‑1 (MQ=255)
                                                                                                                              gCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAAcg                         >  2:7285890/1‑149 (MQ=255)
                                                                                                                              gCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTGAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACt                         <  2:7696106/150‑1 (MQ=255)
                                                                                                                              gCCGCACCGGCCGCTGCGCGGGCGGGCCCGCCCCCCGGGGGATCTCGCAGGGCCTTCGGGCCGTGGTGAACCGGGGGGGGCGGGCGCGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACt                         <  1:2788064/150‑1 (MQ=255)
                                                                                                                                cTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATAGCATTGCCTTAGGGCCGTGGTGAACCAGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGc                       >  2:4017687/1‑150 (MQ=255)
                                                                                                                                  cccgggCCTCGTCGGGGGGTTCCCGGCACCCCCGGGTTTACCGCTTGGCCTTCGGGCGGGGGTGACCCGGGGGGGGGGGGCGCGCCCGCGACGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCa                     <  1:4644962/146‑1 (MQ=255)
                                                                                                                                   cCCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCAGGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCAc                    >  2:4233523‑M1/2‑150 (MQ=255)
                                                                                                                                       GCCTCATCGCGGGCTTGCCCGCTCCCACTGGGAGATCGCATGGCCTTAGGGCCGTGCTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCGGAGACCAACTGCCACctct                <  2:2995218‑M1/150‑1 (MQ=255)
                                                                                                                                          TCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTccc             <  1:5607222‑M1/150‑1 (MQ=255)
                                                                                                                                          TCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTccc             <  1:3287660‑M1/150‑1 (MQ=255)
                                                                                                                                          TCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTccc             <  2:1848616‑M1/150‑1 (MQ=255)
                                                                                                                                            TTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTCCCCa           <  2:1229208‑M1/150‑1 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgggggcgggcgagcccgcgaagaggccggtgaggcgatacaaaatcagtgcaggcgcttgcggcaaaccagctgagccagcttggcagcatccg      <  2:8012534‑M2/150‑97 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgggggcgggcgagcccgcgaagaggccggtgaggcgatacaaaatcagtgcaggcgcttgcggcaaaccagctgagccagcttggcagcatccg      <  2:8998904‑M2/150‑97 (MQ=255)
                                                                                                                                                 GGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTCCCCACCGcc      <  1:8470014‑M1/150‑1 (MQ=255)
                                                                                                                                                     ttGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTgtgggggcgggcgagcccgcgaagaggccggtgaggcgatacaaaatcagtgcaggcgcttgcggcaaaccagctgagccagcttggcagcatccgggcg  <  2:7691250‑M2/150‑101 (MQ=255)
                                                                                                                                                     TTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTCCCCACCGCCCCCg  <  1:694660‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGTTCGATGTAACGCCGGCCGACCTGATCGATGTGATCGTGACCGAACGCGGCATCGTCGAGCGCCCGGATGCTGCCAAGCTGGCTCAGCTGGTTTGCCGCAAGCGCCTGCACTGATTTTGTATCGCCTCACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACTGGGATATCGCATGGCCTTAGGGCCGTGGTGAACCTGTGGGGGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGCAATACAAAAACCAGGTTTCAAGCGCGGCTGAGACCAACTGCCACCTCTCCCCACCGCCCCCG  >  NC_002947/1970633‑1970931

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: