Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,263,710 C→T L436L (CTG→CTA)  tdcC ← threonine/serine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,7100CT96.8% 85.6 / ‑4.0 31L436L (CTG→CTA) tdcCthreonine/serine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/11);  total (20/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 = 7.71e-01

AATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTCAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGT  >  NC_000913/3263470‑3263952
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
aaTGCACCTTCTTAGCTGTTGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACTGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATGCATTCATCTCTTTGCTCATCCGGAGGTACGGGTGAAAACATTTTGTATGCGATGTGTAGGATGGTca                                                                                                                                                                                                                                            <  2:179182/251‑1 (MQ=255)
aaTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTc                                                                                                                                                                                                                                             >  1:198603/1‑250 (MQ=255)
           gTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑Tc                                                                                                                                                                                                                                >  1:31634/1‑251 (MQ=255)
                                             ccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGt                                                                                                                                                                                             >  1:338129/1‑251 (MQ=255)
                                                     ggATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGGCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGCCCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTcacc                                                                                                                                                                                     >  2:118884/1‑248 (MQ=255)
                                                                 aTTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCt                                                                                                                                                                         >  1:92504/1‑251 (MQ=255)
                                                                                     cGGCAATACCTGACATTAATACTTCACAGTCGATGGCATCGAGAACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTAATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCgg                                                                                                                                                     >  1:371892/1‑251 (MQ=255)
                                                                                                                aGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAAt                                                                                                                          >  1:234097/1‑251 (MQ=255)
                                                                                                                 gTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATc                                                                                                                         <  1:96208/251‑1 (MQ=255)
                                                                                                                          aTCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAACCAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGccc                                                                                                                >  1:382825/1‑251 (MQ=255)
                                                                                                                               gCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGc                                                                                                           <  1:20022/251‑1 (MQ=255)
                                                                                                                                  cGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACTGCGCTTTACGGATGGCATACATCGGCAACAGGCACATCAGGGATGCGATAATCGGTGCGCCCATGGCTTc                                                                                                        <  1:406598/251‑1 (MQ=255)
                                                                                                                                          ttAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTAATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGt                                                                                                >  1:423643/1‑251 (MQ=255)
                                                                                                                                               cGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAgg                                                                                           <  1:89182/251‑1 (MQ=255)
                                                                                                                                                 aaGACGAACCACAGTTAATAACCCAAACAACCGGAAATTCATTCAACTCTTTTCTCATCCTGAGATACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGAAGGCGCTTTACGGATGGCATACATCGGCAACAGGCCCAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGtacaggat                                                                                         >  2:248446/1‑244 (MQ=255)
                                                                                                                                                  aGACGAACCACAGTTAATAACCAAAACAACCGACAATACATTCACATCTTTTCTCAACCTGAGTTACGGATTAAAACAGTCTGTATACGATGTTTAGGATGCTAAGCA‑ACCAAATCACGGTAAACAAA‑AAGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATg                                                                                        >  1:452253/1‑251 (MQ=255)
                                                                                                                                                                       cAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCATGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGcaaca                                                                   <  1:400372/251‑1 (MQ=255)
                                                                                                                                                                                        cattcatCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCat                                                  >  1:156121/1‑251 (MQ=255)
                                                                                                                                                                                            catCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATCCGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACAACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGCTAATCGGTGCGCCCATGGCTGCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACTACACAGGTGGAGCCCatgatg                                              >  1:436500/1‑251 (MQ=255)
                                                                                                                                                                                                 ctTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCAtgatga                                             <  1:152726/247‑1 (MQ=255)
                                                                                                                                                                                                        cATCCTGAGTTACGGATTAAAGCAGTTGGTATACGATGGTTAGGATGGTCAGCAGGCCAA‑TCACGGTAA‑CAAACACGTTATCCAGGCGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCGCATGGCTTCAATCAGGTCAAGGATGGTCGTGTTGGCGTAGGCAACAACCCAGGTGGTGCCCATGATGAAGATCATg                                     <  2:47594/248‑1 (MQ=255)
                                                                                                                                                                                                            cTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAgg                                                                                           >  1:233828/1‑190 (MQ=255)
                                                                                                                                                                                                            cTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAgg                                                                                           <  2:233828/190‑1 (MQ=255)
                                                                                                                                                                                                            cTGAGTTACGGATTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGAt                                >  1:54420/1‑249 (MQ=255)
                                                                                                                                                                                                                       aTTAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCACCAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGttt                    >  1:253183/1‑250 (MQ=255)
                                                                                                                                                                                                                        ttAAAACAGTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTAc                  >  1:16657/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCAGTTTGTATACGATGGTTATGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACGCGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCGTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCGTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGc             <  2:33790/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTTTGTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGac          >  1:87249/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTATACGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCAGGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCgtctct       <  2:338129/250‑6 (MQ=255)
                                                                                                                                                                                                                                       taCGATGTTTAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGAACACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATAGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACttt     >  1:359802/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aCGATGTTCAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATCGTGTTCAGTTTACCCAGCGACACTTTAGt  >  2:267004/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
AATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTTCTCATCCTGAGTTACGGATTAAAACAGTTTGTATACGATGTTCAGGATGGTCAGCAGACCAA‑TCACGGTAA‑CAAACACGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGT  >  NC_000913/3263470‑3263952

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