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

CTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGC  >  NC_000913/4320166‑4320654
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
cTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGCGCGGCTCGCTTACCCCCAGCTTCGCCAACgtgg                                                                                                                                                                                                                                                >  1:261323/1‑251 (MQ=255)
                  cgcCAGCAGGCTGAAAACGTGCGGCGACGGCTGCGGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCgcagca                                                                                                                                                                                                                              <  2:13206/251‑1 (MQ=255)
                      ggaagGCTGAAAACGTGCGGCGCCGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGaa                                                                                                                                                                                                                          <  1:226767/248‑1 (MQ=255)
                                aaaCGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCCGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCCAGAAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGGGTTCGAGACGCATCCCGGTGGTGGCGAGAGGCTCGATTACCGCCCGCTTCTCCGACGTGGTGCGGGAGGCGCGCCGGAGGAAAATCGcgta                                                                                                                                                                                                                 >  2:213851/1‑247 (MQ=255)
                                   cgggcggcgACGGCTGCTGTTCGCTGTCGGCGGTGGGCAGCGCCGGCGCGTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACg                                                                                                                                                                                                             <  2:217095/248‑1 (MQ=255)
                                      gcggcgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGtt                                                                                                                                                                                                           >  1:146993/1‑250 (MQ=255)
                                        ggggacggcgtctgTTCGCTGTCGGCGGTGGTGAGCGTCGGCGCTTCGCCGTTTGCCGGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGAAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTGATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCTGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGTAAAATCGCTTCAACGTTATc                                                                                                                                                                                                        <  2:48465/240‑1 (MQ=255)
                                         gcggaggCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcg                                                                                                                                                                                                       <  2:47663/246‑1 (MQ=255)
                                         gcgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGGGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcg                                                                                                                                                                                                       >  1:73242/1‑251 (MQ=255)
                                           gACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATcgcc                                                                                                                                                                                                     <  1:170265/251‑1 (MQ=255)
                                           gACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCCCCAACGGGATGCGGGAGGGGCGCCGCAGGGAAATCGCTTCAACGTTTTcgcc                                                                                                                                                                                                     >  2:125227/1‑251 (MQ=255)
                                              ggCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTTCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCt                                                                                                                                                                                                  <  1:169815/251‑1 (MQ=255)
                                              ggCTGCGGTCCGTTGTCGGCGGTGGGCAGCGTCGGCGCTGCGCCGGTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGTCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCt                                                                                                                                                                                                  <  2:87383/251‑1 (MQ=255)
                                               gctgcggtgcgctggaggcggTGATCGGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAACGAGCAGGCGAGGGGTGTAGTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAAGCGTTTCAACGTTATCGCCGCTg                                                                                                                                                                                                 <  2:45950/235‑1 (MQ=255)
                                                  gctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTCTGCCAGCACGGGAAAATCGAGCCAGCGCTGGGGGCCGTCGTAAGTTGCGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGCAACGCCGTTCGCGAAGCCTCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTCATGCGGTAGGCCCGCCGCCGGACAATCGCTTCAACGTTCTCCCCGCTgggc                                                                                                                                                                                              >  1:242445/1‑249 (MQ=255)
                                                   cgtttCGCTGTCGGCGGTGGTCAGCGTCGGGGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCt                                                                                                                                                                                             <  1:226847/248‑1 (MQ=255)
                                                     gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATTGGTGTAGTCGTATGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTg                                                                                                                                                                                            <  1:164780/250‑1 (MQ=255)
                                                            gcaggcggGGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGTTGTAGGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcagcg                                                                                                                                                                                    <  2:230354/248‑1 (MQ=255)
                                                            gTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcag                                                                                                                                                                                      <  1:244018/249‑1 (MQ=255)
                                                                 cggTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCCGGCGATGGGTGTAGACGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGAGGAGGTAGGCGCGCAGCAGGAAAATCGCTTCAAAGTTATCGCCGCTGGACTGTTTcagcgccagc                                                                                                                                                                               >  2:49211/1‑251 (MQ=255)
                                                                  ggtggtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGAGGGGTGTAGTCGTCGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCCGGACAATCGCTTCACCGTTATCCCCGCTGGCCTGGTTcagcgccagcg                                                                                                                                                                              >  2:34160/1‑251 (MQ=255)
                                                                     ggtCAGCGTCGGCGCTTCGCCGTTGGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGGGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAAGCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGcgccg                                                                                                                                                                           <  2:25660/251‑1 (MQ=255)
                                                                            gTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTc                                                                                                                                                                    >  1:8408/1‑251 (MQ=255)
                                                                                      cGCCGTTTGCAAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCa                                                                                                                                                          <  1:49547/251‑1 (MQ=255)
                                                                                       gCCGTTCGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGCAGTCGTAGGTTGGGACAAGCAACTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTCTCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCa                                                                                                                                                          >  2:236041/1‑250 (MQ=255)
                                                                                         cGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCTCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATg                                                                                                                                                       >  2:81836/1‑251 (MQ=255)
                                                                                                               gAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGAt                                                                                                                                 >  1:76217/1‑251 (MQ=255)
                                                                                                                ggCAGGCGATGGGGTAAGTGGTAGGTTGGGCAAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCGCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGAAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATc                                                                                                                                <  2:88198/250‑1 (MQ=255)
                                                                                                                  cAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCt                                                                                                                               >  1:89378/1‑250 (MQ=255)
                                                                                                                   gggCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGGCATCACGCGATCTAc                                                                                                                             >  2:67027/2‑251 (MQ=255)
                                                                                                                           ggTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGt                                                                                                                     <  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)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
CTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGC  >  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: