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,5250GT97.8% 137.5 / NA 47A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/17);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTTGC  >  NC_000913_3_vch_tpiA/2679393‑2679671
                                                                                                                                     |                                                                                                                                                  
acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGt                                                                                                                                    <  1:17631/149‑1 (MQ=255)
 caacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGAG‑TGCCGCTAAAATACCCCAGTCATGAGGGTTTAGTGAACGCCGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTGATCGCCGta                                                                                                                                   >  2:728032/1‑148 (MQ=255)
  aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCa                                                                                                                                  <  1:394752/149‑1 (MQ=255)
   acGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTTCAGTTACTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCCTCCGGCCAATCCGGATTAATCGCCATCaa                                                                                                                                 <  2:642527/149‑1 (MQ=255)
           ttCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgcgc                                                                                                                         >  2:337660/1‑149 (MQ=255)
               aaTTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATTGCCGTCAATTCAGCGCGAGc                                                                                                                     >  2:436869/1‑149 (MQ=255)
               aaTTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCATGGCGAGc                                                                                                                     >  2:233381/1‑149 (MQ=255)
               aaTTATTCTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                                                                                                     >  2:523705/1‑149 (MQ=255)
                  tattTTGTGATGAGCGGCAGCGCG‑TGCCGCTAAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCACCTACGGATGAATCGCCGTCAATTCAGCGCGAGCAGa                                                                                                                  <  2:181923/149‑1 (MQ=255)
                         tgatgaGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCatat                                                                                                           >  2:62698/1‑149 (MQ=255)
                          tttgaGCGGCA‑CGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTcgtgtt                                                                                                          >  1:739760/3‑144 (MQ=255)
                          gaggaGGGGCAGTGCG‑TGCCTCTGAAATATTACAGTTAATAGGGTTTAGTGAACGTGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt                                                                                                          <  2:505847/146‑1 (MQ=255)
                             gaGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGATTTCCGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg                                                                                                       <  1:585898/149‑1 (MQ=255)
                                cGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATCTTGCGCaa                                                                                                    >  2:766256/1‑149 (MQ=255)
                                 ggCAGTGGG‑TGCCGCTCAAATATTTCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                                   <  2:67894/149‑1 (MQ=255)
                                      cgcg‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc                                                                                              <  1:408948/149‑1 (MQ=255)
                                       gcg‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                                                                                                                  <  2:454754/128‑1 (MQ=255)
                                       gcg‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                                                                                                                  >  1:454754/1‑128 (MQ=255)
                                       gcg‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                             <  2:502707/149‑1 (MQ=255)
                                                   aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCgg                                                                                  >  1:547218/1‑149 (MQ=255)
                                                     atatTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAc                                                                                >  2:405589/1‑149 (MQ=255)
                                                        ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATGACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCAAAATGCGCAACCTCCCCGAGTACCGAGGACGcc                                                                             >  2:384092/1‑149 (MQ=255)
                                                            aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCa                                                                         >  1:240730/1‑149 (MQ=255)
                                                            aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCa                                                                         >  2:76257/1‑149 (MQ=255)
                                                                      gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGACTAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCACTTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCCAATCCCGCGAACGCCAGCATCCccatcct                                                               >  2:62038/1‑145 (MQ=255)
                                                                      gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATATTGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCTGCCTCCAGCATCGccagcac                                                               >  2:266234/1‑148 (MQ=255)
                                                                        ttAGTGAACGCGGTTAAAAATAAAAGGCATCACGATGCTGAAGAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCa                                                             >  1:624890/1‑149 (MQ=255)
                                                                              aaagcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAtt                                                       >  2:817052/4‑149 (MQ=255)
                                                                              aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCCtt                                                       >  2:543554/1‑149 (MQ=255)
                                                                              aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTACCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGACAtt                                                       >  1:346413/1‑149 (MQ=255)
                                                                               aCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGACGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATAGCCACCATCAGGCAttt                                                      >  2:318313/1‑149 (MQ=255)
                                                                                   ggTAAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCCTTCGGCCAATCCGGATTAATCGCCGTCAATTGAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTCCt                                                  >  1:127190/1‑149 (MQ=255)
                                                                                     ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCttt                                                <  1:62698/149‑1 (MQ=255)
                                                                                           aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgctg                                          >  2:800403/1‑149 (MQ=255)
                                                                                                  gCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgttgt                                   >  2:96058/1‑149 (MQ=255)
                                                                                                       aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATc                              >  1:627846/1‑149 (MQ=255)
                                                                                                       aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATc                              >  1:627845/1‑149 (MQ=255)
                                                                                                          gtggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGcc                           <  1:96058/149‑1 (MQ=255)
                                                                                                          gtggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGcc                           <  1:438551/149‑1 (MQ=255)
                                                                                                               gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACa                      <  1:405589/149‑1 (MQ=255)
                                                                                                               gAATAATGCCGCTTTCGGCCAATCCGGACTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACa                      <  1:428340/149‑1 (MQ=255)
                                                                                                                         gCTTTCGGCCAACCCGGATTAATCGCCGTCAATTCAGCTCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa            <  2:815706/149‑1 (MQ=255)
                                                                                                                             tCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCa        <  1:817052/149‑1 (MQ=255)
                                                                                                                              cGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAg       <  2:240730/149‑1 (MQ=255)
                                                                                                                                gCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGGGTCCCGCGGACGCCAGCATCGCCAGCATCAGGCATTTTTCTTTATCCTGCTGGTGTAAATCGCCATCCTCCTACAACAAGTCACGt     >  2:40582/1‑149 (MQ=255)
                                                                                                                                gCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGTCAAcaagtcaagt     >  2:208343/1‑149 (MQ=255)
                                                                                                                                   aaGCAGGATTAATCGCCGGCAATTGAGCGCGGGCAGACTCATATTGCGCAACTTTCCCGAGTCCGGCGGACGCCAGCATCTCCACCATCAGGCATTTTTCTTTCTGCTGCTGTTGTAAATCGCCAGACAGCGACAACAAGTCAGGTTGc  <  2:814550/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
ACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTTGC  >  NC_000913_3_vch_tpiA/2679393‑2679671

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