Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 331,524 C→G 7.7% intergenic (+28/‑847) 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,5240CG 7.7% 169.8 / 4.2 52intergenic (+28/‑847)betT/yahAcholine transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base G (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 = 9.87e-01

TTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATG  >  NC_000913/331224‑331823
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
ttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCACGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTc                                                                                                                                                                                                                                                                                                             >  1:535293/1‑301 (MQ=255)
  ggTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTcag                                                                                                                                                                                                                                                                                                           >  2:274727/1‑299 (MQ=255)
          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)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
TTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTCTCCTACCTCTTGATACGTTATATCTATACGGTTAAGCCCTTAGTATCTATTGATGATTACCAGACAACAGATAATAAAAAAGAAAAGAACTATTGCAGCCCAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATG  >  NC_000913/331224‑331823

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: