Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,336 A→T 6.8% intergenic (+58/‑99) ppsR → / → aroH PEP synthase kinase and PEP synthase pyrophosphorylase/3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,3360AT 6.8% 236.1 / 6.7 89intergenic (+58/‑99)ppsR/aroHPEP synthase kinase and PEP synthase pyrophosphorylase/3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base T (3/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

Only 100 of 103 total aligned reads displayed.

GATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACG  >  NC_000913/1788188‑1788485
                                                                                                                                                    |                                                                                                                                                     
taTCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      <  2:1726710/149‑1 (MQ=255)
gATCCCGTTGATTAACAGTACCAATTAATCGGTAGAAGAGATTGCCAACAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGTACTAGTGCATTAGCTTATTTTTTTGGTATCATGCTAACCACCCGGCGAGGGGTGac                                                                                                                                                      <  1:817564/150‑1 (MQ=255)
gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGCGCCTCGATATCATGGGCGTTAGTCGCCGAATGTACTAGAGAACTAGGGCAGTAGCTTATTTTTGTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      >  1:857944/1‑150 (MQ=255)
gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      <  1:763904/150‑1 (MQ=255)
gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      >  2:1626630/1‑150 (MQ=255)
gATACCGTGTATTAAGAGTAACAATTATTCGGTAGAATAGATTGCCACCTAGATCCTCTTTATCAAGCGACTTAGACGCCGAATGTACTAGAGTTTTAGTGAATAAGTTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      <  1:288893/150‑1 (MQ=255)
  tCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacac                                                                                                                                                    <  1:1365860/150‑1 (MQ=255)
   cccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacaca                                                                                                                                                   <  2:804740/150‑1 (MQ=255)
    ccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacacac                                                                                                                                                  >  2:279206/1‑150 (MQ=255)
    ccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGACACCAAGATCCTCGATAGCATGGGCCTTAGTCGCAGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCCCCCGGCGAGGTGTGaaaaac                                                                                                                                                  >  2:1118545/1‑145 (MQ=255)
     cGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACAcc                                                                                                                                                 >  2:1644922/1‑150 (MQ=255)
      gTGGATAAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCt                                                                                                                                                <  2:864810/150‑1 (MQ=255)
        ggATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGtca                                                                                                                                                     >  2:184587/1‑143 (MQ=255)
        ggATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGtca                                                                                                                                                     <  1:184587/143‑1 (MQ=255)
        ggATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACTCACCTCg                                                                                                                                              >  2:250766/1‑150 (MQ=255)
         gATTAACAGTAGCAATTATTCGGTAGAAGAGATTGCCACCAAGTTCCTCTATATCTTGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATTCTAACCACCCGGCGAGGTGTGTCACACCgcgc                                                                                                                                             <  1:278092/150‑1 (MQ=255)
         gATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGc                                                                                                                                             <  2:1533827/150‑1 (MQ=255)
          aTTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGAGTGACCCACCACGCa                                                                                                                                            >  2:1570723/1‑150 (MQ=255)
             aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccggg                                                                                                                                         <  2:125501/150‑5 (MQ=255)
              aCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTg                                                                                                                                        <  1:1487046/150‑1 (MQ=255)
                 gTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaaa                                                                                                                                     >  1:1546108/1‑150 (MQ=255)
                 gTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATAACATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGAAACACATCGCACTTGaaa                                                                                                                                     >  1:628326/1‑150 (MQ=255)
                  tACCAATTAATCGGTAGAAGAGATTTCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTACTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAAt                                                                                                                                    <  1:55834/150‑1 (MQ=255)
                   aCCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATc                                                                                                                                   <  1:1644111/150‑1 (MQ=255)
                    ccAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                                  >  1:1620899/1‑150 (MQ=255)
                      aaTTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagc                                                                                                                                >  2:700722/1‑150 (MQ=255)
                          attCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCg                                                                                                                            >  1:1412792/1‑150 (MQ=255)
                          attCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCg                                                                                                                            <  2:1295120/150‑1 (MQ=255)
                                 agaagaGATTGCCACCAAGATCCTCGATTACATGGGCCTTTGTCGCCGTATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGtt                                                                                                                     <  2:628326/150‑1 (MQ=255)
                                 agaagaGATTGCCACCAAGATACTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGtt                                                                                                                     <  1:659230/150‑1 (MQ=255)
                                  taagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGttt                                                                                                                    <  1:376028/149‑1 (MQ=255)
                                       gaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGt                                                                                                               >  2:1398541/1‑150 (MQ=255)
                                        aTTGCCACCAAGATCCTCGATATCATGGGCCCTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTg                                                                                                              <  2:959515/150‑1 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGCGGTTTGACACCCCACGAACTTGAAATCAGCAGCGAATTGTTTATCGTGATg                                                                                                           >  1:1424340/1‑150 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                                                                           >  2:186454/1‑150 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                                                                           >  2:6290/1‑150 (MQ=255)
                                             caccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGCACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATgcg                                                                                                         >  1:848792/1‑150 (MQ=255)
                                             caccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGAGATTGGTTTATCGTGATgcg                                                                                                         >  2:1344238/1‑150 (MQ=255)
                                               ccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                                                                       >  2:717758/1‑150 (MQ=255)
                                                    aTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACt                                                                                                  <  1:1404840/150‑1 (MQ=255)
                                                               cATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTc                                                                                       >  2:295540/1‑150 (MQ=255)
                                                                 tGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGCCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGCTGCGCATCACTTCCCGGCAGTccg                                                                                     >  1:1663406/1‑149 (MQ=255)
                                                                 tGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCt                                                                                     <  2:43726/150‑1 (MQ=255)
                                                                    gCCTTAGTCGCCGAATGTACTAGTGTAGTAGTGCTTTTGCTTTTTTTTTTGTTATCATGCAAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGcc                                                                                  <  1:770310/150‑1 (MQ=255)
                                                                     ccTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCg                                                                                 <  2:608748/150‑1 (MQ=255)
                                                                       ttAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTa                                                                               <  2:164715/150‑1 (MQ=255)
                                                                          gTCGCCGAATGTACTAGAGAAATAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGAGAGGTGTCACACACCTAGCACTTGAAATCAGCATCTATTGGTTTATCTTGAATCTCATAACTTAACGGCAGTCCTGCCGTACaa                                                                            >  1:1053299/1‑150 (MQ=255)
                                                                          gTCGCCGAAAGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGaa                                                                            <  1:1232294/150‑1 (MQ=255)
                                                                             gCCGATTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATCCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGca                                                                         <  2:1333434/150‑1 (MQ=255)
                                                                               cGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGcaac                                                                       <  1:985704/150‑1 (MQ=255)
                                                                                 aaTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGcaacaa                                                                     <  1:279206/150‑1 (MQ=255)
                                                                                      aCTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTc                                                                <  2:1620899/150‑1 (MQ=255)
                                                                                      aCTAGAGAACTAGTGCATAAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTc                                                                <  1:891923/150‑1 (MQ=255)
                                                                                       cTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATTCGCATCAATTCCCGGCAGTCCAGCCGTATAAGAAACAAATTTCt                                                               >  1:1423047/1‑150 (MQ=255)
                                                                                       cTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt                                                               >  2:494077/1‑150 (MQ=255)
                                                                                       cTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt                                                               >  2:1372000/1‑150 (MQ=255)
                                                                                        tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTg                                                              <  1:917953/150‑1 (MQ=255)
                                                                                        tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGGGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTg                                                              <  1:1733543/150‑1 (MQ=255)
                                                                                        tAGAGAACTAGTGCATTAGATTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACATCGAAATTGAAATAAGCAGCTATTGTTTTATCTTGATGCGAATCACTTCCAGGAAGTCATGCCGTAGAAGCAACAAAtttatt                                                              >  1:30337/1‑147 (MQ=255)
                                                                                         agagAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                             >  1:439498/1‑150 (MQ=255)
                                                                                            gtaCTAGTGCATTAGATTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGGGTGACACACCTCGCACTTGAAATCTGCTGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          <  2:1377406/148‑1 (MQ=255)
                                                                                             aaCTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACt                                                         <  2:726601/150‑1 (MQ=255)
                                                                                                tAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGt                                                      <  2:1412792/150‑1 (MQ=255)
                                                                                                 aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTa                                                     >  2:700395/1‑150 (MQ=255)
                                                                                                  gTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATAAGAAGCGATTGGATTATCGTGATGCGCATCACTTACCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTaa                                                    >  1:62032/1‑150 (MQ=255)
                                                                                                     cATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGa                                                 >  1:892048/1‑150 (MQ=255)
                                                                                                         aGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCCCTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATAACTTCCCGGCAGTCCTTGCGTAGAAGGAACAAATTTATGAGACTTTTTATGAAAAg                                             >  1:1313730/1‑150 (MQ=255)
                                                                                                           cTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACTCACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCCTCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGaa                                           >  1:323686/1‑150 (MQ=255)
                                                                                                           cTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGaa                                           >  2:931411/1‑150 (MQ=255)
                                                                                                           cTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCGGCCGTAGAAGCAACAAATTGATGAGACTTGTAATGAACAGac                                           >  1:328615/1‑149 (MQ=255)
                                                                                                             tATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACt                                         >  2:1358880/1‑150 (MQ=255)
                                                                                                              aTTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTg                                        <  2:1641967/150‑1 (MQ=255)
                                                                                                                       ttATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTcc                               <  1:1177528/150‑1 (MQ=255)
                                                                                                                         aTCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGt                             >  1:662642/1‑150 (MQ=255)
                                                                                                                           cATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTAc                           <  2:1608383/150‑1 (MQ=255)
                                                                                                                             tGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGACCTGCCGAAGAAGCAACAAATTTCTGAGACTTGTAAAGAACAGAACTGACGAACTCCGTAcag                         >  1:149203/1‑148 (MQ=255)
                                                                                                                              gCTAACCCCCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTgc                        <  1:893390/150‑1 (MQ=255)
                                                                                                                                 aaCCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGcgcg                     <  1:41225/150‑1 (MQ=255)
                                                                                                                                    caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAt                  >  2:1679137/1‑150 (MQ=255)
                                                                                                                                    caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAt                  >  2:582710/1‑150 (MQ=255)
                                                                                                                                     accCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTATATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAtt                 >  1:386871/1‑150 (MQ=255)
                                                                                                                                      cccGGCGCTGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg                <  1:126587/150‑1 (MQ=255)
                                                                                                                                      cccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg                <  1:295540/150‑1 (MQ=255)
                                                                                                                                      cccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg                <  2:1401637/150‑1 (MQ=255)
                                                                                                                                      cccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg                <  1:887897/150‑1 (MQ=255)
                                                                                                                                       agtgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTga               >  1:78799/4‑150 (MQ=255)
                                                                                                                                        gtgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag              >  2:114115/3‑150 (MQ=255)
                                                                                                                                        cGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag              <  1:1455736/150‑1 (MQ=255)
                                                                                                                                        cGGCGAGGGGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTAGCGTGATGCGCATCACTTCCCGGCAGGCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTACTGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag              <  1:1118545/150‑1 (MQ=255)
                                                                                                                                         ggCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGCACTGACGAACTCCGTACTGCGCGTATTgaga             >  2:79580/1‑150 (MQ=255)
                                                                                                                                          gCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTCCTGCGCGTATTGagag            >  1:1601455/1‑150 (MQ=255)
                                                                                                                                           cGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGc           <  1:1398541/150‑1 (MQ=255)
                                                                                                                                             aGGTGTGACACACCTCGCACTTGAAATCAGAAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTACTGCCGTAGCAGCAACAAATTTCTGAGACTTGTTATGAACATAACTGACGAACTCCGTACTGCGAGGATTGAGAGACt         >  1:564980/1‑150 (MQ=255)
                                                                                                                                               gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTgg       <  1:1372000/150‑1 (MQ=255)
                                                                                                                                               gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTgg       >  2:71535/1‑150 (MQ=255)
                                                                                                                                               gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTgg       >  2:639138/1‑150 (MQ=255)
                                                                                                                                               gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTgg       <  2:848792/150‑1 (MQ=255)
                                                                                                                                                 gtgACACACCTCGCACTTGAAATTAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTa     >  1:998016/1‑150 (MQ=255)
                                                                                                                                                   gACACACCTCGCAGTTGAAATCCGCATCGATTGGGTTATCGTGATGCGCATCACTTCCCGGCATTCCTGCAGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAAc   <  1:1570723/150‑1 (MQ=255)
                                                                                                                                                    acacacCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAAcg  <  1:717758/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACG  >  NC_000913/1788188‑1788485

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