Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,247 A→G 5.4% intergenic (+366/‑212) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,2470AG 5.4% 1779.4 / 49.3 583intergenic (+366/‑212)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (326/193);  new base G (6/24);  total (341/242)
Fisher's exact test for biased strand distribution p-value = 7.72e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 592 total aligned reads displayed.

CGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTT  >  NC_000913/375099‑375395
                                                                                                                                                    |                                                                                                                                                    
       gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                              >  2:8082698/1‑150 (MQ=37)
        tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAACCGCCttc                                                                                                                                             >  1:7879395/1‑149 (MQ=37)
          ggttttcgCATCAGCGCACGACTGCCGGATGCGGCTTGAACTCCTTATCGGGCCTACGGTTGGCTCGAGAATTTGTAGGTCTGAGAAGACGAGTTTGCGTCGCTTCAGGCATCTGCGCACGACTGCAGGATGCGGCGTAAACGCCTTATc                                                                                                                                           <  2:7327794/143‑1 (MQ=11)
            atctcgCATCTGAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                         <  1:3312073/145‑1 (MQ=255)
            aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                         <  1:3043511/150‑1 (MQ=37)
              cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                       <  1:3806082/150‑1 (MQ=37)
                 gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGACTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                    >  2:5729354/1‑150 (MQ=37)
                 gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                    >  1:1397325/1‑150 (MQ=37)
                 gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAAGTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                    >  1:7223486/1‑150 (MQ=37)
                 gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGGGTGGCGCGGGAAGTTGTAGGCCTGATAAGACGCGTGAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAGGCGGCGTAGACGCAGTATCCGGCCTa                                                                                                                                    >  2:5877865/1‑150 (MQ=11)
                  cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTAc                                                                                                                                   >  1:815899/1‑150 (MQ=37)
                   aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                                  <  2:371979/150‑1 (MQ=37)
                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                               <  1:8225696/150‑1 (MQ=38)
                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                               >  1:5049084/1‑150 (MQ=38)
                        cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGATAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATgg                                                                                                                             <  1:4956555/150‑1 (MQ=255)
                               cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                      >  1:363874/1‑150 (MQ=39)
                                  ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGATGGCGCGGGAAtt                                                                                                                   >  1:3700573/1‑150 (MQ=255)
                                  ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAtt                                                                                                                   >  1:1713451/1‑150 (MQ=39)
                                  ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAtt                                                                                                                   >  1:6856177/1‑150 (MQ=39)
                                    ggATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                 <  2:7577050/150‑1 (MQ=39)
                                      aTGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                               <  1:6697499/150‑1 (MQ=39)
                                       tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                                              >  1:5780517/1‑150 (MQ=39)
                                       tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                                              <  2:4105239/150‑1 (MQ=39)
                                         cggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                            >  1:4397581/1‑150 (MQ=39)
                                           gcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                                          <  2:7900380/150‑1 (MQ=39)
                                           gcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                                          >  2:1971696/1‑150 (MQ=39)
                                             gTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAATTGTAGGCCTGa                                                                                                        >  1:4685998/1‑150 (MQ=255)
                                                aaCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATaa                                                                                                     >  1:1984075/1‑150 (MQ=39)
                                                  cGCCTTATCCGGCCTTCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCTGGCATCGGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                   <  2:1593558/150‑1 (MQ=255)
                                                  cGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                   <  2:1494399/150‑1 (MQ=39)
                                                   gCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAATCGCCGTATCCGGCCTACGGATGGCGCGGGAAGTTGTAGGCCTGATATGAc                                                                                                  >  2:3995874/1‑150 (MQ=255)
                                                     cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                >  1:4066528/1‑150 (MQ=39)
                                                           cGGCCTCCGGGTGGCGCGAGTATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGTCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                                          <  2:4399454/150‑1 (MQ=255)
                                                           cGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                                          <  1:7320399/150‑1 (MQ=21)
                                                            ggCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                         >  2:5127353/1‑150 (MQ=39)
                                                             gCCTGCGGTTGGCGCGAGAATTTGTAGGCCTGATCCGACGCGTTCGAGTCGCATCAGGCATCAGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                        <  2:337282/150‑1 (MQ=11)
                                                              ccTACGGGTGGCGCGAGAATTTGTAGGCCAGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                                                       <  1:6342475/150‑1 (MQ=255)
                                                               cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCg                                                                                      >  1:5155850/1‑150 (MQ=39)
                                                                tACGGTTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                                                     <  2:3023828/150‑1 (MQ=37)
                                                                tACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                                                     >  2:7808806/1‑150 (MQ=39)
                                                                  cGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAGTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                                                   >  1:4566359/1‑150 (MQ=255)
                                                                    ggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 >  2:5366500/1‑150 (MQ=38)
                                                                     gTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                                >  2:2755246/1‑150 (MQ=38)
                                                                     gTGGCGCGAGAATGTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                                <  2:3829074/150‑1 (MQ=38)
                                                                         cgcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                            <  2:6574294/150‑1 (MQ=37)
                                                                         cgcgAGAATTTGTAGGCCTGATAAGACGAGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                            <  2:7110595/150‑1 (MQ=37)
                                                                          gcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCagc                                                                           >  2:4645770/1‑148 (MQ=25)
                                                                           cgAGAATGTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCt                                                                          <  2:6798053/150‑1 (MQ=11)
                                                                              gaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGca                                                                       <  2:2491141/150‑1 (MQ=11)
                                                                              gaATGTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGGGCACGAGTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGca                                                                       <  1:1007210/150‑1 (MQ=255)
                                                                                 tttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                                                    >  1:2779776/1‑150 (MQ=37)
                                                                                 tttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                                                    <  2:236995/150‑1 (MQ=37)
                                                                                     tAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACt                                                                <  1:6130700/150‑1 (MQ=17)
                                                                                          cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCgg                                                           >  1:2199314/1‑150 (MQ=34)
                                                                                            gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                         >  1:803379/1‑150 (MQ=34)
                                                                                            gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                         >  1:3950155/1‑150 (MQ=34)
                                                                                            gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAGTTGTAGGCCTGATCAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                         >  2:4611621/1‑150 (MQ=21)
                                                                                             aTAAGACGCGTTAGCGTCGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCCGCGGCACACGACTGCCGGaag                                                        >  2:2831730/1‑148 (MQ=11)
                                                                                                aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                     >  2:6884388/1‑150 (MQ=35)
                                                                                                aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                     >  1:3754436/1‑150 (MQ=35)
                                                                                                aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                     >  2:4923216/1‑150 (MQ=35)
                                                                                                  aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                   >  1:5242098/1‑150 (MQ=37)
                                                                                                    ggggtAGCGTCGTATCCGGCTTCTGCGCTCGAGTGCCGGCGGGGGGGTTCACGCGTTATCCGGCCTACGGCTGGCGGGGGAATTTGTCGGCCTGATATGACGCGTTAGCGTCGCATCTGGCATCTGCACACGACTGCCGGATGCGATaaa                                                 <  1:5990402/146‑1 (MQ=255)
                                                                                                    gcgagaGCGTCGCATCACGCATCAGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaaa                                                 <  2:5539252/145‑1 (MQ=255)
                                                                                                    gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                                          <  2:6171527/93‑1 (MQ=33)
                                                                                                    gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaaa                                                 >  1:5579746/1‑150 (MQ=38)
                                                                                                      gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAGGATGGCGCTGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAACAGGCATCTGCACACGACTGCCGGATGCGATAAACg                                               >  1:2513448/1‑150 (MQ=255)
                                                                                                       ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              <  2:1713513/150‑1 (MQ=40)
                                                                                                       ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              <  2:6965720/150‑1 (MQ=40)
                                                                                                       ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              >  2:3588979/1‑150 (MQ=40)
                                                                                                        tAGCGTCGCATCAGGCATCTGCGCACGTCTGCCGGTTGGGGCGTAAAGGCGTTATCCGGACTACGTATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTc                                             <  1:3451455/150‑1 (MQ=14)
                                                                                                           cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTg                                          <  2:8325154/150‑1 (MQ=255)
                                                                                                            gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                         >  1:2883025/1‑150 (MQ=255)
                                                                                                            gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                         >  1:3193663/1‑150 (MQ=255)
                                                                                                              cacATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTcc                                       >  1:5438376/3‑150 (MQ=255)
                                                                                                               gCATCAGGCATCGGCGCTCGTCTGCCGGATGCGGCTTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGGGTTAGCGTCGCATCAGGCATCTGCTCACGACTGCCGGATGCGATAAACGTCTTGTCCg                                      <  1:8340084/150‑1 (MQ=14)
                                                                                                                cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAATACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                     <  1:1519011/150‑1 (MQ=255)
                                                                                                                cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                     <  2:4807724/150‑1 (MQ=255)
                                                                                                                 aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGc                                    <  1:911499/150‑1 (MQ=255)
                                                                                                                  tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                   >  1:310450/1‑150 (MQ=255)
                                                                                                                    aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                                 >  2:3109492/1‑150 (MQ=255)
                                                                                                                    aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                                 >  1:2070165/1‑150 (MQ=255)
                                                                                                                    aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                                 <  2:7196384/150‑1 (MQ=255)
                                                                                                                         tccgcgcACGACTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCTGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                            <  2:11060/147‑1 (MQ=255)
                                                                                                                           tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCTGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCg                          <  1:4815098/150‑1 (MQ=255)
                                                                                                                               cACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCg                      >  1:1107118/1‑150 (MQ=255)
                                                                                                                                 cgccTGCCGGATGCGGCGTAAACGCCGTAGCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTa                    <  1:4427850/147‑1 (MQ=255)
                                                                                                                                 cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTa                    <  2:6452977/150‑1 (MQ=255)
                                                                                                                                  gACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGATAAACGACTAGTCCGGCCTACAGTTCGCCCGGAg                   >  2:3215787/1‑150 (MQ=17)
                                                                                                                                   accgCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                  <  1:4288049/147‑1 (MQ=255)
                                                                                                                                   aCTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGCGCGGGATTTTGTAGGCCTGATAAGACGCGATCGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                  <  1:128069/150‑1 (MQ=255)
                                                                                                                                      gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg               <  1:2828480/150‑1 (MQ=255)
                                                                                                                                       gcGGTTGGGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt              <  1:6535198/149‑1 (MQ=255)
                                                                                                                                       ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt              >  1:18663/1‑150 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt              >  1:7655018/1‑150 (MQ=18)
                                                                                                                                         ggATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa            <  2:2666380/150‑1 (MQ=255)
                                                                                                                                           aTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAtt          >  2:2208690/1‑150 (MQ=255)
                                                                                                                                               ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt      >  2:7528776/1‑150 (MQ=255)
                                                                                                                                                gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATa     >  2:1279940/1‑150 (MQ=255)
                                                                                                                                                gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATa     >  1:2352043/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTT  >  NC_000913/375099‑375395

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: