Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,400,868 C→G 5.1% intergenic (+123/‑119) aroQ‑III → / → PP_3004 type II 3‑dehydroquinate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,400,8680CG 5.1% 532.0 / 13.4 158intergenic (+123/‑119)aroQ‑III/PP_3004type II 3‑dehydroquinate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (72/78);  new base G (3/5);  total (75/83)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 159 total aligned reads displayed.

GCGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAACAAG  >  NC_002947/3400721‑3401017
                                                                                                                                                   |                                                                                                                                                     
 cgtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatg                                                                                                                                                    <  2:4538661/150‑1 (MQ=255)
  gtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatgg                                                                                                                                                   <  2:4694179/150‑1 (MQ=255)
  gtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatgg                                                                                                                                                   >  1:1191359/1‑150 (MQ=255)
  gtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatgg                                                                                                                                                   <  1:5111061/150‑1 (MQ=255)
   tgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatgga                                                                                                                                                  <  1:3639603/150‑1 (MQ=255)
    gcggcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatggat                                                                                                                                                 <  2:488478/150‑1 (MQ=255)
        cgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaa                                                                                                                                             <  1:1842649/150‑1 (MQ=255)
         gATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                            <  2:5081700/150‑4 (MQ=255)
         gATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaaa                                                                                                                                            <  1:100726/150‑1 (MQ=255)
               cATGGTCTGAGCAGGTCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGt                                                                                                                                      >  2:3188278/1‑150 (MQ=255)
                aTGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTa                                                                                                                                     <  1:1760951/150‑1 (MQ=255)
                aTGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTa                                                                                                                                     <  2:4667249/150‑1 (MQ=255)
                   gTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCg                                                                                                                                  >  1:5006896/1‑150 (MQ=255)
                   gTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCg                                                                                                                                  <  2:1329267/150‑1 (MQ=255)
                   gTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCg                                                                                                                                  >  1:2564663/1‑150 (MQ=255)
                   gTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCg                                                                                                                                  >  1:2564690/1‑150 (MQ=255)
                      tGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACt                                                                                                                               >  2:2790484/1‑150 (MQ=255)
                        aGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTc                                                                                                                             >  1:5066420/1‑150 (MQ=255)
                         gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTcc                                                                                                                            <  1:4130629/150‑1 (MQ=255)
                         gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTcc                                                                                                                            <  1:2470730/150‑1 (MQ=255)
                         gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTcc                                                                                                                            <  2:885004/150‑1 (MQ=255)
                           aGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGt                                                                                                                          <  1:351593/150‑1 (MQ=255)
                           aGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGt                                                                                                                          <  1:7534/150‑1 (MQ=255)
                                  aGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCa                                                                                                                   <  2:3482795/150‑1 (MQ=255)
                                    gggTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATa                                                                                                                 <  1:4985050/150‑1 (MQ=255)
                                      gTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                            >  2:3722697/1‑118 (MQ=255)
                                      gTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                            <  1:3722697/121‑4 (MQ=255)
                                      gTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGt                                                                                                               >  1:3798042/1‑150 (MQ=255)
                                      gTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGt                                                                                                               >  1:3798052/1‑150 (MQ=255)
                                       tCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTa                                                                                                              <  1:4852663/150‑1 (MQ=255)
                                         ggATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGTGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACg                                                                                                            >  1:738066/1‑150 (MQ=255)
                                         ggATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACg                                                                                                            >  1:3026194/1‑150 (MQ=255)
                                          gATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACgg                                                                                                           >  1:2518982/1‑150 (MQ=255)
                                             ctcCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTcc                                                                                                        >  1:1970958/1‑150 (MQ=255)
                                                cGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGcc                                                                                                     >  1:4654247/1‑150 (MQ=255)
                                                cGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGcc                                                                                                     >  1:1178509/1‑150 (MQ=255)
                                                      cgTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaac                                                                                               <  2:4193857/150‑1 (MQ=255)
                                                       gTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaaca                                                                                              >  1:1430957/1‑150 (MQ=255)
                                                        tCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaacaa                                                                                             <  2:412188/150‑1 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  2:386911/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  2:20041/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  1:4585787/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  2:386950/1‑150 (MQ=255)
                                                                cgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCg                                                                                     <  1:1220912/150‑1 (MQ=255)
                                                                   aGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGtt                                                                                  >  2:2706194/1‑150 (MQ=255)
                                                                    gCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTc                                                                                 <  1:802324/150‑1 (MQ=255)
                                                                     cAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCg                                                                                >  1:4559635/1‑150 (MQ=255)
                                                                      aGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGt                                                                               <  2:4939788/150‑1 (MQ=255)
                                                                      aGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGt                                                                               <  1:2756096/150‑1 (MQ=255)
                                                                       ggTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGtt                                                                              <  2:68813/150‑1 (MQ=255)
                                                                       ggTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGtt                                                                              >  1:3290278/1‑150 (MQ=255)
                                                                                 ggCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAAcc                                                                    >  2:2878177/1‑150 (MQ=255)
                                                                                      gcggcgACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAg                                                                                                                                    >  2:279058/1‑81 (MQ=255)
                                                                                           gACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCGGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCg                                                          <  2:790273/150‑1 (MQ=255)
                                                                                             cGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCggg                                                        <  2:2743966/150‑1 (MQ=255)
                                                                                              ggCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCgggg                                                       >  2:4384798/1‑150 (MQ=255)
                                                                                                cgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGct                                                     <  2:2478372/150‑1 (MQ=255)
                                                                                                     gAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCAt                                                >  1:1310898/1‑150 (MQ=255)
                                                                                                     gAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCAt                                                >  1:459943/1‑150 (MQ=255)
                                                                                                         tGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGAcc                                            <  2:4674325/150‑1 (MQ=255)
                                                                                                            ccGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg                                         >  1:495015/1‑150 (MQ=255)
                                                                                                            ccGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg                                         <  2:459943/150‑1 (MQ=255)
                                                                                                             cGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTggg                                        <  1:1067004/150‑1 (MQ=255)
                                                                                                             cGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTggg                                        >  2:1690766/1‑150 (MQ=255)
                                                                                                              gTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                       >  2:810327/1‑150 (MQ=255)
                                                                                                              gTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                       <  2:4592703/150‑1 (MQ=255)
                                                                                                              gTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                       <  1:1745733/150‑1 (MQ=255)
                                                                                                                aGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAgg                                     >  2:1470851/1‑150 (MQ=255)
                                                                                                                aGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAgg                                     >  2:1470877/1‑150 (MQ=255)
                                                                                                                    gtggtgGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGa                                 >  2:4057235/1‑150 (MQ=255)
                                                                                                                     tggtggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGaa                                <  2:5011103/150‑1 (MQ=255)
                                                                                                                     tggtggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGaa                                <  1:2252338/150‑1 (MQ=255)
                                                                                                                      ggtggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTcatga                               <  2:495015/150‑5 (MQ=255)
                                                                                                                      ggtggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAAt                               <  1:3530650/150‑1 (MQ=255)
                                                                                                                       gtggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATg                              >  1:1640146/1‑150 (MQ=255)
                                                                                                                          gATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                           <  2:2986726/150‑1 (MQ=255)
                                                                                                                          gATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                           <  2:2059626/150‑1 (MQ=255)
                                                                                                                            tGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGt                         <  2:1191410/150‑1 (MQ=255)
                                                                                                                            tGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGt                         <  2:1191359/150‑1 (MQ=255)
                                                                                                                              tCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCg                       <  1:3622187/150‑1 (MQ=255)
                                                                                                                                     ggctttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAAc                >  1:2039789/7‑150 (MQ=255)
                                                                                                                                     gcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAAc                >  1:2941051/1‑150 (MQ=255)
                                                                                                                                        tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  2:2074163/4‑150 (MQ=255)
                                                                                                                                        tttccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  1:1580003/4‑150 (MQ=255)
                                                                                                                                        tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  2:1238717/1‑150 (MQ=255)
                                                                                                                                        tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  2:2871525/1‑150 (MQ=255)
                                                                                                                                         aCCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATc            >  2:472500/1‑150 (MQ=255)
                                                                                                                                         aCCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATc            >  1:4397694/1‑150 (MQ=255)
                                                                                                                                          cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCa           <  2:4654247/150‑1 (MQ=255)
                                                                                                                                          cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCa           <  1:2790484/150‑1 (MQ=255)
                                                                                                                                          cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCa           <  1:4143643/150‑1 (MQ=255)
                                                                                                                                           ccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAg          >  2:4992929/1‑150 (MQ=255)
                                                                                                                                           ccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAg          >  2:2080099/1‑150 (MQ=255)
                                                                                                                                             aTCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTc        >  2:657475/1‑150 (MQ=255)
                                                                                                                                             aTCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTc        <  2:3026194/150‑1 (MQ=255)
                                                                                                                                              tCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTca       >  2:2871033/1‑150 (MQ=255)
                                                                                                                                                 aGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaaca    >  2:2247453/1‑150 (MQ=255)
                                                                                                                                                 aGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaaca    >  2:2223252/1‑150 (MQ=255)
                                                                                                                                                  gCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGTCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaacaa   <  2:1430957/150‑1 (MQ=255)
                                                                                                                                                   cATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAACAAg  >  1:4431792/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAACAAG  >  NC_002947/3400721‑3401017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: