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% 1304.0 / 41.9 426intergenic (+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 (225/151);  new base T (14/13);  total (249/177)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 435 total aligned reads displayed.

GCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTT  >  NC_000913/375098‑375395
                                                                                                                                                     |                                                                                                                                                    
gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                      >  2:1029493/1‑150 (MQ=35)
 cgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTAGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaa                                                                                                                                                     >  1:5645774/1‑150 (MQ=25)
  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCGTTATCCGGCCTACGGGTGGCGCGTGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaaa                                                                                                                                                    >  1:1160570/1‑150 (MQ=255)
   gtAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAAc                                                                                                                                                   <  1:5989533/149‑1 (MQ=37)
      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAGGCGGCGTAAACGcc                                                                                                                                                >  2:1478007/1‑150 (MQ=37)
             aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                         >  2:877498/1‑150 (MQ=37)
                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                      >  1:3349908/1‑150 (MQ=37)
                  ccagcggcgcCCGACTGCCGGTTGCGGCGTGCCCGCCTCCGCCGGCCTCCCGGGGGCGCGAGACTTTGTCGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                    <  2:5628699/144‑1 (MQ=11)
                    aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                                  >  1:2750372/1‑150 (MQ=37)
                    aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAACCGGCCTACGGGTGGAGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCg                                                                                                                                  >  1:499131/1‑150 (MQ=255)
                       tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                               >  1:5671391/1‑150 (MQ=38)
                             cgacccccGGCTGCGGCGAGCACGCCTTATCCGGCCTACGGCTGGCGCGAGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGcgcg                                                                                                                         <  2:5393145/144‑1 (MQ=11)
                             cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGcgcg                                                                                                                         <  2:2447595/150‑1 (MQ=39)
                               aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCggg                                                                                                                       <  1:4860656/150‑1 (MQ=39)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGaa                                                                                                                     >  2:4991830/1‑150 (MQ=39)
                                  gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGAGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAt                                                                                                                    <  2:412/150‑1 (MQ=37)
                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAtt                                                                                                                   >  2:3508610/1‑150 (MQ=39)
                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGACttt                                                                                                                  >  1:4156418/1‑150 (MQ=255)
                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                 <  2:750581/150‑1 (MQ=39)
                                      gaggcggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTCGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                                <  2:2246147/147‑1 (MQ=255)
                                        tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGAGAGAATTTGTAGGACTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                                              <  1:4336421/150‑1 (MQ=21)
                                          cggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                            <  1:6152691/150‑1 (MQ=39)
                                            gcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGACTTTGTAGGCCt                                                                                                          >  2:1764440/1‑150 (MQ=255)
                                               tGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                       >  1:4123806/1‑150 (MQ=39)
                                               tGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                       <  2:2360541/150‑1 (MQ=21)
                                                 aaCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATaa                                                                                                     >  1:6712936/1‑150 (MQ=39)
                                                 aaCGCCTTATCCGGCCTACGGGTGGCGAGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATAAGGCGTTTACGCCGCATCCGGcag                                                                                                                                     >  1:3121342/1‑116 (MQ=11)
                                                    gCCTTATGAGGGCTACGGTTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAc                                                                                                  <  2:4403925/150‑1 (MQ=255)
                                                      cTTATCCGGCCTACGGTTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTTATAAGAcgc                                                                                                <  2:5173965/150‑1 (MQ=255)
                                                      cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                <  2:4784431/150‑1 (MQ=39)
                                                       ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                               >  1:5131074/1‑150 (MQ=39)
                                                            cGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGAGGGACTTTGTAGGCCTGATAAGACGCGTTAGc                                                                                          >  2:4023061/1‑150 (MQ=255)
                                                            cGGCCTACGGGTGGCGAGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                                          <  2:3457630/150‑1 (MQ=37)
                                                                cgaCGGGTGGCGCGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCg                                                                                      <  1:1283716/148‑1 (MQ=12)
                                                                    gggTGGCGCGAGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                                                  <  1:3897421/150‑1 (MQ=255)
                                                                     ggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 <  2:1530061/150‑1 (MQ=38)
                                                                     ggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 >  1:3075272/1‑150 (MQ=38)
                                                                     ggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 <  1:1451888/150‑1 (MQ=38)
                                                                     ggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 >  1:3739319/1‑150 (MQ=38)
                                                                     ggTGGCGCGAGAATTTGAAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                 >  1:4764860/1‑150 (MQ=38)
                                                                      gTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                                >  1:1978752/1‑150 (MQ=38)
                                                                          cgcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                            <  2:2750372/150‑1 (MQ=37)
                                                                           gcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                           >  1:4631337/1‑150 (MQ=35)
                                                                            cgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCt                                                                          >  1:974485/1‑150 (MQ=35)
                                                                            cgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCt                                                                          >  1:5316132/1‑150 (MQ=35)
                                                                                caTTTGTAGGCCTGATAAGACGCGTTTGCGTCGGATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcac                                                                      <  1:1546857/149‑1 (MQ=255)
                                                                                  tttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                                                    <  1:5192765/150‑1 (MQ=37)
                                                                                    tGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGTTTGGCGCGGGATTTTGTAGGCCTGATAAGACGCGTTTGCGTCGCATCAGGCATCTGCACACGa                                                                  <  1:2378724/150‑1 (MQ=255)
                                                                                       aGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAATCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTGAGCGTCGCATCAGGCAGCTGCACACGACTg                                                               >  2:6527208/1‑150 (MQ=255)
                                                                                         gCCTGATACGACGCGTTAGCGTCGCATCTGGCTTCTGCGCACGACTGCCGGATGCGGCTTAAACGCCTTATCGGGCCTACGGATGGCGGGGGAATTTTTAGGCCTTATAAGACGGGTTAGGGCCGCATCTGGCATCTGCACACGACTGcc                                                             <  1:2554394/150‑1 (MQ=11)
                                                                                         gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGcc                                                             >  1:3182447/1‑150 (MQ=35)
                                                                                          ccTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                                            >  1:3524488/1‑150 (MQ=35)
                                                                                            tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAttt                                                                                                                  >  1:909486/1‑94 (MQ=33)
                                                                                            tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGa                                                          <  2:2116002/150‑1 (MQ=34)
                                                                                            tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGa                                                          >  1:2871183/1‑150 (MQ=34)
                                                                                             gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                         >  1:4000340/1‑150 (MQ=34)
                                                                                                 aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                                                     >  2:4704422/1‑150 (MQ=35)
                                                                                                  gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGAt                                                    >  1:1729530/1‑150 (MQ=35)
                                                                                                  gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGAt                                                    >  1:2039443/1‑150 (MQ=35)
                                                                                                   aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACt                                                                <  2:2647201/137‑1 (MQ=34)
                                                                                                       gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGc                                                                                                                                                 <  2:5955122/52‑1 (MQ=35)
                                                                                                        ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              >  2:3068382/1‑150 (MQ=40)
                                                                                                        ctAGCGGCGCATCCGGCATCCCCGCCCGACTGCCGGACGCGGCGTACCCGCCTATTCCGGCCTCCGGCTCGCGCGGGAATTTGTAGGCCTGAGAAGACGCGTTAGATTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              <  2:1927408/149‑1 (MQ=14)
                                                                                                           gCGACGCATCTGTCAACTGCGCTCCCCTGCCGGATGCGGCTCACACGCCTTATCCGGCCTACGGATGGCGAGCGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCtt                                           <  2:4173252/150‑1 (MQ=14)
                                                                                                           gCGACGCATCAGGCACCTGCGCACGACTGCCGTAGGGGGGGTAAACGGCTTATGGGGCCTAGGGGTGGGGCGGGCCTTTGTAGGCCTTATAAGACGCGTTTGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCtt                                           <  1:6493413/150‑1 (MQ=255)
                                                                                                                  aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAACTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGc                                    <  2:3337628/150‑1 (MQ=255)
                                                                                                                     aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                                 >  2:4645173/1‑150 (MQ=255)
                                                                                                                     aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTa                                 >  1:1928489/1‑150 (MQ=255)
                                                                                                                      ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTAc                                >  1:5863780/1‑150 (MQ=255)
                                                                                                                      cgCATCCCCGCACCACTGCCCCAACCGGAGATAACGCGTGATCGGTGCTCCGTATGGCGCGCTACCTTGAAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTAc                                <  1:4297048/149‑1 (MQ=9)
                                                                                                                       gCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATACACGTCTTGTCCGGCCTACa                               >  2:6250786/1‑150 (MQ=255)
                                                                                                                       gCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACa                               >  1:542811/1‑150 (MQ=255)
                                                                                                                          tCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttc                            <  1:1024522/150‑2 (MQ=255)
                                                                                                                           cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTc                           <  2:4313055/150‑1 (MQ=255)
                                                                                                                            tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCg                          >  2:2088458/1‑150 (MQ=255)
                                                                                                                            tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCg                          >  1:2706036/1‑150 (MQ=255)
                                                                                                                             gcgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc                         <  1:217227/150‑1 (MQ=255)
                                                                                                                              agcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGGGGGAATTTGTAGGACTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCTCACGACTGCCGGATGCGATAAGCGACAAGTCCGGCCTACAGTTCGcc                        >  2:1646368/2‑150 (MQ=14)
                                                                                                                               gcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                        >  2:3169456/1‑85 (MQ=33)
                                                                                                                               gcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGccc                       >  2:5489230/1‑150 (MQ=255)
                                                                                                                                cACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCg                      >  1:6628709/1‑150 (MQ=255)
                                                                                                                                 aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGt                     >  2:6095258/1‑150 (MQ=255)
                                                                                                                                 aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGt                     <  1:4704400/150‑1 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCa                >  2:3559515/1‑150 (MQ=255)
                                                                                                                                       gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAGTTCGCCCGTAGGCAg               >  2:3730900/1‑150 (MQ=18)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTATTAATTATAATAATCTTTATCTTAGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt              >  2:5908839/1‑150 (MQ=18)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGt              >  2:1363752/1‑150 (MQ=255)
                                                                                                                                         cGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTc             >  1:6223044/1‑150 (MQ=255)
                                                                                                                                         cGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTc             <  1:2497092/150‑1 (MQ=255)
                                                                                                                                          gtaTGCGTAGTAAACCCCTTATCCGCCCTTCGGATGTCGCGTGAATTTGTACTCCTCATAAGACGAGTTAGCGTCGAATCAGGAATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCATACATTTCGCCCGTAGGCAGTCa            <  2:1160570/148‑1 (MQ=255)
                                                                                                                                           gATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAt           >  2:2950801/1‑150 (MQ=255)
                                                                                                                                            ataaggcgTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAtt          >  1:1093784/5‑150 (MQ=18)
                                                                                                                                              aaggcgTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaa        >  1:340223/3‑150 (MQ=255)
                                                                                                                                               cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa       <  1:4065385/150‑1 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGAAAAACTTCATGTACGCCCTACAATTCTCCCGTAGGAAGTAATTAAAt      >  2:6697124/1‑150 (MQ=14)
                                                                                                                                                 gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCg                          >  1:806175/1‑129 (MQ=255)
                                                                                                                                                 gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATa     >  2:5211061/1‑150 (MQ=255)
                                                                                                                                                 gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATa     >  1:4915254/1‑150 (MQ=255)
                                                                                                                                                  cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                              >  1:2104475/1‑108 (MQ=255)
                                                                                                                                                   gTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt   >  1:45743/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTT  >  NC_000913/375098‑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: