Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,543 G→A L41F (CTT→TTT)  phnI ← carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,5430GA98.4% 198.2 / NA 65L41F (CTT→TTT) phnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (32/31);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCG  >  NC_000913/4320298‑4320787
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
gTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGATTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAactgct                                                                                                                                                                                                                                                 >  1:97506/1‑251 (MQ=255)
   ggTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGtt                                                                                                                                                                                                                                              >  1:263649/1‑251 (MQ=255)
         gCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCt                                                                                                                                                                                                                                        <  1:3609/251‑1 (MQ=255)
          ccAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTg                                                                                                                                                                                                                                       >  1:176713/1‑251 (MQ=255)
              gcagcTGGCCGCCGGGAATGTCTTTAAAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGGGGGGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCg                                                                                                                                                                                                                                   <  2:72142/251‑1 (MQ=255)
                 gcggcccgccgGGAATGTCTTGATAACCGGCGGAATTCCGGCGTTCGAGACCCACCCCGGTGGTGTCGGGCGGCTCGCTTCCCGCCAGCTTCGCCAACGGGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCTCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACg                                                                                                                                                                                                                                <  2:134656/246‑1 (MQ=255)
                            ggAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCggg                                                                                                                                                                                                                     >  2:13333/1‑251 (MQ=255)
                             gAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATACCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCCGGACAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGGTAAACTGCTGTTCAATCTGGGCGACGATCAGTTCGGGc                                                                                                                                                                                                                    >  2:246080/1‑251 (MQ=255)
                              aaTGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGc                                                                                                                                                                                                                    >  1:101384/1‑250 (MQ=255)
                               atggcTTTATAAACGACGGAGAGACGGCGTTCGAGACGCACCCCGGTGGTCGCGAGCGGGCCGCTTACCGCCAGCTTCGCCAAAGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACGGCTGTGCAATCTGGGCGACGCTCAGTTCGGGCaa                                                                                                                                                                                                                  <  2:90868/247‑1 (MQ=255)
                                     ttATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCCGTTCGGGCAAATCGGt                                                                                                                                                                                                            >  1:236633/1‑251 (MQ=255)
                                        taAACGACGGAGACCCGGCGTTCGAGACGAATCACGGTCGTGTCGAGAGGCTCGCTTACCGCCAGATTCGCCAACGTGATGAGGTAGGCGCGCAGCAGGAAACTCGCTTCACCGTTATCGCCGCTGGCCTGTTTCAGCGACAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCCATCTACCGAGCGGTTAAACTGCTGTTCAATCTGGGCGCCGCTCAGTTCGGGCAACTCGGTATc                                                                                                                                                                                                         >  2:175239/1‑251 (MQ=255)
                                         aagcgccggAGATGCGGCGTTCGCGACGCGTCCCGGTGGTAGCGAGCGGCTCGCATCCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTGTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCAGGGCGACGCTCAGTTCGGGCAAATCGGTATCg                                                                                                                                                                                                        <  2:6146/245‑1 (MQ=255)
                                          aaCGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGc                                                                                                                                                                                                       <  1:13333/251‑1 (MQ=255)
                                            cgtgggAGATACGGCGTCCGAGACGCCTCCCGGGGGTGTCGAGCGGCTCGCTTACCGCCAGCTCCGCCAACGTTATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACTTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCt                                                                                                                                                                                                     <  2:261323/247‑1 (MQ=255)
                                              acggAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCACCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCg                                                                                                                                                                                                   >  2:138083/2‑251 (MQ=255)
                                                   gatgcGGCGTTCGAGACCCCCCCGGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTc                                                                                                                                                                                              <  1:144629/247‑1 (MQ=255)
                                                   gaTACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTc                                                                                                                                                                                              >  1:50928/1‑251 (MQ=255)
                                                         gCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGctct                                                                                                                                                                                        >  1:229729/1‑251 (MQ=255)
                                                                 gaCGAATCCCGGTAGTGTCGAGCCGCTCGCTTACCGCCAGCTTAGCCAACGGGCTGCGGGAGGCGCGCAGCAGGACAATCGCTTCCCAGTTCTAGCCGCTGGACTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCCATGCCGCCTTCGGTCATCACGAGATCTACCGAGAGGTTAAACTGCTGTTCACTCTGGGCGACGCTCCGTTCGGGCAAATCGGTTTCGCCTCGGCGTCGGCTCTCTTGCAgg                                                                                                                                                                                >  2:42145/1‑251 (MQ=255)
                                                                    gCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCg                                                                                                                                                                             >  1:13990/1‑251 (MQ=255)
                                                                           ggtggcGTCGAGCGCCTCGCTTACCACCAGCTTCGCACCCGCGATGCGGTAGGAGCGCAGCAGGACAATCGCTTCACCGTTATCGCCGCTGGCCTGTTTAAGCGCCAGAGCCGCCAGTTCGCGGGCGGCCATGCCGCCTTCGGTCATCACGCACTCTACCGCGAGGATAAACTGCTGGTCAATATGGGCGACGCCCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGATCTCTTGCAGGGCGTGggcgg                                                                                                                                                                      >  2:180803/1‑251 (MQ=255)
                                                                              ggtgTCGGGGGGCTCGCTTACCGCCAGCTTCGCAAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTAATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGt                                                                                                                                                                   <  1:229492/251‑1 (MQ=255)
                                                                                           ccGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTtctc                                                                                                                                                      <  2:76217/250‑1 (MQ=255)
                                                                                                      cAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTc                                                                                                                                            >  1:68463/1‑250 (MQ=255)
                                                                                                          ttCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGc                                                                                                                                       >  1:63744/1‑251 (MQ=255)
                                                                                                                aacgtgaggCGGTAGGCGCGCAGCAGGAAAATCGCTCCACCGTTATCGCCGCTGGCCCGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTa                                                                                                                                 <  2:146993/243‑1 (MQ=255)
                                                                                                                     gatgCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATc                                                                                                                             <  1:34160/248‑1 (MQ=255)
                                                                                                                       tgcggtcggaacgcAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCCGCGCCGCCAGTTCGCGGTCGGCAATAACGCCTTCGGTCATCACGCGATCTAACGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGCCTCTCTTGCAGGGCGTGGGCGGCGTCGCTCGCCTTATCGCCCCCTTTCACGGCCACGTACATccgc                                                                                                                          >  2:124205/12‑248 (MQ=255)
                                                                                                                         cGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGcac                                                                                                                        >  2:172487/1‑251 (MQ=255)
                                                                                                                           gTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCAGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCCCGGCAACGGACATCAGCacac                                                                                                                      >  1:150559/1‑251 (MQ=255)
                                                                                                                                           aaaaTCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGAGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTAACGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCCACGTACATCAGCACACCTCCACATGAGTGGtt                                                                                                      >  2:216002/1‑251 (MQ=255)
                                                                                                                                            aaaTCGCTTCAACGTTAGCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCGACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTc                                                                                                     <  1:138083/251‑1 (MQ=255)
                                                                                                                                            aaaTCGCTCCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCGTCGGTCATCACGCGATCTAACGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTGCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGGCGATCGCCTTCTCGCCCCCTTTCACGGCAACGGACAGCAGCACACCTCCACATGAGTGGTTc                                                                                                     <  2:54085/251‑1 (MQ=255)
                                                                                                                                              aTCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCCACGTACATCAGCACACCTCCACATGAGTGGTTcgc                                                                                                   >  2:89616/1‑251 (MQ=255)
                                                                                                                                                      aaCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGCTCTACCGCGGGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGAATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCAAATCAGTGGCTCGCGGAATAGc                                                                                           >  2:186565/1‑251 (MQ=255)
                                                                                                                                                      aaCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGc                                                                                           <  1:213851/251‑1 (MQ=255)
                                                                                                                                                         gTTATCCCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCCGTTCGCGCCCGGCAATGCCGCCTTCGGTCATCCCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGcc                                                                                          >  2:41003/1‑249 (MQ=255)
                                                                                                                                                            aTCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAACCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGCTCGCCTTCTCGCCCCCTTTCCCGGCACCGTACATCAGCACACCTCCCCATGAGTGGTTCGCGGAATAGCcagc                                                                                       >  2:218126/1‑249 (MQ=255)
                                                                                                                                                            aTCGCCGCGGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAAGCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGGCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCcagcag                                                                                     <  2:87862/251‑1 (MQ=255)
                                                                                                                                                               cccgcGGGCCTTTTTCATCGCCGGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTGCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGgcg                                                                                  <  2:13990/250‑1 (MQ=255)
                                                                                                                                                                 cgcTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGACTAGCCAGCAGGCGCt                                                                                >  2:93528/1‑251 (MQ=255)
                                                                                                                                                                 cgcTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCt                                                                                >  2:47742/1‑251 (MQ=255)
                                                                                                                                                                              cagcgccagcgcCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTGCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCa                                                                    <  1:263611/250‑1 (MQ=255)
                                                                                                                                                                                ccgcccgcgccGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGa                                                                 <  2:196486/250‑1 (MQ=255)
                                                                                                                                                                                        gccgccAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCCGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTACGCGGACTAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCg                                                         >  2:21690/1‑251 (MQ=255)
                                                                                                                                                                                        gccgccAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCg                                                         >  1:166991/1‑251 (MQ=255)
                                                                                                                                                                                        cccgccAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCg                                                         <  2:263649/250‑1 (MQ=255)
                                                                                                                                                                                         ccgccAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACGCCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCg                                                         >  1:26196/1‑250 (MQ=255)
                                                                                                                                                                                            ccAGTTCGCGGTCGGCAATGCCGCCTGGGGGCGTCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAAGCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGGTCGCCACAGGTCAGGATCAGGTCGATGc                                                     <  2:190674/251‑1 (MQ=255)
                                                                                                                                                                                                  cgcgGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCcgcgcg                                               >  1:14225/1‑251 (MQ=255)
                                                                                                                                                                                                        cGGCAATGCCGCCTTCGGGCATCACGCGATCGACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGTGGGAAc                                          <  1:216002/250‑1 (MQ=255)
                                                                                                                                                                                                           cAATGCCGCCTTCAGTCATCACGCGAGTGTCCGTGAGGTTAAGCTTCTGGTCAACCTCGGCGACGCTCAGTTGGGGCAAATCGGTATCGCTGCGGCGGCGGCTCTCTTGCAGGGCTGGGGCGGCGACGAGCGCCTTCTCGCACGCTTTCACGGCAACGTACAACAGCACACCTCCACATGAGTGTTTCTCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGATTGGGTACGGAt                                      <  2:64147/251‑1 (MQ=255)
                                                                                                                                                                                                                  gccTTCGGTCATCACGCGGCCTACAGCGAGGTAAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCACGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGCCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGAGGAATAGCCAGCAGGCGCACGCCACAGGTCAGGATCAGGTCGATGCCGCGCGGGAACGGATGCGGgcg                               >  2:236132/1‑251 (MQ=255)
                                                                                                                                                                                                                    cTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCCATCTGTGCGACGCTCAGTTCGGGCAAATCGGTCTCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTACACATGAGAGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGGTGACGAGAGGGAACGGATGCGGGCGCt                             >  2:63204/1‑251 (MQ=255)
                                                                                                                                                                                                                      tCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTc                            <  1:113355/250‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTCATCACGCGATCTACCTCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGGACATCAGCACACCGCCACAGGAGTGGGTCGCGGAATAGCCAGCAGGCGCGCGCCACAGGTCAGTATCAGGTCGATGCCGAGCGGGAGTGGAGGCGGGGGCTCg                           <  2:263464/250‑1 (MQ=255)
                                                                                                                                                                                                                        ggTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGt                         <  1:172487/251‑1 (MQ=255)
                                                                                                                                                                                                                         gTCATCACGCGAGCTACCGCGAGGTTAAACTGCTGGGCAAGCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCAGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCTTCGATCTCCTTCTCGCCCCCGTTCACGGCAACGTACATCAGCACACCGCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGTAACGGATGCGGGCGCTCGGTg                        <  1:159756/251‑1 (MQ=255)
                                                                                                                                                                                                                              caCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTc                   <  1:81836/251‑1 (MQ=255)
                                                                                                                                                                                                                              caCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTGCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCGCTTGCAGGGCGTGGGCGGCGTCGATCGCCGTCTCTCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCATGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGtt                    <  1:238525/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aCCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAg             <  1:58065/247‑1 (MQ=255)
                                                                                                                                                                                                                                         aCCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATg         <  1:125227/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCAc      >  1:71735/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggAGGTTAAACTGCTGGTCAATCTGGGCGACGCTCAGTTCGGGCAAAGCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACt     <  2:101384/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ggTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCg  >  2:184227/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCG  >  NC_000913/4320298‑4320787

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