Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 331,529 T→C 7.7% intergenic (+33/‑842) betT → / → yahA choline transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913331,5290TC 7.7% 148.8 / 3.4 52intergenic (+33/‑842)betT/yahAcholine transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAA  >  NC_000913/331234‑331829
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
tGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATAgc                                                                                                                                                                                                                                                                                                         >  2:457601/1‑301 (MQ=255)
ggggCGAAGGGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGgc                                                                                                                                                                                                                                                                                                         <  1:320458/300‑1 (MQ=255)
    ggAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGGAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCAGAAAACCGCGCCATGTCTGCATATGGCGCTg                                                                                                                                                                                                                                                                                                     <  1:66292/300‑1 (MQ=255)
        gagCAAAATTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCGACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTtc                                                                                                                                                                                                                                                                                                 <  2:535966/301‑1 (MQ=255)
          gCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTACCTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTtctc                                                                                                                                                                                                                                                                                               >  2:153882/1‑301 (MQ=255)
           aaaaaCTGTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTcc                                                                                                                                                                                                                                                                                              <  2:371307/300‑1 (MQ=255)
             aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTa                                                                                                                                                                                                                                                                                            >  2:103809/1‑301 (MQ=255)
               attttGTGTTTTAGTTTTGGCCGCAGCCATATTCGGTGCCGGGCTTTACTTACCGCGCACGCAGCGGTAAATCGACCTACGACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCGTGACCAGTACGAGCGGCACCTTAACTTTATTCAGCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATAGGGCGCTGTTTCTCCTAcc                                                                                                                                                                                                                                                                                          <  1:150957/299‑1 (MQ=255)
                cTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCAAGTCTGCATATGGCCCTGTTTCTCCTAccc                                                                                                                                                                                                                                                                                         >  1:211638/1‑300 (MQ=255)
                           ggttttGGCGGCAGCATTTTTCGGTGCGGTGCTTTCCCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAATCCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGtt                                                                                                                                                                                                                                                                              <  1:624970/298‑1 (MQ=255)
                                 ggCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCt                                                                                                                                                                                                                                                                        <  2:211881/301‑1 (MQ=255)
                                          aaTATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGAAGTTCCCGGACGCGTGATTGTTATTGCATACAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATAAGTTATATCTATACGGttt                                                                                                                                                                                                                                                               >  1:569348/1‑300 (MQ=255)
                                                    gCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGt                                                                                                                                                                                                                                                     <  2:509484/301‑1 (MQ=255)
                                                     ccGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTACGCCCTTactt                                                                                                                                                                                                                                                    >  1:543729/1‑298 (MQ=255)
                                                           tttACCTACCGCGGACGGGGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCGGTTAGAAGGCAGACAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTGAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTAtt                                                                                                                                                                                                                                              <  1:177480/301‑1 (MQ=255)
                                                                 aacggcggaggcatgtgTAAATCGACCTGTTACGTGCTGTAAACCTTCCTGGTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCGTCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCGTAGTATCGATTgatgat                                                                                                                                                                                                                                        <  1:112458/285‑1 (MQ=255)
                                                                                ggTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTCTATCTATACGGTTAAGCCATTAGTATCTATTGATGATTACCAGacaacaaaa                                                                                                                                                                                                                         >  1:158038/1‑298 (MQ=255)
                                                                                         tcttccggtctgCTGGAACCCTTCCTGTTAGAAGGCATCCAGGGCAACGACCTTATGTACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATaaaaaa                                                                                                                                                                                                                <  1:518734/290‑1 (MQ=255)
                                                                                                    gCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATaaaaaaagatcggaaga                                                                                                                                                                                                     >  2:519406/1‑290 (MQ=255)
                                                                                                       ggAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATACAACAGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACACGACAACAAATAATAAAAAAAAAAAAAACAAttt                                                                                                                                                                                                  >  1:421898/1‑300 (MQ=255)
                                                                                                       ggAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGAGTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTg                                                                                                                                                                                                  <  1:103685/301‑1 (MQ=255)
                                                                                                          aaCCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGGGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAg                                                                                                                                                                                               <  1:457037/301‑1 (MQ=255)
                                                                                                                   gTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTCGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAAca                                                                                                                                                                                      >  2:328956/1‑300 (MQ=255)
                                                                                                                           gCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTg                                                                                                                                                                              <  2:455937/301‑1 (MQ=255)
                                                                                                                                    gCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGc                                                                                                                                                                     <  2:123975/301‑1 (MQ=255)
                                                                                                                                           ccTGAGGGACGTCAGCTAAGTGCAGGTGATCTCCGATATTCTTGACCGGTACGAGCGGCACCTTAGTTTTATTCTTCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTGCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGGTATAGCTATACGGTTAAGCACTTAGTATCTATTGATGATTGCCAGAGAACAGATAATAAAAAAGAAAAGAACTATGGCAGCCCAAAAGCGACAATTGGGCTGTTGCGAATGtt                                                                                                                                                              <  1:123786/301‑1 (MQ=255)
                                                                                                                                                  cccgtcgctcttccgagtttgagtaCCGATATTCTTGACCAGTACGAGCGGCACCTTAACGTTATTCATCGCCATCGTGAGGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTTATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAg                                                                                                                                                       <  1:278400/277‑1 (MQ=255)
                                                                                                                                                       cAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAg                                                                                                                                                  >  2:25305/1‑301 (MQ=255)
                                                                                                                                                        aGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGt                                                                                                                                                 <  2:570062/301‑1 (MQ=255)
                                                                                                                                                         gCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTAGGGCTGTTGCGAATTTTAAATAAGTTTAGta                                                                                                                                                >  2:608674/1‑300 (MQ=255)
                                                                                                                                                           aaaGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCtt                                                                                                                                              >  2:470260/1‑301 (MQ=255)
                                                                                                                                                                   ggTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCAAAAACCTACATTTGGGCTGATGCGAATGTTCAATAAGTTTAGTATtaattaatt                                                                                                                                      >  1:92094/1‑300 (MQ=255)
                                                                                                                                                                     tGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAgagatcggaagagcacacgt                                                                                                                                    >  2:278703/1‑282 (MQ=255)
                                                                                                                                                                     tGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTa                                                                                                                                    >  2:619087/1‑301 (MQ=255)
                                                                                                                                                                       aTCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTaaa                                                                                                                                  >  2:629772/1‑301 (MQ=255)
                                                                                                                                                                             gATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGc                                                                                                                            <  2:544402/301‑1 (MQ=255)
                                                                                                                                                                                        ccAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGa                                                                                                                 >  2:480107/1‑301 (MQ=255)
                                                                                                                                                                                          aGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGaaa                                                                                                               >  2:377860/1‑301 (MQ=255)
                                                                                                                                                                                             ttttctttttttttatctgttgtctggtaatcatcaatagatactaagggcttaaccgtatagatataacgtatcaagaggtgggagaaacaGCGCCATATGCAGACATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGa                                                                                                            >  1:58600‑M1/93‑301 (MQ=255)
                                                                                                                                                                                               ttcttttttattatctgttgtctggtaatcatcaatagatactaagggcttaaccgtatagatataacgtatcaagaggtgggagaaacaGCGCCATATGCAGACATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGATGATCATTTGAAAtgaaa                                                                                                          <  2:58668‑M1/211‑3 (MQ=255)
                                                                                                                                                                                                             cTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATgaga                                                                                            >  2:393772/1‑301 (MQ=255)
                                                                                                                                                                                                              tttATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGaa                                                                                           <  1:25254/301‑1 (MQ=255)
                                                                                                                                                                                                                cgccttctggtcgttcttccgatctggCCATTGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAt                                                                                         <  1:478786/276‑1 (MQ=255)
                                                                                                                                                                                                                 aTTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAtg                                                                                        >  2:54300/1‑301 (MQ=255)
                                                                                                                                                                                                                     tggtcCTTCGTGAAGGGCCGGTCCTGAGTGTTATGTTTCCGGATGGGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCACTGTTTCTCCTACCTCTTTATACGGTATATCTATACGGTGAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATGTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTa                                                                                    <  1:274424/298‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAtagatcggaagagcacacgtctgaaa                                                                >  2:479387/1‑276 (MQ=255)
                                                                                                                                                                                                                                                           tccctacacgacgctctttcgatctaaacaGCGCCATATGCAGACATGGCTCGGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTAtt                                              >  1:549264‑M1/31‑301 (MQ=255)
                                                                                                                                                                                                                                                                      aTTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAAt                                   >  2:101423/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       ttGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAAtt                                  >  2:603284/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    aaacaGCGCCATATGCAGACATGGCGCGGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAActgaaattaacgcaaccattagttattagatcggaagagcacacatctgaaa                     <  2:549937‑M1/296‑26 (MQ=255)
                                                                                                                                                                                                                                                                                                 tCTGCATATGGCGCTGTTTCTTCTACCTCTTGATTCGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATg        <  1:377364/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 tCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATg        <  2:422461/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       tatGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCTTGAGAAATGTGTATCGTAAATCAACTGAAATTAAAGAAACCATTTGTTATTAAAGGTTTAATTATCTGTATGAGATATTTTATTGAATGTTTTaa  >  1:624602/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
TGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAA  >  NC_000913/331234‑331829

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