Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA2,679,5250GT95.8% 72.6 / NA 24A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (15/8);  minor base C (0/1);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTT  >  NC_000913_3_vch_tpiA/2679390‑2679669
                                                                                                                                       |                                                                                                                                                
ttaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGc                                                                                                                                     >  2:345036/1‑149 (MQ=255)
  acgtacGGGTTATTTTAGGTATTTTGTGAGGAGCGGCAGCGCGTGACGCTCAAATAGTAGAGTGAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                                                                                                   <  2:153682/145‑1 (MQ=255)
   acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGt                                                                                                                                  <  1:441406/149‑1 (MQ=255)
                cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                                     <  1:38487/149‑1 (MQ=255)
                cTAAATAGTTTTTGATGAGCGGCAGCGCATGCCGCTCAAATATTGCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCAGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                                     >  2:177687/1‑149 (MQ=255)
                                     gCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCCATCCGGATTAATCGCGGTCAGTTCAGCGCGAGCAGACTCATCTTGCGCGAcc                                                                                                >  1:129921/1‑149 (MQ=255)
                                     gCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGCCTCATATTGCGCAAcc                                                                                                >  1:129918/1‑149 (MQ=255)
                                      cAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCt                                                                                               >  1:381336/1‑149 (MQ=255)
                                       aGCGCGTGCCGCTAAAATATTAAAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTGCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc                                                                                              <  1:43690/149‑1 (MQ=255)
                                          gcgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                           <  1:210830/149‑1 (MQ=255)
                                               ccGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAACGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTc                                                                                      >  1:12118/1‑149 (MQ=255)
                                                                   aTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGcc                                                                                                                                    >  1:20175/1‑83 (MQ=255)
                                                                   aTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGcc                                                                                                                                    <  2:20175/83‑1 (MQ=255)
                                                                    tAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGCCTCATATTGCGCAACCTCCCCGAGTACCGAGGACGCCCGCATCGcca                                                                 >  2:222728/1‑149 (MQ=255)
                                                                        gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCAt                                                             >  1:322474/1‑149 (MQ=255)
                                                                                      gTTAAAAATAAAAGGCATCACGGCGGTGAATAATGCCGCTTCCGGCCAACCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCtt                                               <  1:58218/149‑1 (MQ=255)
                                                                                        tAAAGATAAAGGGCTTCACGGTGGTGAATTATGTGGCTTTCGGGCAATCGGGATTAATCGTCGTCAATTTAGGGCGAGCAGGTTCGTATTGCGCAGGGTCGCCGTGTCCCGCGGACGCCAGCATCGGCACCATCAGGCATTTTTCtttg                                             <  2:89910/149‑2 (MQ=255)
                                                                                          aaaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtg                                           >  2:480385/1‑149 (MQ=255)
                                                                                                   ggCATCACGGTGGTGAATAATGCCGCTTTCAGCCAATCCGGACTAAGCGCCGTCAATTCAGCGCGAGCAGACTCCTATTGCGCAACCTCCCCGAGTCCCGCGGACGCCCGCATCGCCACCATCAGGCATTTTCCTCTATGCTGCtgttg                                  >  2:349974/1‑149 (MQ=255)
                                                                                                             tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCa                        <  1:345036/149‑1 (MQ=255)
                                                                                                               gtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATa                      >  1:349401/1‑149 (MQ=255)
                                                                                                                  aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCTCCATACAg                   >  1:153035/1‑149 (MQ=255)
                                                                                                                         ccGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGacaac            >  1:270885/1‑149 (MQ=255)
                                                                                                                                   ccAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGGCGCCTGCATCTCCACCATCAGGCATTTTTTTTTATGCTGCTGTTGCAAGTCGCCATTCAGCGGCACAAAGTCAGGtt  >  2:3/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                
TTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTT  >  NC_000913_3_vch_tpiA/2679390‑2679669

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