Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,336 A→T 16.2% 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,3360AT16.2% 235.4 / 36.5 98intergenic (+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 (6/10);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 103 total aligned reads displayed.

AGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAA  >  NC_000913/1788238‑1788435
                                                                                                  |                                                                                                   
aGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                    <  1:2003034/100‑1 (MQ=255)
aGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                    >  2:763159/1‑100 (MQ=255)
aGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                    <  1:1416129/100‑1 (MQ=255)
  aTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacac                                                                                                  >  2:108681/1‑100 (MQ=255)
      tCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCt                                                                                              <  1:1802749/100‑1 (MQ=255)
         atatCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGc                                                                                           <  1:2653368/100‑1 (MQ=255)
          tatCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCa                                                                                          >  2:2517510/1‑100 (MQ=255)
          tatCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACAACCCGGCGAGGTGTGACAAACCTCGCa                                                                                          >  1:1252696/1‑100 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTTGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                         <  1:454577/100‑3 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                         <  2:48398/100‑3 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                         <  2:2243976/100‑3 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                         <  2:1441141/100‑3 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCAAGGTGTGTCACACCTCGccg                                                                                         <  1:1652814/100‑3 (MQ=255)
           atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTAGTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                         <  2:713169/100‑3 (MQ=255)
            tCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgg                                                                                        <  2:2335755/100‑4 (MQ=255)
             cATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccggg                                                                                       <  1:1954619/100‑5 (MQ=255)
               tGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgggtg                                                                                     <  2:1236095/100‑7 (MQ=255)
               tGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGa                                                                                     <  2:842372/100‑1 (MQ=255)
                gggCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaa                                                                                    <  1:571293/100‑1 (MQ=255)
                 ggCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTTTTTTTTTTCTAACCACCCAGCGAGTTGTGTTACACCTCGCACTTGtaa                                                                                   >  1:2503301/1‑100 (MQ=255)
                  gCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgggtggtt                                                                                  <  2:933868/100‑10 (MQ=38)
                   ccTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATc                                                                                 >  2:2190616/1‑100 (MQ=255)
                     ttAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcag                                                                               <  1:1752047/100‑1 (MQ=255)
                      tAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagc                                                                              <  2:2383788/100‑1 (MQ=255)
                        gTCGCCGAATGTACTAGAGATCTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCAACCGGCGAGGTGTGACACACCTCGCTCTTGAAATcagca                                                                             <  1:1796368/99‑1 (MQ=255)
                         tctccGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagc                                                                           <  1:2801235/97‑1 (MQ=255)
                         tctccGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagc                                                                           <  1:2801715/97‑1 (MQ=255)
                          cGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCg                                                                          <  1:1133587/100‑1 (MQ=255)
                          cGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCg                                                                          >  2:675298/1‑100 (MQ=255)
                           gCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcggc                                                                           >  2:1288405/1‑98 (MQ=255)
                            ccGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagcaga                                                                        >  2:2801245/1‑97 (MQ=255)
                            ccGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagcaga                                                                        >  2:2801725/1‑97 (MQ=255)
                             cGAATGTACTAGAGAACTAGTGCATTAACTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagc                                                                           >  1:798315/1‑96 (MQ=255)
                               aaTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTgg                                                                     <  2:352003/100‑1 (MQ=255)
                                aTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGt                                                                    <  1:1876241/100‑1 (MQ=255)
                                aTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGt                                                                    <  2:834443/100‑1 (MQ=255)
                                aTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTCATCATGCTAACCACCCGGCTAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGt                                                                    <  1:1582808/100‑1 (MQ=255)
                                      tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCg                                                              >  2:46870/1‑100 (MQ=255)
                                      tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCg                                                              <  2:176227/100‑1 (MQ=255)
                                      tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCg                                                              <  2:2853960/100‑1 (MQ=255)
                                       agagAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCg                                                              >  2:1834002/1‑99 (MQ=255)
                                       agagAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGt                                                             >  2:2170322/1‑100 (MQ=255)
                                        gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTg                                                            >  2:1492194/1‑100 (MQ=255)
                                        gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTg                                                            >  2:2554474/1‑100 (MQ=255)
                                           aaCTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACATACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                         <  1:1592209/100‑1 (MQ=255)
                                             cTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATgcg                                                       >  2:1164292/1‑100 (MQ=255)
                                               aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                     >  2:12265/1‑100 (MQ=255)
                                                 tGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATc                                                   <  1:1352376/100‑1 (MQ=255)
                                                 tGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATc                                                   <  2:1132418/100‑1 (MQ=255)
                                                  gCATTAGCTTATTTTTTTTTTATTATTTTAACCACCAAGCGGGATATGACATACCGCTCATGTTAAAGCAGCAGCGATGGGTTTATCTTGGTgcgc                                                      >  1:2424542/1‑96 (MQ=255)
                                                   cATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAc                                                 <  1:1833995/100‑1 (MQ=255)
                                                   cATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAc                                                 <  2:2567073/100‑1 (MQ=255)
                                                        gCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTccc                                            <  2:1846201/100‑1 (MQ=255)
                                                        gCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTccc                                            <  1:46869/100‑1 (MQ=255)
                                                         cTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCg                                           >  1:57391/1‑100 (MQ=255)
                                                          ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCgg                                          >  2:2844797/1‑100 (MQ=255)
                                                          ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGGTTATCGTGATGCGCATCACTTCCCgg                                          >  2:503083/1‑100 (MQ=255)
                                                           tATTTTTTTGTTCTCATGCTAACCACCCGGCGAGGTGTGAGATACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAcgataaacc                                         <  1:2190672/100‑9 (MQ=255)
                                                           tATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGc                                         <  2:396601/100‑1 (MQ=255)
                                                              ttttttGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAg                                       >  1:1931881/1‑99 (MQ=255)
                                                                   tGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGc                                 <  2:2789105/100‑1 (MQ=255)
                                                                   ggTTAGCATGATAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGc                                 >  2:695263/2‑100 (MQ=255)
                                                                    gTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGcc                                >  1:1092376/1‑100 (MQ=255)
                                                                       atgaTGCTAACCACCCGGCGAGGTGTGACACACCGCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCGGCCg                               <  1:1288403/95‑1 (MQ=255)
                                                                       aTCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTa                             <  2:764351/100‑1 (MQ=255)
                                                                       aTCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTa                             >  2:1023979/1‑100 (MQ=255)
                                                                        tCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCggga                            <  2:2220061/100‑4 (MQ=255)
                                                                        tCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCggga                            <  2:2306260/100‑4 (MQ=255)
                                                                         cATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGa                           <  2:501939/100‑1 (MQ=255)
                                                                          aTGCTACCCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCgggaag                          <  1:2317552/100‑6 (MQ=255)
                                                                          aTGCTAACCACCCGGCGATGCGTCACACACCTCGCACTTGAAATCAGCAGCGAGTGGTTTACCGTGATGCGCATTACTTCCCGGCAGTCCTGCGGTAGaa                          >  1:2603688/1‑100 (MQ=255)
                                                                          aTGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCCCTTCCCGGCAGTCCTGCCGTAGaa                          >  1:162186/1‑100 (MQ=255)
                                                                          aTGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCgggaag                          >  2:1587553/1‑95 (MQ=255)
                                                                           tGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAg                         <  2:1254729/100‑1 (MQ=255)
                                                                             cTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAGATCAGCAGCGATTGGTTTAGCGTGATGCGCATCACTTCCCGGCAGGACTGCCgggaagtga                       <  1:1587549/100‑9 (MQ=255)
                                                                             cTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGca                       <  2:1611000/100‑1 (MQ=255)
                                                                               aacccaccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCACCACTCCCCGGCAGTCCTGACGTAGAAGcaac                     <  1:705683/94‑1 (MQ=255)
                                                                                accaccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGcaaca                    >  1:659163/1‑100 (MQ=255)
                                                                                 ccaccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGcaacaa                   >  2:1495704/1‑100 (MQ=255)
                                                                                  caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACaaa                  >  2:1884861/1‑100 (MQ=255)
                                                                                  caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACaaa                  >  2:1811939/1‑100 (MQ=255)
                                                                                    cccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGTTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAAt                 >  1:1263412/1‑99 (MQ=255)
                                                                                    cccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAAtt                >  2:506767/1‑100 (MQ=255)
                                                                                     ccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAAttt               >  2:1838510/1‑100 (MQ=255)
                                                                                       tgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt             >  1:129347/2‑100 (MQ=255)
                                                                                       tgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGAAGTCCTGCCGTAGAAGCAACAAATTTc              >  1:370589/2‑99 (MQ=255)
                                                                                        gCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt             >  2:1703585/1‑99 (MQ=255)
                                                                                        gCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTg            >  1:339562/1‑100 (MQ=255)
                                                                                         cGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga           >  2:2020317/1‑100 (MQ=255)
                                                                                          gagaggtgACACTCCTCGCACTGTAGTTCAGCAGCGGTTGTTTTATCGTGATGCGGATCATTTCCTGGCAGTCCTGCCGTAGAAGCATCAAATTTCTg            <  1:809730/93‑1 (MQ=255)
                                                                                          gAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTgag          >  1:1329238/1‑100 (MQ=255)
                                                                                            ggTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc        >  2:2675694/1‑100 (MQ=255)
                                                                                            ggTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc        >  1:2381217/1‑100 (MQ=255)
                                                                                             gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACt       >  2:2818595/1‑100 (MQ=255)
                                                                                              tgtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACtt      >  2:1918415/1‑100 (MQ=255)
                                                                                               gtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg     <  1:1640536/100‑1 (MQ=255)
                                                                                               gtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg     >  1:2749026/1‑100 (MQ=255)
                                                                                               gtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCAGCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg     <  1:70653/100‑1 (MQ=255)
                                                                                                 gACACACCTCGCACTTGAAATCAGCAGCGATTGATTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCAGTAGAAGCAACAAATTTCTGAGACTTGTa   >  1:1742369/1‑100 (MQ=255)
                                                                                                  acacacCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTaa  <  2:2173071/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAA  >  NC_000913/1788238‑1788435

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