Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,306 A→T 5.0% intergenic (+91/‑74) pheA → / → PP_1770 bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,3060AT 5.0% 895.6 / 21.2 266intergenic (+91/‑74)pheA/PP_1770bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (161/91);  new base T (9/5);  total (170/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCG  >  NC_002947/1976157‑1976455
                                                                                                                                                     |                                                                                                                                                     
tCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTa                                                                                                                                                       <  1:5502846/150‑1 (MQ=255)
tCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTa                                                                                                                                                       <  1:4236117/150‑1 (MQ=255)
    ccAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGccc                                                                                                                                                   >  2:3223392/1‑150 (MQ=255)
      aggaggCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCcgc                                                                                                                                                 >  2:8116998/1‑150 (MQ=255)
      aggaggCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCcgc                                                                                                                                                 >  2:3503750/1‑150 (MQ=255)
           gCTGTGGCGCTGATGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGTGGCTGCCTTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAgg                                                                                                                                            <  1:6423902/150‑1 (MQ=255)
              gtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgct                                                                                                                                         <  1:5771580/150‑2 (MQ=255)
                ggCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgc                                                                                                                                       >  2:2665825/1‑150 (MQ=255)
                          gTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGccc                                                                                                                             <  1:8214397/150‑1 (MQ=255)
                             atGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAGCATGAGATGCCGATATCTGTGGGAGCGGGCTACCCAGCGAAGGGCCGCGCCGCCGCCCCaa                                                                                                                          >  1:1067538/2‑150 (MQ=39)
                                 gCTCTTATCCGACGGCGGTGCTTTGATTCCGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGGGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                                                      <  1:4815337/150‑1 (MQ=255)
                                  ctctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                                                                                                     >  2:2566192/1‑150 (MQ=255)
                                   tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGccc                                                                                                                    >  2:7113903/1‑150 (MQ=255)
                                     ttATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTc                                                                                                                  >  1:6046048/1‑150 (MQ=255)
                                       aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                                                                <  2:605743/150‑1 (MQ=255)
                                       aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                                                                <  2:1199666/150‑1 (MQ=255)
                                       aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                                                                >  1:8361078/1‑150 (MQ=255)
                                        tCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcaca                                                                                                               <  2:8265119/150‑1 (MQ=255)
                                         ccGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAg                                                                                                              <  1:6599069/150‑1 (MQ=255)
                                         ccGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAg                                                                                                              <  2:8767081/150‑1 (MQ=255)
                                         ccGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAg                                                                                                              <  1:9031888/150‑1 (MQ=255)
                                           gAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGc                                                                                                            >  1:7220727/1‑150 (MQ=255)
                                               gCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGTGAGCATGTGAGGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGACGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAcc                                                                                                        <  1:6342792/150‑1 (MQ=255)
                                                  gTGCTTTGATTCCGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAg                                                                                                     <  1:3618183/150‑1 (MQ=255)
                                                  gTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGATCAAGTGATGACGATCTCTGTGGGAGCGGGCAAGCCCCCGAAGGGCCCCGCCGCCGCCCCAAAGGCCCACAAAGGCTACCAAg                                                                                                     >  1:2523633/1‑150 (MQ=37)
                                                      tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                 >  2:6811050/1‑150 (MQ=255)
                                                      tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                 >  2:61996/1‑150 (MQ=255)
                                                      tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                 >  2:6989623/1‑150 (MQ=255)
                                                        tGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                                                                                               <  2:6356854/150‑1 (MQ=255)
                                                        tGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                                                                                               <  1:7015581/150‑1 (MQ=255)
                                                         gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGTGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTg                                                                                              <  2:7626512/150‑1 (MQ=255)
                                                         gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTg                                                                                              <  2:3002326/150‑1 (MQ=255)
                                                         gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCCCGCCGCGGCCCCACTGGCCCGCACAGGCTACCGAGATTTCTg                                                                                              >  1:1173811/1‑150 (MQ=255)
                                                           ttCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGaa                                                                                            >  2:4360742/1‑150 (MQ=255)
                                                           ttCAGGGGCCTTGGGGCTGCCTGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGCGCTTGTGTTGCCGATCTCTGTGGGAGCGGGCTTGCCCGCGATGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGaa                                                                                            <  2:7551135/150‑1 (MQ=255)
                                                             ccggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATTTGATGCCTATCTCTGTGGGTGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACa                                                                                          <  1:3463143/148‑1 (MQ=255)
                                                               ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                        <  2:6923364/150‑1 (MQ=255)
                                                               ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                        <  2:6251271/150‑1 (MQ=255)
                                                               ggggCATTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                        <  1:8584781/150‑1 (MQ=255)
                                                                  gCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGtt                                                                                     <  2:228291/150‑1 (MQ=255)
                                                                     ttGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGc                                                                                  <  1:2718679/150‑1 (MQ=255)
                                                                     ttGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGGTCGc                                                                                  >  1:6780212/1‑150 (MQ=255)
                                                                      aagggCTGCAAGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCa                                                                                 >  2:1852254/3‑150 (MQ=255)
                                                                       agggCTGCAAGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCAc                                                                                >  2:2004661/2‑150 (MQ=255)
                                                                              ccTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGt                                                                         <  2:4116865/150‑1 (MQ=255)
                                                                                ttGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCCATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTaa                                                                       >  2:8287437/1‑150 (MQ=255)
                                                                                 tGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCCATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTaaa                                                                      >  2:5620611/1‑150 (MQ=255)
                                                                                    aGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCCCCGGCTACCTAGATTTCTGAACAGGGTTCGCACCCGGGGTAAAAgc                                                                   >  1:3985552/1‑150 (MQ=255)
                                                                                    aGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgc                                                                   >  2:6556717/1‑150 (MQ=255)
                                                                                     gcccttcgcgggctagcccgctcccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCATGCCCGCAAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgca                                                                  >  2:513066‑M2/47‑150 (MQ=255)
                                                                                     gCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgca                                                                  <  2:8607841/150‑1 (MQ=255)
                                                                                     gCCCTTCGCGGGCAGGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgca                                                                  <  1:5627914/150‑1 (MQ=255)
                                                                                      cccTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCGCACAGGCGACCGAGATTTCTGAACAGGGTTCGCACCCCTGGTAAATgcag                                                                 >  1:1614750/1‑150 (MQ=255)
                                                                                       ccttcgcgggctagcccgctcccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCATGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgcagc                                                                >  2:1618023‑M2/45‑150 (MQ=255)
                                                                                       ccTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGACGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgcagc                                                                >  2:1397285/1‑150 (MQ=255)
                                                                                            gcgGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCaaaaa                                                           <  1:2387087/150‑1 (MQ=255)
                                                                                                  aTGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAaaaccaaac                                                     <  1:8822160/150‑1 (MQ=255)
                                                                                                    gCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGtt                                                                                     >  1:5226192/1‑116 (MQ=255)
                                                                                                      ccgctcccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgc                                                 >  1:7129527‑M2/30‑150 (MQ=255)
                                                                                                      ccGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCGAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAAACAAACCcgc                                                 >  1:4066921/1‑150 (MQ=255)
                                                                                                      ccGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCGCACAGGCTACCGAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgc                                                 >  1:5043815/1‑150 (MQ=255)
                                                                                                      ccGCTCACACAGTGAGCATTTGATGACGATATCTGTGGGATCTGGCTAGACCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGCGTTCGCACCAGAGGAAAAAGCAGCAAAAACCAAAACcgc                                                 >  2:7444665/1‑150 (MQ=255)
                                                                                                          tCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGccc                                             >  2:8995892/1‑150 (MQ=255)
                                                                                                          tCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAAACCGCGccc                                             >  1:923854/1‑150 (MQ=255)
                                                                                                            ccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCATGCCCGCGAAGGGCTGCAAGGCAGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcat                                           >  2:7573007‑M2/24‑150 (MQ=255)
                                                                                                            ccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCATGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTc                                                                                    <  2:3974313‑M2/86‑1 (MQ=255)
                                                                                                            ccacagagatcggcatcacatgcTCACTGTGGGAGCGGGCATGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTc                                                                                    >  1:3974313‑M2/24‑109 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcat                                           >  2:4602024/1‑150 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcat                                           >  2:6859729/1‑150 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcag                                           >  1:8938549/1‑149 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcaa                                           >  1:1921071/1‑149 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGGGGTAAATGCAGCAAAAACCAAACCCGCGCCcat                                           >  1:1389990/1‑150 (MQ=255)
                                                                                                            ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAACTAGATTTGTTTACATGGTTCTCACCAGTGGTTAATGCATCAAAAAAAAAAAACGCGCCcat                                           >  1:48797/1‑150 (MQ=255)
                                                                                                             cacagagatcggcatcacatgcTCACTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatc                                          >  1:864648‑M2/23‑150 (MQ=255)
                                                                                                             cacaGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatc                                          >  2:4933977/1‑150 (MQ=255)
                                                                                                               caGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAGTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatcag                                        >  2:5366800/1‑149 (MQ=255)
                                                                                                                  tGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGa                                     <  1:1879440/150‑1 (MQ=255)
                                                                                                                    aGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGAcc                                   >  2:2750951/1‑150 (MQ=255)
                                                                                                                        tgtgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCt                               <  2:8361078/150‑1 (MQ=255)
                                                                                                                         gtgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCaaaaa                                                           <  1:1250566/121‑1 (MQ=255)
                                                                                                                         gtgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtg                              >  2:5960584/1‑150 (MQ=255)
                                                                                                                           gATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggt                            <  1:632445/150‑1 (MQ=255)
                                                                                                                           gATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggt                            <  2:4859267/150‑1 (MQ=255)
                                                                                                                            ctGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtc                           <  2:2127926/149‑1 (MQ=255)
                                                                                                                              gCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtcgt                         >  2:552140/1‑150 (MQ=255)
                                                                                                                                 gATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtcgtggt                      <  2:238873/150‑1 (MQ=255)
                                                                                                                                     tctGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTggtcggt                  <  1:7632475/150‑1 (MQ=255)
                                                                                                                                      ctGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTggtcggtc                 >  2:3718130/1‑150 (MQ=255)
                                                                                                                                      ctGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTggtcggtc                 >  2:387242/1‑150 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                        >  2:2087558/1‑78 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                        <  1:2086571/78‑1 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCt                >  2:867113/1‑150 (MQ=255)
                                                                                                                                         tgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATTCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAg              >  1:1536665/1‑150 (MQ=255)
                                                                                                                                         tgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAg              >  2:4520865/1‑150 (MQ=255)
                                                                                                                                         tgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAg              <  1:2151814/150‑1 (MQ=255)
                                                                                                                                         tgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAg              <  1:3558739/150‑1 (MQ=255)
                                                                                                                                           ttcgCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGc            >  1:6502601/4‑150 (MQ=255)
                                                                                                                                            tcgCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGcc           >  1:4953964/3‑150 (MQ=255)
                                                                                                                                            gAGCGGGCATGCCCGCGAAGGGCTGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGcc           >  1:8743800/1‑150 (MQ=255)
                                                                                                                                                 ggCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATc      <  2:52582/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCG  >  NC_002947/1976157‑1976455

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: