Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,704 A→G 13.3% S107P (TCG→CCG)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,7040AG13.3% 73.8 / 4.9 30S107P (TCG→CCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (9/17);  new base G (2/2);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

TCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGA  >  NC_000913/1608566‑1608844
                                                                                                                                          |                                                                                                                                            
tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                                                                     <  3:128561/148‑1 (MQ=255)
     gCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgta                                                                                                                                <  3:175585/148‑6 (MQ=255)
              cGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                                                                      <  1:294821/149‑1 (MQ=255)
                                 tGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                                                                    <  3:312514/148‑1 (MQ=255)
                                 tGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                                                                    <  4:120069/148‑1 (MQ=255)
                                      tAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCa                                                                                                                                      <  3:104323/109‑1 (MQ=255)
                                      tAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCa                                                                                                                                      >  4:104323/1‑109 (MQ=255)
                                          gAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                                                          <  1:61544/149‑1 (MQ=255)
                                              cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                                                      <  2:264217/149‑1 (MQ=255)
                                                   cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccc                                                                                  >  3:165629/1‑148 (MQ=255)
                                                         ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                                                            >  4:43132/1‑148 (MQ=255)
                                                             tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                                                        <  3:391515/148‑1 (MQ=255)
                                                                 tGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgcg                                                                    <  3:184828/148‑1 (MQ=255)
                                                                             gTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                                                                   <  3:54115/105‑1 (MQ=255)
                                                                             gTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                                                                   >  4:54115/1‑105 (MQ=255)
                                                                                   cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt                                                  >  4:371176/1‑148 (MQ=255)
                                                                                     tAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAc                                                >  3:11743/1‑148 (MQ=255)
                                                                                          aGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCa                                          <  1:76021/149‑1 (MQ=255)
                                                                                           gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCa                                          <  4:11743/148‑1 (MQ=255)
                                                                                              ggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                                                                    >  4:27248/1‑87 (MQ=255)
                                                                                              ggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGATTTTACCTCGCGCAGTAACGACCGTGa                                                                                                    <  3:27248/87‑1 (MQ=255)
                                                                                                    aGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGaa                                <  1:268662/149‑1 (MQ=255)
                                                                                                                 cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGaaaaa                    <  4:93250/148‑1 (MQ=255)
                                                                                                                   ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatg                  <  4:165629/148‑1 (MQ=255)
                                                                                                                    gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTCTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatgg                 <  4:253005/148‑1 (MQ=255)
                                                                                                                              ttACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGaca      >  2:182691/1‑149 (MQ=255)
                                                                                                                                 acattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACAcc    >  4:9158/4‑148 (MQ=255)
                                                                                                                                   cTTGCCGAATCCGGCAATGTCGTTTTTGCCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCga  >  4:248262/1‑148 (MQ=255)
                                                                                                                                    ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCatta                                       >  1:29975/1‑110 (MQ=255)
                                                                                                                                    ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCatta                                       <  2:29975/110‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGA  >  NC_000913/1608566‑1608844

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