Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,247 A→T 5.6% 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,2470AT 5.6% 723.0 / 11.8 253intergenic (+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 (129/88);  new base T (4/10);  total (137/116)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

GCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAG  >  NC_000913/375098‑375393
                                                                                                                                                     |                                                                                                                                                  
gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                    >  1:3099851/1‑150 (MQ=35)
 cgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaa                                                                                                                                                   <  2:375282/150‑1 (MQ=14)
  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGACTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaca                                                                                                                                                  >  2:2050216/1‑148 (MQ=25)
  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaaa                                                                                                                                                  >  1:1329022/1‑150 (MQ=37)
  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaaa                                                                                                                                                  >  1:3471378/1‑150 (MQ=37)
       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCt                                                                                                                                             >  2:2635152/1‑150 (MQ=37)
           gcctCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATCAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                         <  2:1744739/147‑1 (MQ=255)
            cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                        <  1:55272/150‑1 (MQ=37)
            cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                        <  2:3082109/150‑1 (MQ=37)
             aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGAAAACGCCTTATCCg                                                                                                                                       >  1:4115490/1‑150 (MQ=37)
             aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGGTGGCGCGAGAATTGGTAGGCCTGATCAGACGCGTTAGCGTCGCATCAGGCAGCTGCGCACGACTGCCGGAGGCGGCGGCCACGCCGTATCCg                                                                                                                                       >  2:854018/1‑150 (MQ=11)
              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                      <  1:1271606/150‑1 (MQ=37)
                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                    >  2:1621673/1‑150 (MQ=37)
                    aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                                <  1:2852625/150‑1 (MQ=37)
                    aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                                <  2:3540639/150‑1 (MQ=37)
                       tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                             >  2:716320/1‑150 (MQ=38)
                       tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                             <  1:2949038/150‑1 (MQ=38)
                          gcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgc                                                                                                                          >  2:3901530/1‑150 (MQ=39)
                           cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                         >  1:4247266/1‑150 (MQ=39)
                           cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                         >  1:409550/1‑150 (MQ=39)
                           cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                         >  2:3811479/1‑150 (MQ=39)
                            aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcgc                                                                                                                        <  2:210541/150‑1 (MQ=39)
                              gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCgg                                                                                                                      >  1:1966083/1‑150 (MQ=39)
                              gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGAATTATAAGGACTAAGGATggagatt                                                                                                                      >  2:2309633/1‑145 (MQ=14)
                               aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCggg                                                                                                                     <  2:1685558/150‑1 (MQ=21)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGac                                                                                                                   >  1:3698354/1‑149 (MQ=39)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGaa                                                                                                                   >  1:2298342/1‑150 (MQ=39)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGaa                                                                                                                   >  1:163881/1‑150 (MQ=39)
                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGCGCGGGAAt                                                                                                                  >  2:4266674/1‑150 (MQ=14)
                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGACtt                                                                                                                 >  2:2158611/1‑150 (MQ=255)
                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                               >  1:2041323/1‑150 (MQ=39)
                                      gATGCGGCGTGAACGCCTTATCCGTCCTGCGGCTTGCGCGTGAATTTGTAGGCCAGTTAAGACTCGTTAGCGTCGCAACTGGCTTCTGCTCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                              <  2:1059085/150‑1 (MQ=11)
                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGTTGGCGTGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                              <  2:2423052/150‑1 (MQ=17)
                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                              >  2:1794111/1‑150 (MQ=39)
                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                              <  1:1371364/150‑1 (MQ=39)
                                       aTGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                             <  1:3586507/150‑1 (MQ=39)
                                            gcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAGTTGTAGGCCt                                                                                                        >  2:3147399/1‑150 (MQ=255)
                                              gTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGa                                                                                                      >  2:2616270/1‑150 (MQ=39)
                                               tGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                     >  2:3720970/1‑150 (MQ=39)
                                                gAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                    <  2:3460646/150‑1 (MQ=39)
                                                  aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTATTAGGCCTGATAAGAAGCGTTATCGTCGCATCAGGCATCTGCGCTCGAATTCCGGATGCGGCATAAACGCCTTATCCGGCATACGGATGGCGCGGGAATTTGTATGCCTGATAAg                                                                                                  <  1:3027881/150‑1 (MQ=14)
                                                    gCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAc                                                                                                >  1:1519327/1‑150 (MQ=39)
                                                      cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                              >  1:1231820/1‑150 (MQ=39)
                                                       ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                             >  2:1213673/1‑150 (MQ=255)
                                                       ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                             >  2:1213742/1‑150 (MQ=39)
                                                       ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                             >  1:1113671/1‑150 (MQ=39)
                                                       ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                             >  1:3642876/1‑150 (MQ=39)
                                                          tCCGGCCTACGGTTGGCGAGAGAATTTGTAGGCCTGATTAGACGCGTTATCGTAGAATCAGGCATCTGCGGACGAATTCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgctaac                                                                                          <  2:2250646/150‑5 (MQ=11)
                                                          tCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgctaac                                                                                          <  2:4013790/150‑5 (MQ=255)
                                                          tCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                                                          <  1:26413/150‑1 (MQ=39)
                                                                cTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCg                                                                                    <  2:591711/150‑1 (MQ=21)
                                                                 tACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                                                   >  1:3253565/1‑150 (MQ=255)
                                                                  aCGGGTGGCGCGAGAATTTGTAGGCCTGAAAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGAAAACGCCATATCCGGACTAAGGAAGGCGCGGGAATTTGTAGGCCTGAAAAGACCCGTTAGCGTCGCa                                                                                  >  1:1311975/1‑150 (MQ=14)
                                                                    gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                                                >  1:3136770/1‑150 (MQ=39)
                                                                      gTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                              >  1:3413999/1‑150 (MQ=38)
                                                                           gcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGTCATc                                                                         >  1:2960790/1‑150 (MQ=25)
                                                                           gcgAGAATTTGTAGGCCCGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                         >  1:3887373/1‑150 (MQ=25)
                                                                               gaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCAGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGca                                                                     >  2:3056312/1‑150 (MQ=25)
                                                                                 aattctAGGCCTGAGAAGACGCGGTTGCTTCTGAACTGGCATCCGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcaca                                                                   <  1:2429526/145‑1 (MQ=11)
                                                                                 aTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcaca                                                                   <  1:77837/150‑1 (MQ=37)
                                                                                   ttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACg                                                                 >  1:2751053/1‑150 (MQ=37)
                                                                                          ccTGATAAGACGCGTTATCGTCGCATCTGGCCTCTGCGCCCGACTGCCGTATGCGGCGTCACCTCATTATCCGGCCTACGGTTGGGGCGGGAATTTGTAGGCCTTATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                                          <  1:610222/150‑1 (MQ=11)
                                                                                           cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCgg                                                         >  1:4147104/1‑150 (MQ=34)
                                                                                           cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCgg                                                         <  2:645390/150‑1 (MQ=34)
                                                                                             gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                       >  1:3411776/1‑150 (MQ=34)
                                                                                                aaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCg                                                    >  2:2526878/1‑150 (MQ=34)
                                                                                                 aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                   <  2:4103658/150‑1 (MQ=35)
                                                                                                 aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                   <  1:3901530/150‑1 (MQ=35)
                                                                                                 aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCATTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                   <  2:3588721/150‑1 (MQ=25)
                                                                                                  gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGAt                                                  >  1:2469524/1‑150 (MQ=35)
                                                                                                     gcgTTTGCGTCGCCTGAGGCATCTGCGCTCGACTGCCGGTAGCGGCGTCACCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaaa                                               <  2:3317788/150‑1 (MQ=11)
                                                                                                     gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                        >  1:3081884/1‑61 (MQ=35)
                                                                                                     gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                              <  2:501866/103‑1 (MQ=32)
                                                                                                      cgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                              <  1:4018079/150‑1 (MQ=39)
                                                                                                       gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                          <  1:2729482/121‑1 (MQ=32)
                                                                                                       gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACAACTGCCGGATGCGATAAACg                                             >  1:3547993/1‑150 (MQ=255)
                                                                                                                 gaTCTGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                   <  2:3386523/149‑1 (MQ=255)
                                                                                                                  aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGc                                  >  1:1100015/1‑150 (MQ=255)
                                                                                                                  aTCAGGCACCTGCGCACGACTGCCGGATTCGGCGTAAACGCCTTATGCGGCCTACGGCTGGCGCGTGAATTTGTAGGACTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGc                                  <  2:3660578/150‑1 (MQ=17)
                                                                                                                     tggCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                               <  1:649510/149‑1 (MQ=255)
                                                                                                                     aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                               >  2:4036009/1‑150 (MQ=255)
                                                                                                                       gCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACa                             >  2:1378280/1‑150 (MQ=255)
                                                                                                                            cgcgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCg                        <  1:902302/149‑1 (MQ=255)
                                                                                                                             gcgcTCGACTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc                       <  1:3794136/150‑1 (MQ=255)
                                                                                                                             gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                   >  2:1926943/1‑138 (MQ=255)
                                                                                                                             gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                   <  1:1926943/138‑1 (MQ=255)
                                                                                                                             gcgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc                       <  1:1869600/150‑1 (MQ=255)
                                                                                                                             cagcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc                       <  2:1583502/148‑1 (MQ=11)
                                                                                                                               gcACGACTCGCGGCTGCTGCGTTCTCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATGGCGccc                     <  2:3012465/150‑1 (MQ=17)
                                                                                                                                 aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGt                   >  2:3675861/1‑150 (MQ=255)
                                                                                                                                   gACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAg                 >  1:2547117/1‑150 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGTCGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTTCATTTCGCCCGTAGGct              >  1:2612219/1‑149 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCa              >  1:479498/1‑150 (MQ=255)
                                                                                                                                        cgggATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGATTTTGTTGGCCTGATAAGACGCGTTAGCGTCGCATCTGGCTTCTGCACACGACTGCCGGTTGCGATTAACGTCTTGTCGGGCCTACATTTCGCCCGTAGGCAGt            <  2:1311975/148‑1 (MQ=18)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt            >  1:2140597/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt            >  2:1515381/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt            >  1:1704396/1‑150 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTc           >  2:3989685/1‑150 (MQ=18)
                                                                                                                                          ggATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGTTAAGACGCGTTGGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGGCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa          >  1:321887/1‑150 (MQ=255)
                                                                                                                                                 gcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc                       >  2:287808/1‑130 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAG  >  NC_000913/375098‑375393

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: