Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,012,153 C→T 9.3% P51P (CCC→CCT)  pqiA → paraquat‑inducible, SoxRS‑regulated inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,012,1530CT 9.3% 134.0 / 4.9 43P51P (CCC→CCT) pqiAparaquat‑inducible, SoxRS‑regulated inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base T (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTA  >  NC_000913/1012015‑1012287
                                                                                                                                          |                                                                                                                                      
tcatGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCa                                                                                                                                       <  1:254384/140‑1 (MQ=255)
      cgcgAAGCACATCCTGGGCTCGCAGTGTGACCTGCAGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAACGCGGCATGTACCCGGTGTGGCACAACGTTAACCGTGGCGCGGGATGCCCCTCGGCAGCGTCCCACCGCCt                                                                                                                                 <  2:1097283/140‑1 (MQ=255)
        cgAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCTACCGCCTAt                                                                                                                               >  1:827551/1‑140 (MQ=255)
        cgAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCTACCGCCTAt                                                                                                                               >  1:827573/1‑140 (MQ=255)
             cacaTCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGtt                                                                                                                          <  1:1125591/140‑1 (MQ=255)
                    tgtgCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCa                                                                                                                   <  1:668982/140‑1 (MQ=255)
                    tgtgCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCACCTATGCGTTGGCTGCa                                                                                                                   <  1:224191/140‑1 (MQ=255)
                     gtgCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCAc                                                                                                                  <  1:124867/140‑1 (MQ=255)
                     gtgCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCAc                                                                                                                  <  1:124992/140‑1 (MQ=255)
                        cTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGt                                                                                                               >  2:367498/1‑140 (MQ=255)
                        cTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGt                                                                                                               >  1:296395/1‑140 (MQ=255)
                         tCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGtt                                                                                                              <  1:1045245/140‑1 (MQ=255)
                         tCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGtt                                                                                                              <  1:558284/140‑1 (MQ=255)
                          cGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTc                                                                                                             >  2:191495/1‑140 (MQ=255)
                             aGTGTGACATGCTGGTGGCGTTACCGCGCATTGAGCATGGTCAGAGAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGGCCCACCGCCGATGCGTTGGCTGCACTGGTCGtg                                                                                                          >  2:310257/1‑140 (MQ=255)
                              gtgtGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCAtgc                                                                                                         >  2:772151/1‑140 (MQ=255)
                                gtgACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCAtgctg                                                                                                       >  1:45364/1‑140 (MQ=255)
                                       gCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCAtgctgttgctgt                                                                                                >  2:1010269/1‑140 (MQ=255)
                                         tggtggCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTcc                                                                                              >  2:93164/1‑140 (MQ=255)
                                             ggCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACt                                                                                          >  2:299471/1‑140 (MQ=255)
                                                   aCCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTc                                                                                    >  1:6167/1‑140 (MQ=255)
                                                         ccTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTtg                                                                              >  1:1040017/1‑140 (MQ=255)
                                                                   ggTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCTACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGaa                                                                    <  2:827573/140‑1 (MQ=255)
                                                                   ggTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCTACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGaa                                                                    <  2:827551/140‑1 (MQ=255)
                                                                   ggTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGaa                                                                    <  1:510257/140‑1 (MQ=255)
                                                                   ggTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGaa                                                                    <  1:1083348/140‑1 (MQ=255)
                                                                         aaaGCGGCATGTCCCGGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGc                                                              >  2:557449/1‑140 (MQ=255)
                                                                              ggCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAg                                                         >  2:1034788/1‑140 (MQ=255)
                                                                                    tCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCa                                                   <  2:204083/140‑1 (MQ=255)
                                                                                      cccGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGt                                                 <  2:834274/140‑1 (MQ=255)
                                                                                      accGGTGTGGCACAACTTTAACCGTGGCGTGGTATGCCCCTCGGCAGCGTCCCACAGCCTATGCGTTTGCTGCACTGTTCATGCTGTTGCTGTCAAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGt                                                 <  2:834265/139‑1 (MQ=255)
                                                                                        cGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGa                                               <  1:329654/140‑1 (MQ=255)
                                                                                          gtgtGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGaaa                                             <  2:242792/140‑1 (MQ=255)
                                                                                              ggCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAattac                                         >  2:526913/1‑140 (MQ=255)
                                                                                              ggCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAattac                                         >  1:989885/1‑140 (MQ=255)
                                                                                                     cGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCtttt                                                                               <  2:573609/95‑1 (MQ=255)
                                                                                                     cGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCtttt                                                                               >  1:573609/1‑95 (MQ=255)
                                                                                                     cGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTg                                  >  2:1156341/1‑140 (MQ=255)
                                                                                                                        tGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCttt               <  2:450384/140‑1 (MQ=255)
                                                                                                                        tGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCttt               <  2:597833/140‑1 (MQ=255)
                                                                                                                               cGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTGTTGTGAATATGAACGTTGCAGGAGTTACAAGTGAAATTACATTACTGGAACTTCCCGGCGTGCTTTTTTCTGa        >  2:87837/1‑140 (MQ=255)
                                                                                                                                  cAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGa     <  1:187471/140‑1 (MQ=255)
                                                                                                                                     cGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTa  >  1:1033627/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCATGGTCAGAAAGCGGCATGTCCCCGGTGTGGCACAACGTTAACCGTGGCGTGGGATGCCCCTCGGCAGCGTCCCACCGCCTATGCGTTGGCTGCACTGTTCATGCTGTTGCTGTCCAACTTGTTTCCTTTTGTGAATATGAACGTTGCAGGAGTTACCAGTGAAATTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTA  >  NC_000913/1012015‑1012287

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