Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 512,793 G→A S73N (AGC→AAC)  ybaT → putative transporter YbaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913512,7930GA96.3% 86.7 / ‑2.9 27coding (218/1293 nt)ybaTputative transporter YbaT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/11);  minor base . (1/0);  total (16/11)
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.55e-01

TTTAAGGGCTGATCATGATGAACACGGAAGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAGCAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGGAGCACCTTATTTTGCTCTACGCGTTGGGGATCATTGCGGTGATGACGCTTTTCAACTCCTTAAGCAACCATGCGGT  >  NC_000913/512562‑513033
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
tttAAGGGCTGATCATGATGAACACGGAAGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCCTTATCGAc                                                                                                                                                                                                                                >  1:185278/1‑251 (MQ=255)
              atgatgAACACGGAAGGGAATAAAGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCAATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGAATTGCTAATGGGAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTAGCGATGTTTTCCGGTTTTGACTATGCGCGTCTAGGGGCGAGCTCTACCCACAATGGCGGCATTATATACTTCTTTCTTCGCgg                                                                                                                                                                                                                  >  2:222784/1‑251 (MQ=255)
                     acacGGAAGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGGCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGCCTTCTTTCGTCGCGGATTAGGc                                                                                                                                                                                                           >  1:114594/1‑251 (MQ=255)
                           aaGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGGTGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTCTCCCAACAATTGCGGCGATTTCGACTTTTTTCGTCTCGGGTTTGGCAACGGc                                                                                                                                                                                                     >  2:113510/1‑251 (MQ=255)
                            aGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCg                                                                                                                                                                                                    >  1:10301/1‑251 (MQ=255)
                                    cctgTAACAAACCCCTCGGTCTATGGAACGTCGTTTCCTTCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCg                                                                                                                                                                                            <  2:70480/248‑1 (MQ=255)
                                     cGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAAGCAACGGCGTCTTTTCGc                                                                                                                                                                                           >  1:171753/1‑251 (MQ=255)
                                            aaaCCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCTTCGCGGATTAGGCAACGGCGTCTTTTTGCTTgcgcc                                                                                                                                                                                    >  1:205632/1‑250 (MQ=255)
                                               ccTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCtctc                                                                                                                                                                                 >  1:79675/1‑251 (MQ=255)
                                                     ggtagaTGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACt                                                                                                                                                                           <  2:43682/246‑1 (MQ=255)
                                                     ggTCTATGGAACGTCTTTTCCATCGGCTGTGGGGCAATGGGGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACt                                                                                                                                                                           <  2:78803/251‑1 (MQ=255)
                                                                               cagtGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCa                                                                                                                                                 <  1:29486/248‑1 (MQ=255)
                                                                                                gggCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCGGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGt                                                                                                                                <  1:219921/251‑1 (MQ=255)
                                                                                                          cTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTTCTGTACCTGTTGACGCTGGCGGGGAACAATGCCCGGGTCGCCCGTGGTTTTTgcg                                                                                                                     >  2:1208/1‑251 (MQ=255)
                                                                                                            tCGCGCTGCTGGGGCAGGCTGCTTTGCTAATGGACGCCTCGACCTGGGTCTCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGTGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTg                                                                                                                        <  2:184194/246‑1 (MQ=255)
                                                                                                                    ctgGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCCTGGTCGCCCGGGGCTTTTGCGCTTATGCCgg                                                                                                           >  2:69953/1‑250 (MQ=255)
                                                                                                                    ctgGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGGTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGGCGCCCGTGCTTTTGGCGCTTATGCCgg                                                                                                           >  1:280057/1‑250 (MQ=255)
                                                                                                                      ggggCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGGGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTTTGCCgggc                                                                                                         >  2:339622/1‑248 (MQ=255)
                                                                                                                       gggCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTGCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCa                                                                                                        <  2:33589/251‑1 (MQ=255)
                                                                                                                           aGGCTGCATTGCTAATAGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGAAGTTTTACGGTTATGCCTATGCGCGTATGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGGGATCTCGTTACTGTCCCTGTTGACGCCGGCGGTGAGCATCGCCCGGGGGTCCCGTGGTTTTGGGGCTTATGCCGTggcgttt                                                                                                    >  2:342914/1‑245 (MQ=255)
                                                                                                                                gCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCa                                                                                               <  1:170171/251‑1 (MQ=255)
                                                                                                                                                                              gTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGGAGCACCTTATTTTGCTCTACGCGTTgggg                                                 <  1:285122/251‑1 (MQ=255)
                                                                                                                                                                                        gATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGGAGCACCTTATTTTGCTCTACGCGTTGGGGATCATTGCgg                                       >  1:327277/1‑251 (MQ=255)
                                                                                                                                                                                              ttCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTAGGCCGTGCAGTTTTTGCATGAAGTCAGCCAGGAGGAGTACCTTATTTTGCTCTACGCGTTGGGTAACGTTGCGTtgatg                                  <  2:54663/250‑1 (MQ=255)
                                                                                                                                                                                                      aTGCCTATGCGCGTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGGTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGTAGCTCCTTATTTTGCTCTACGCGTTGTGGATCATTGCGGTGATGACGCTTTTc                         <  2:151072/251‑1 (MQ=255)
                                                                                                                                                                                                                 cgTCTGGGGGCGAGCTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCCGTTTTTGCCTGAAGGCAGCCAAGAGGAGCACCTTATTTTGCTCTACGCGTTGGGGGTCATTGCGGTGGTGACGCTTTTCAACTCCTTa                >  2:95811/1‑249 (MQ=255)
                                                                                                                                                                                                                               cTATCCCAACAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGGAGCACCTTATTTTGCTCTACGCGTTGGGGATCATTGCGGTGATGACGCTTTTCAACTCCTTAAGCAACCATGCGGt  >  1:37298/1‑249 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
TTTAAGGGCTGATCATGATGAACACGGAAGGTAATAACGGTAACAAACCTCTCGGTCTATGGAACGTCGTTTCCATCGGCATTGGGGCAATGGTGGGGGCGGGGATCTTCGCGCTGCTGGGGCAGGCTGCATTGCTAATGGAAGCCTCGACCTGGGTCGCCTTTGCTTTTGGCGGTATTGTGGCGATGTTTTCCGGTTATGCCTATGCGCGTCTGGGGGCGAGCTATCCCAGCAATGGCGGCATTATCGACTTCTTTCGTCGCGGATTAGGCAACGGCGTCTTTTCGCTGGCGCTCTCGTTACTGTACCTGTTGACGCTGGCGGTGAGCATCGCCATGGTCGCCCGTGCTTTTGGCGCTTATGCCGTGCAGTTTTTGCATGAAGGCAGCCAGGAGGAGCACCTTATTTTGCTCTACGCGTTGGGGATCATTGCGGTGATGACGCTTTTCAACTCCTTAAGCAACCATGCGGT  >  NC_000913/512562‑513033

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