Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_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_pae_tpiA2,679,5250GT95.8% 61.6 / NA 24A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (13/10);  minor base C (0/1);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGT  >  NC_000913_3_pae_tpiA/2679384‑2679663
                                                                                                                                             |                                                                                                                                          
tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAtt                                                                                                                                     >  1:178483/1‑149 (MQ=255)
 ttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTa                                                                                                                                    >  2:376666/1‑149 (MQ=255)
  tttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                                                                                                   >  2:276720/1‑149 (MQ=255)
      ttaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGc                                                                                                                               <  1:287339/149‑1 (MQ=255)
       taACAACGGCTTATTCTAATTATTTTGTGATGAGCTGCAGCGCGTGCCGCTCTAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGcc                                                                                                                              <  1:1180/149‑1 (MQ=255)
         acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGACGCTTTCGGCCACTCCGGATTAATCGCCGt                                                                                                                            >  2:315011/1‑149 (MQ=255)
                        aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                                                                                             >  1:508225/1‑149 (MQ=255)
                                          ggCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAAAAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                           >  1:363255/1‑149 (MQ=255)
                                                       gCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGGGCGCAGGCTCATATTTAGCAACCTCCCCGAgcacc                                                                              >  1:191055/1‑145 (MQ=255)
                                                             atatTACAGTTAATGGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAc                                                                        <  1:356995/149‑1 (MQ=255)
                                                                ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                     >  2:105610/1‑149 (MQ=255)
                                                                                 tAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAg                                                    <  1:320579/149‑1 (MQ=255)
                                                                                      aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGCCGCCAGCATCGCCCCCATCAGGCAtt                                               >  1:294052/1‑149 (MQ=255)
                                                                                         gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAttttt                                            <  1:185528/149‑1 (MQ=255)
                                                                                           ggTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                          >  2:63221/1‑149 (MQ=255)
                                                                                             ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATc                                                      >  1:423048/1‑135 (MQ=255)
                                                                                             ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATc                                                      <  2:423048/135‑1 (MQ=255)
                                                                                                 aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCACGCATGTTTCTTTAtgc                                    >  1:267832/1‑149 (MQ=255)
                                                                                                               aCGGTGGTGAATAATGCCGCTTTCGGCCAACGCTGACTAATCGCCGGCAAGTCAGCGCGAGCCGATTAACCCTACGCAACGCCCCCTATTCCCGCGGACGCCAGCACCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATc                      <  2:215372/149‑1 (MQ=255)
                                                                                                                   tggggAGTAAGGCGGCTTGCGGGTAATCCGGATGAATCGCCGTCAATTCAGCGTCAGCAGACTCATCTTGCGCAACCCCCCCGAGTCCCGCGGACGCCCGCATCGCGACCATCAGGCAGTTTTCTTTATGCTGCTGTTGTAAATCGCCa                  <  2:396781/149‑1 (MQ=255)
                                                                                                                     gtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATa                <  1:391632/149‑1 (MQ=255)
                                                                                                                           aatGGCGCTTTCGGCCAATGCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGa          <  1:129391/149‑1 (MQ=255)
                                                                                                                             tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGGGCAACCTCGCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTACGCACCTGATGTAAATCGCCATACAGCGtca        >  1:10529/1‑149 (MQ=255)
                                                                                                                                   tgtccgCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAAGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGt  <  2:346930/144‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGT  >  NC_000913_3_pae_tpiA/2679384‑2679663

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