Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,795,396 C→A Q382K (CAG→AAG)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,795,3960CA96.7% 92.7 / NA 30Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (13/16);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795268‑2795539
                                                                                                                                |                                                                                                                                               
tACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGCTAAACCAGTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGCTGCGGATTCGCAAAGtt                                                                                                                             >  2:94045/1‑149 (MQ=255)
   tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTCCGCGCCTGCGAAGGTATTTAAATTTCTGATCGGCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGCTTCAAAGTTTCGGAGGCGCAAAATTCTg                                                                                                                          >  2:81411/1‑149 (MQ=255)
   tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                          >  1:156403/1‑149 (MQ=255)
     tcCACCGGCGAGGCATTCTTAACGGTGGTGGGAAACTATTACGGGCCTACGAAGGTGTTTAAATATCTGAACTCCAGCTCCGGCGCTATCGCCCTCCTGGTTTATATAGTCATGGCGTGTTGAAAGTTTCGGTAGCTCAAAACTCTGCg                                                                                                                        >  2:75038/1‑149 (MQ=255)
            gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                                 <  1:274465/149‑1 (MQ=255)
                          aCGGTGGGGGTGAACTATTCCGCGCCTGCGAATGGATTCAAATTTCAGATCGACAGCTCCGGCGCTAGCGCCCGGCTGGGTTAGTTAGTCATCGCCGTTTCAAAGGTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAGTCg                                                                                                   <  2:181189/149‑1 (MQ=255)
                           cGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                  <  2:156403/149‑1 (MQ=255)
                           cGGTGGTGGTGAACTATTACGCGCCTGCGAAACTATTTGTGGTCTTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                  <  1:377441/149‑1 (MQ=255)
                                       aCTATTACGCGCCTTCGAAGGGATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                      <  2:126138/149‑1 (MQ=255)
                                        cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     <  1:494946/149‑1 (MQ=255)
                                                  cttGCGAAGGCATGTAAATGTTTCATCGATGGCGCCGGTGCTTCCGCCCTGCTTGTTTATTTAGTCATCGCCGTTTCAAATTTGCCGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                           <  2:2193/147‑1 (MQ=255)
                                                    tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     <  2:546866/137‑1 (MQ=255)
                                                    tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     >  1:546866/1‑137 (MQ=255)
                                                       gAAGGTATTTATATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTCCGCCTGTTGCGCTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCa                                                                      <  1:128393/149‑1 (MQ=255)
                                                       gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTACTGGTTTATTTAGTCATCGCCGATTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGACGGAAACGAAATGCACTTCCCCAAGTGCCTATATCCGTGGCtac                                                                      >  2:433025/1‑147 (MQ=255)
                                                         aGGTATTTAAATCTCTGATTGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGCCacaca                                                                    <  2:421852/149‑2 (MQ=255)
                                                                tcaaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTCAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAg                                                             <  1:205728/147‑1 (MQ=255)
                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGCTGCGCAAAATTCTGAGAGCAGAAGGAAGCGAAATTCCGTTGCGAATGTGGCTTTATCCGCGGCTCACCTGGCTAGTCAt                                                         >  2:128943/1‑149 (MQ=255)
                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaag                                                                                                                <  1:504854/94‑1 (MQ=255)
                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaag                                                                                                                >  2:504854/1‑94 (MQ=255)
                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                         >  2:143866/1‑149 (MQ=255)
                                                                          aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCTTTGCTTAGCTGGCGAGTAATTTGCtt                                                   >  1:191544/1‑149 (MQ=255)
                                                                                      ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                       <  1:292641/149‑1 (MQ=255)
                                                                                      ggCGCTATCGCCCTGCTGGTTTATTTAGTCAACGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                       <  1:61626/149‑1 (MQ=255)
                                                                                           tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCCCCTGCGAATGGAGCAAAACGCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  <  1:268982/149‑1 (MQ=255)
                                                                                                   tgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTGATTACCTTTGTGTTGGTAGTGATg                          >  1:364422/1‑149 (MQ=255)
                                                                                                        gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     >  1:452877/1‑149 (MQ=255)
                                                                                                               taGGCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg              <  1:471928/149‑1 (MQ=255)
                                                                                                                     tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa        >  1:355392/1‑149 (MQ=255)
                                                                                                                           tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTa  >  2:66657/1‑149 (MQ=255)
                                                                                                                                |                                                                                                                                               
TACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795268‑2795539

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: