Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,344 G→A A107V (GCC→GTC)  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,3440GA94.9% 152.5 / NA 60A107V (GCC→GTC) phnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (30/26);  minor base C (3/0);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGC  >  NC_000913/4320097‑4320590
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gCGTGAGCATGGGTAAACCGGCGGCGTGCGGGTGATGTCAGCCGGCTGTGCGCCGCTATCCGCTTCAATCTTCGCCAGCCCCTGACGCGCCAGCAGTCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCTGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCAGGAATGTCTTTATAAAcgacgg                                                                                                                                                                                                                                                     <  1:129805/251‑1 (MQ=255)
      gCAGGGGTAAACCGGCGGCGTGCGGGTTATTTCATCCGGCTGTGCTCCGCTATCCTCTACAAACTTCGCCAGCCCCAGACGGGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGTTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATAc                                                                                                                                                                                                                                               <  1:6245/251‑1 (MQ=255)
          gggTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCAAGCCCCTGAAGCGCCAGCAGGCTGAACACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCCGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAACTCGAGCCGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGTCGGCCGCCGGGAATGTACTTATAAAcgacgggtgtaaggcg                                                                                                                                                                                                                                           >  2:263611/1‑241 (MQ=255)
          gggTAAACCGGAGGCGTGCGGGTGATGGCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCg                                                                                                                                                                                                                                           <  1:64748/251‑1 (MQ=255)
              aaaCCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGTCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCg                                                                                                                                                                                                                                       <  1:258080/251‑1 (MQ=255)
              aaaCCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGCGATACGGCGTTCg                                                                                                                                                                                                                                       >  1:88198/1‑251 (MQ=255)
                      ggcgTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCCGGGTAAAATCGAGCAGGCGCTGGGGGTAGTCGTCGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGGCTTTATAAACGACGGGGATACGGCGTTCGAGAAGCAt                                                                                                                                                                                                                               >  2:226767/1‑251 (MQ=255)
                         gTGCGGGAGATGTCATCCGGCTGTGCGACCCTATCCTCTTCAACCTTCGCCCGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGGGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGCAGGGTGTAAGCGTAAGTTGGGCCAAGAAGCTGGCCGCAGGGAAGGTCTTTATAAAAGACGGAGATAAGGCGTTCGAGACGCATccc                                                                                                                                                                                                                            >  2:113355/1‑251 (MQ=255)
                            cGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCggg                                                                                                                                                                                                                         >  1:196486/1‑250 (MQ=255)
                             gggTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGGGATACGGCGTTCGAGACGCATCCCggtg                                                                                                                                                                                                                        >  2:144629/1‑251 (MQ=255)
                                   tGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGCTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCg                                                                                                                                                                                                                  >  2:42391/1‑251 (MQ=255)
                                       aTCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGGCCGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTGGGGTGGGCCCAACAACTGGCCGCCGGGAATGGCTTTATAAACGAAGGGGATAAGGCGTTCGGGACGCCTCCCGGGGgtgtggaacg                                                                                                                                                                                                              >  2:58065/1‑245 (MQ=255)
                                           ggCTGTGAGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCAAGCAGGCTGAACACGTGAGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGAAAAATCGAGCCGGCACTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAAGGTCTTTATAAACGAAGGAGATACGGCGTTAGAGACGCATCCCGGCGGGGTCGAGCGGCTc                                                                                                                                                                                                          >  1:42306/1‑251 (MQ=255)
                                            gCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAAAGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCg                                                                                                                                                                                                         >  1:37095/1‑251 (MQ=255)
                                                       cTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGcccg                                                                                                                                                                                              >  1:190674/1‑249 (MQ=255)
                                                             tctTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGGAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGc                                                                                                                                                                                        >  1:87862/1‑251 (MQ=255)
                                                                    actccGCCATCCCCTGACGCGCCAGCAGGCTGAAAACGGGCGGCGTCGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCTAGCAGGGGATGGGTGTAGTCGTAGGTTGGGCCAAGGAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGTCGCAGCCCTGTGGTGGCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACgtg                                                                                                                                                                                 <  1:238396/247‑1 (MQ=255)
                                                                     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)
                                                                                                                                                                                                                                                    aaCGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGc  <  1:13333/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGC  >  NC_000913/4320097‑4320590

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