Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,416 G→A T83I (ACC→ATC)  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,4160GA98.2% 172.3 / NA 57T83I (ACC→ATC) phnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (27/28);  total (29/28)
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

CTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGC  >  NC_000913/4320166‑4320654
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
cTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGCGCGGCTCGCTTACCCCCAGCTTCGCCAACgtgg                                                                                                                                                                                                                                                >  1:261323/1‑251 (MQ=255)
                  cgcCAGCAGGCTGAAAACGTGCGGCGACGGCTGCGGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCgcagca                                                                                                                                                                                                                              <  2:13206/251‑1 (MQ=255)
                      ggaagGCTGAAAACGTGCGGCGCCGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGaa                                                                                                                                                                                                                          <  1:226767/248‑1 (MQ=255)
                                aaaCGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCCGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCCAGAAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGGGTTCGAGACGCATCCCGGTGGTGGCGAGAGGCTCGATTACCGCCCGCTTCTCCGACGTGGTGCGGGAGGCGCGCCGGAGGAAAATCGcgta                                                                                                                                                                                                                 >  2:213851/1‑247 (MQ=255)
                                   cgggcggcgACGGCTGCTGTTCGCTGTCGGCGGTGGGCAGCGCCGGCGCGTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACg                                                                                                                                                                                                             <  2:217095/248‑1 (MQ=255)
                                      gcggcgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGtt                                                                                                                                                                                                           >  1:146993/1‑250 (MQ=255)
                                        ggggacggcgtctgTTCGCTGTCGGCGGTGGTGAGCGTCGGCGCTTCGCCGTTTGCCGGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGAAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTGATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCTGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGTAAAATCGCTTCAACGTTATc                                                                                                                                                                                                        <  2:48465/240‑1 (MQ=255)
                                         gcggaggCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcg                                                                                                                                                                                                       <  2:47663/246‑1 (MQ=255)
                                         gcgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGGGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcg                                                                                                                                                                                                       >  1:73242/1‑251 (MQ=255)
                                           gACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcgcc                                                                                                                                                                                                     <  1:170265/251‑1 (MQ=255)
                                           gACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCCCCAACGGGATGCGGGAGGGGCGCCGCAGGGAAATCGCTTCAACGTTTTcgcc                                                                                                                                                                                                     >  2:125227/1‑251 (MQ=255)
                                              ggCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTTCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCt                                                                                                                                                                                                  <  1:169815/251‑1 (MQ=255)
                                              ggCTGCGGTCCGTTGTCGGCGGTGGGCAGCGTCGGCGCTGCGCCGGTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGTCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCt                                                                                                                                                                                                  <  2:87383/251‑1 (MQ=255)
                                               gctgcggtgcgctggaggcggTGATCGGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAACGAGCAGGCGAGGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAAGCGTTTCAACGTTATCGCCGCTg                                                                                                                                                                                                 <  2:45950/235‑1 (MQ=255)
                                                  gctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTCTGCCAGCACGGGAAAATCGAGCCAGCGCTGGGGGCC‑GTCGTAAGTTGCGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGCAACGCCGTTCGCGAAGCCTCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTCATGCGGTAGGCCCGCCGCCGGACAATCGCTTCAACGTTCTCCCCGCTgggc                                                                                                                                                                                              >  1:242445/1‑249 (MQ=255)
                                                   cgtttCGCTGTCGGCGGTGGTCAGCGTCGGGGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCt                                                                                                                                                                                             <  1:226847/248‑1 (MQ=255)
                                                     gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATTGGTGTA‑GTCGTATGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTg                                                                                                                                                                                            <  1:164780/250‑1 (MQ=255)
                                                            gcaggcggGGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGTTGTA‑GGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcagcg                                                                                                                                                                                    <  2:230354/248‑1 (MQ=255)
                                                            gTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcag                                                                                                                                                                                      <  1:244018/249‑1 (MQ=255)
                                                                 cggTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCCGGCGATGGGTGTA‑GACGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGAGGAGGTAGGCGCGCAGCAGGAAAATCGCTTCAAAGTTATCGCCGCTGGACTGTTTcagcgccagc                                                                                                                                                                               >  2:49211/1‑251 (MQ=255)
                                                                  ggtggtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGAGGGGTGTA‑GTCGTCGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCCGGACAATCGCTTCACCGTTATCCCCGCTGGCCTGGTTcagcgccagcg                                                                                                                                                                              >  2:34160/1‑251 (MQ=255)
                                                                     ggtCAGCGTCGGCGCTTCGCCGTTGGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGGGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAAGCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGcgccg                                                                                                                                                                           <  2:25660/251‑1 (MQ=255)
                                                                            gTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTc                                                                                                                                                                    >  1:8408/1‑251 (MQ=255)
                                                                                      cGCCGTTTGCAAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCa                                                                                                                                                          <  1:49547/251‑1 (MQ=255)
                                                                                       gCCGTTCGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGCA‑GTCGTAGGTTGGGACAAGCAACTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTCTCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCa                                                                                                                                                          >  2:236041/1‑250 (MQ=255)
                                                                                         cGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCTCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATg                                                                                                                                                       >  2:81836/1‑251 (MQ=255)
                                                                                                               gAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGAt                                                                                                                                 >  1:76217/1‑251 (MQ=255)
                                                                                                                ggCAGGCGATGGG‑GTAAGTGGTAGGTTGGGCAAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCGCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGAAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATc                                                                                                                                <  2:88198/250‑1 (MQ=255)
                                                                                                                  cAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCt                                                                                                                               >  1:89378/1‑250 (MQ=255)
                                                                                                                   gggCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGGCATCACGCGATCTAc                                                                                                                             >  2:67027/2‑251 (MQ=255)
                                                                                                                           ggTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGt                                                                                                                     <  1:196469/251‑1 (MQ=255)
                                                                                                                                  ggcgagggTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTaaact                                                                                                               <  2:76908/245‑4 (MQ=255)
                                                                                                                                     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)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
CTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTA‑GTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGC  >  NC_000913/4320166‑4320654

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