Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,751 C→T A202A (GCG→GCA yaaW ← putative enzyme‑specific chaperone YaaW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,7510CT95.3% 123.4 / ‑4.1 44A202A (GCG→GCAyaaWputative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (28/13);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTAC  >  NC_000913/10509‑10989
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
catcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                                                                                                                                                                        <  1:307868/251‑1 (MQ=255)
                  gAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCAATGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACccgc                                                                                                                                                                                                                      >  1:671098/1‑251 (MQ=255)
                        atTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATCATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGACTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGATTTAACCCCATTTAATTCTGCGCCCACAGGGCCTCACAGACccgcgccgcg                                                                                                                                                                                                                >  2:594424/1‑251 (MQ=255)
                                    aacagTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCTCGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCTCGCAGCAAACCAt                                                                                                                                                                                                     <  1:544018/246‑1 (MQ=255)
                                        gTAGTGTTTTCAGGGCGCGGCCTAATAATCAGCCAGTGGGGCAGTGTCTCCGCTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGAACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCCTTTAATGCTGCGCCCACAGAGCCTCCCAGACCCGACCCGCGCAGCAAACCATGcacc                                                                                                                                                                                                >  2:573100/1‑248 (MQ=255)
                                        gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGccc                                                                                                                                                                                                 >  1:543978/1‑250 (MQ=255)
                                        gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGccc                                                                                                                                                                                                 >  1:428882/1‑250 (MQ=255)
                                                    gggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCCAGCACGCTCAtt                                                                                                                                                                                    >  1:175103/1‑251 (MQ=255)
                                                                  aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCg                                                                                                                                                                      >  1:340001/1‑251 (MQ=255)
                                                                   atCAGCCAGTGGGGCAGTGTCTTCGAGCTTTGGAGGGGAAAATGAAAATTTGCCCCGGGTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGTGCTTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACTCTCATTGCTGCGTGGGTGCGt                                                                                                                                                                     <  2:543978/251‑1 (MQ=255)
                                                                     cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCCCTGCCGCTGACCGCTTTAACCCCATTTATTGCTGCGCACACAGGGCCTCCCAGaccc                                                                                                                                                                                                                        >  2:13605/1‑198 (MQ=255)
                                                                     cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGaccc                                                                                                                                                                                                                        <  1:13605/198‑1 (MQ=255)
                                                                       gCCAGTGGGGCAGTGTCTACGATCTTGTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTaaaa                                                                                                                                                                 <  2:479303/251‑1 (MQ=255)
                                                                        ccAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAACTTTTCCCCGGTTTACGGTGTCAGCCCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGCTTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCACATTTAATGCTGCGCCCACAGGGCATCCCAGACCCCCTCCGCGCAGCAAACCATGCCCAAGCCCGCTCATTGCTACGAGGGGGCGTaatat                                                                                                                                                                >  2:594745/1‑248 (MQ=255)
                                                                          aGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTTCTGCGTGGCGGCGTAAAATGc                                                                                                                                                              >  2:92522/1‑251 (MQ=255)
                                                                                      tCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGGGGGTGCGGAAAATGCGGGTCAATTGGc                                                                                                                                                  >  2:86445/1‑251 (MQ=255)
                                                                                           gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCCCGCTCATTGATGCg                                                                                                                                                                              >  2:233777/1‑218 (MQ=255)
                                                                                           gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCg                                                                                                                                                                              <  1:233777/218‑1 (MQ=255)
                                                                                                             tgaaaatTTTCCCCGGTTTACGGTGTCCGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGAACGCTTTAACCCCATTTAATGCTGCGCCCCCAGGGCCTCCCCGACCCGAGCCGCGCAGCAAACCAAGCCCCAGCACGCTCATTGCTGCGTGGGt                                                                                                                                                                         >  2:27616/1‑205 (MQ=255)
                                                                                                             tgaaaatTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGt                                                                                                                                                                         <  1:27616/205‑1 (MQ=255)
                                                                                                                       ccccGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCt                                                                                                                                                                                                                                 >  1:378857/1‑139 (MQ=255)
                                                                                                                          cGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGACAGGCACTGCCGCTGGCCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATACCCAAGCACGCCAATTGCCGCGTGGGAGCGGCAACTGCGGGTCAATTGCCTGGAAAGCAAATGCGACACCCCtttacccatttatt                                                                                                              >  2:580999/1‑240 (MQ=255)
                                                                                                                           ggTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGCCACGCCTTTTGCCAATAAttt                                                                                                             >  1:477370/1‑251 (MQ=255)
                                                                                                                                               aGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgcc                                                                                         >  1:421635/1‑250 (MQ=255)
                                                                                                                                                      cTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCTCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCAAGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCcgcagcgc                                                                                  >  2:303943/1‑249 (MQ=255)
                                                                                                                                                       tAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCAGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGAGACCCGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGctct                                                                                 >  2:169097/1‑251 (MQ=255)
                                                                                                                                                              ccGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCt                                                                          >  1:562704/1‑251 (MQ=255)
                                                                                                                                                                ggCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTc                                                                         >  1:242494/1‑250 (MQ=255)
                                                                                                                                                                              aTTTGAAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCGGACCGCTTTAACCCCATTTAATGCTTCGCCCACCGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAATCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCAt                                                           <  2:200148/250‑1 (MQ=255)
                                                                                                                                                                                  caaGTACGGCTGGAAGTGTCACGCGATAGGCACGGCCGCTGACCGCTTTAACCCCATCTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCTCGCTCAGTGCTGCGTGGGTGCGTAAAATGCGTGTCAAGGGGCTGGAAAGCAAATGCGACACGCCGTTTGCCAATAATTTGTCTTTCATCAGCCGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGAATCGAc                                                       <  2:242494/248‑1 (MQ=255)
                                                                                                                                                                                         ggCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCa                                               <  1:93534/251‑1 (MQ=255)
                                                                                                                                                                                                attgTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTc                                        >  1:654241/4‑251 (MQ=255)
                                                                                                                                                                                                   tctccgcgATAGGCACTACCTCTAACCGCTTTAACCACATTTACATCTCGCCCTCCAGGGCCGCCCATGCCCGCGCCTCACACAAAACCATGCCCAAGCACGCTCATTGCTGCGTGGTTGCTTAAAATGCGGGTCAATTGGCTTGAAAGCAAATGCGACACGCCTTTTGCCAACTATTTGTCTTCCATCAGCAGCGGCAGCAGCTCGTCCAGCGTAGTCACCCTGGCATCGACCGCGTGCAGAAAcgcca                                      <  2:344785/246‑5 (MQ=255)
                                                                                                                                                                                                     cACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGc                                    <  1:92522/250‑1 (MQ=255)
                                                                                                                                                                                                         cgATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCCGCTCTTCCAGCTCATTCACCATGGCATCGACCGCGTGCAGAAACTCCTGCTTATg                               >  2:9616/1‑251 (MQ=255)
                                                                                                                                                                                                              ggcacTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCCCAGCGCCTCCCAGCCCCGCGCCGAGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCACTTGGCTGGAAAGCAACTGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAtactccc                          >  2:679355/5‑245 (MQ=255)
                                                                                                                                                                                                                  cTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACAATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTACTGCTACTGTTCCTCGTc                      >  2:447729/1‑251 (MQ=255)
                                                                                                                                                                                                                    gCCGCTGACCGCCTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCa                    <  1:17093/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccGCTGACCGCTTTAACGCTAATTAATGCGGCGCCCTCAGGGCCTCCCAGACCCGCGCCTCGCAGCAAACCATGCCCAAGCACGCTCATTGCTTCGTGGGTGCGTAAAATGCGGGTCAATGGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAt                   <  2:489802/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTGAACGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCACGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTTCTCGTCCAtt                  >  2:662177/1‑251 (MQ=255)
                                                                                                                                                                                                                       gCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGGAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAtat                 <  1:575848/251‑3 (MQ=255)
                                                                                                                                                                                                                                tttAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAgg        >  1:278852/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccccATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTa   <  1:594424/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cccATTTAATGCTGCGCCCACAGGGCATCCCAGACCCGCGCCGCGCAGCAAACCACGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCCTCGCCCGCGCGCAGACACTCCAGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTAc  >  2:366859/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
CATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTAC  >  NC_000913/10509‑10989

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