Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,111,596 G→A L179L (CTG→TTG)  yiiQ ← DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,111,5960GA97.0% 100.0 / ‑4.8 33L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/13);  minor base T (0/1);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGC  >  NC_000913_3_pae_tpiA/4111451‑4111736
                                                                                                                                                 |                                                                                                                                            
aaaGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACTGTATCGGGGTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAccc                                                                                                                                           <  2:458790/149‑1 (MQ=255)
          ggCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCg                                                                                                                                 <  1:520950/149‑1 (MQ=255)
                 tagCAGGGCAGTAGGTCTGTCAAACTAGATAGAGATTCATCGGCACAAAGGCTGTGCTTGTTGTCATTGATTCAAACCTTCAAGCGATTCAGATAGCGCCTGCTTCATCGTTTCAAAAGCCAAGGTCATTCCCTTTTCGCCGTTTTCcg                                                                                                                          <  2:109908/147‑1 (MQ=255)
                          agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACAt                                                                                                                 >  1:370418/1‑149 (MQ=255)
                          agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCAAAGATCAACCCCATATCGCCATTATCCGCGACAACAt                                                                                                                 >  1:10274/1‑149 (MQ=255)
                                 cTGTGAAACAGTATAGAGATTCATCGGAACAAAGGCTATGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                          <  1:601023/149‑1 (MQ=255)
                                    tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTTTTTTTTTTTCAAACCTTCAAGCGCTTCAGATAGCACAAGCTTTATCTGTTCAACAGCGAATGTCAACCCATTTTCGCTGTTCTCCGACACACCATCACGCaatct                                                                                                       >  2:192220/1‑145 (MQ=255)
                                     gAAACAGTATAGAGATTCATCGGCACAAAGGCTATGGTGTTGGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                                      >  1:316750/1‑149 (MQ=255)
                                         cAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTTTTCAAACCTTCAAGCGATTCAGCTAGCGCCAGCTTAATCGGTTCAACACCAAAGTTCAACCCCTTTTCGCCGTTGTCCGCGACAACATCACGCAACGCACCtt                                                                                                  >  2:88530/1‑149 (MQ=255)
                                          aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                 <  1:238148/149‑1 (MQ=255)
                                              taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                         <  2:333620/121‑1 (MQ=255)
                                              taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                         >  1:333620/1‑121 (MQ=255)
                                                           gCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                >  2:452472/1‑149 (MQ=255)
                                                           gCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                >  2:233715/1‑149 (MQ=255)
                                                                aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                           <  1:415996/149‑1 (MQ=255)
                                                                 cgcgttggctttTTGTCATTTATTGAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                          <  1:240151/142‑1 (MQ=255)
                                                                   gtgttgctttTTGTCATTTATTCAAACCTACGAGCGATTCAGATGGCGCCAACTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgttt                                                                        <  2:370418/146‑1 (MQ=255)
                                                                   gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTTTAATAACatttgttt                                                                        >  1:299170/1‑149 (MQ=255)
                                                                     tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                     >  1:623040/1‑86 (MQ=255)
                                                                     tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                     <  2:623040/86‑1 (MQ=255)
                                                                                tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                           >  2:388650/1‑149 (MQ=255)
                                                                                          aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTATGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGGCTg                                                 >  1:410258/1‑149 (MQ=255)
                                                                                                      gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                     >  2:446673/1‑149 (MQ=255)
                                                                                                       aTTCAGATAGCGCCAGTTTAATCGGTTTAACAGCGAAGGTCAACCCTTTTTCGCCGTTTTCCGCTACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                    <  1:184772/149‑1 (MQ=255)
                                                                                                       aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                    <  1:446673/149‑1 (MQ=255)
                                                                                                           aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGGttt                                >  1:197888/1‑149 (MQ=255)
                                                                                                               aGCGCCAGCTTAATCGGTTCAACAGAGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTATTAGCGACTGCAGTTTTTTCTGACTTTgtgt                            >  2:526401/1‑149 (MQ=255)
                                                                                                                 cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGgtt                          >  1:236507/1‑149 (MQ=255)
                                                                                                                     aGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCACGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                      >  1:593593/1‑149 (MQ=255)
                                                                                                                               gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg            <  1:330416/149‑1 (MQ=255)
                                                                                                                               gTTCAACAGCGAAGGTCAACCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg            <  1:68183/149‑1 (MQ=255)
                                                                                                                                    aCAGCGAAGGTCAACCCCTTTTCGCCGTTATCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCACCGCCCTTCGCAGCTGCTTCAGTTTTTTCTGGCTTTGTGCTTTGTTCCTTAAGGGGATGAg       >  2:394471/1‑149 (MQ=255)
                                                                                                                                         gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTTATTAATGGGGTGAAAGTGc  >  1:439175/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGC  >  NC_000913_3_pae_tpiA/4111451‑4111736

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: ‑