Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,389,186 T→C N51S (AAC→AGC)  aaeX ← DUF1656 domain‑containing protein AaeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,389,1860TC97.3% 104.4 / NA 37N51S (AAC→AGC) aaeXDUF1656 domain‑containing protein AaeX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/14);  minor base A (1/0);  total (23/14)
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.98e-01

TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACGATAACGGG  >  NC_000913/3389047‑3389325
                                                                                                                                           |                                                                                                                                           
ttAAAAATTGCGATAAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCCCGAAAGACCACTCGCGATATTAAATAAAACAAGCAGCAACAGAGCGCGGTGCTGa                                                                                                                                          >  1:126756/1‑143 (MQ=255)
   aaaaTTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTGTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCGATAGAGCGCGGTGCTGgaca                                                                                                                                       >  1:100842/1‑143 (MQ=255)
   aaaaTTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaaca                                                                                                                                       >  2:82159/1‑143 (MQ=255)
   aaaaTTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAACGCGGTGCGGaaca                                                                                                                                       >  2:218189/1‑143 (MQ=255)
       ttGTGATGAAGTCCAGAAGGACTAATACGACCGTGATGGCCGTGCGGGAGAAGTTTCTTATTAGTGTTTTCACTTCAACCCCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGc                                                                                                                                   <  2:760660/143‑1 (MQ=255)
       ttGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGc                                                                                                                                   >  2:627851/1‑143 (MQ=255)
                  gCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTCGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGAGCGGTGCTGAACAACGCCGGATGCCGGa                                                                                                                        >  2:261381/1‑143 (MQ=255)
                   ccAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGAc                                                                                                                       >  2:435644/1‑143 (MQ=255)
                    cAGAATGACTAAGACTACCGTGATGGCCGTACTGGAGCATCTGCTTATTACTCTTTTCACTTCAACCTCAGACGACCAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACa                                                                                                                      <  2:797833/143‑1 (MQ=255)
                       aaTGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACACCGGATGCCAGACAAAg                                                                                                                   <  1:548768/143‑1 (MQ=255)
                            ctactccGACCGTGACGGCCGTACGGGAGAATTTTCTTAGTAGTGTTTTCACTTCAACCTCAAACGAACACGCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTa                                                                                                              <  1:616150/137‑1 (MQ=255)
                                tACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATa                                                                                                          <  1:399316/143‑1 (MQ=255)
                                 aCGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGGACAACACCGGATGCCAGACAAAGTCGTAGATAc                                                                                                         >  1:283717/1‑143 (MQ=255)
                                 aCGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATAc                                                                                                         >  1:856084/1‑143 (MQ=255)
                                 aCGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGAAGCCAGACAAAGTCGGAGATAc                                                                                                         >  1:218812/1‑143 (MQ=255)
                                         gATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGGCGTAGATACCTGTTGGc                                                                                                 >  1:330672/1‑143 (MQ=255)
                                            ggCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGcaca                                                                                              >  1:519070/1‑143 (MQ=255)
                                             gCCCTACGGGAGCGTGTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGGGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACaa                                                                                             <  2:435484/143‑1 (MQ=255)
                                                    gggAGAATTTTCTTATTAGTGTTTTCACTTCAACCGCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCg                                                                                      <  1:161630/143‑1 (MQ=255)
                                                                  attaGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCCATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCCCAAGTACCCGGCGCACCAGCCAGa                                                                        >  1:100246/1‑143 (MQ=255)
                                                                     attgtTTTCACTTGCACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGaaaa                                                                     <  2:269992/141‑1 (MQ=255)
                                                                       tgtTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAACGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATc                                                                   >  1:435589/1‑143 (MQ=255)
                                                                         ttttCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGc                                                                 <  2:481149/143‑1 (MQ=255)
                                                                              aCTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACACGCAGCAATAGAGCGCGGAGCTGAACAACGCCGGATGCCAGACAAAGGCGTAGATACCTGTTGGCTCAAGTACCCGGCGCACCAGCCACAAAATCggcgggg                                                            >  2:58236/1‑137 (MQ=255)
                                                                                ttCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGAt                                                          >  2:920415/1‑143 (MQ=255)
                                                                                                   tCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAata                                       >  1:874859/1‑143 (MQ=255)
                                                                                                   tCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGACGAACAACGACGGAAGCAAGACGAAGGCATCGATACCTGTTGGCACAAGTACCAGGCGCACCACCCAGCAAAACGCCAGGGATAAAACCAATTGAACAAata                                       >  2:128604/1‑143 (MQ=255)
                                                                                                        atatCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCAGAATGCCAGACAAAGTCGTAGATACCTGTTGGCACACGTACCCGGCGCACCAGCCAGCAAATCGCCCGTAATAAAAGCACTTCAAAACATATCGGt                                  >  2:755269/1‑143 (MQ=255)
                                                                                                              aaTAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACCAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTAAAAAAATAACGAAGGGAAg                            >  2:293045/1‑143 (MQ=255)
                                                                                                                  aaaCAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGACGTAGATACCTGTTGGCACAAGTACCCGACGCACCACCCAGAAAATCGCCAGTGATAACAGCAAGTCCAAAAATATCGGTGGCAAGGACa                        >  2:428870/1‑143 (MQ=255)
                                                                                                                   gaCAAGCAGTAAGAGAGCGCGGTGCTGAACAGCGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAg                       <  2:214928/142‑1 (MQ=255)
                                                                                                                   acccagccgcAATAGCGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATCCCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAg                       <  1:218189/139‑1 (MQ=255)
                                                                                                                   aaCAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGGGATAAAAGCCATTCAAAAAATAACGGCGGGACGGACAg                       >  2:472151/1‑143 (MQ=255)
                                                                                                                         cagcaATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCaaa                 <  1:627851/143‑1 (MQ=255)
                                                                                                                            caATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCAGAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATCAGGGCGAAATAGAGAGTATCGGAGCGAAGAACAGCCCAAacac              >  2:190899/1‑143 (MQ=255)
                                                                                                                               tcgagCGCGGTGCTGAACAACGCCGGATGACAGACAAACTCGTAGATATCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAAcaccac           <  2:837348/141‑1 (MQ=255)
                                                                                                                               tAGAGCGCGGTGCTGATCAACGCCGGATGCCAGACAAAGTCGTAGATAGCTGTGGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAAcaccac           <  2:750540/143‑1 (MQ=255)
                                                                                                                                        gtgcTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCGGCCCCCTGAAAAACGCGCGTGGTGAGGTAACGACAAAAAAGACCGGTGGGCAAGACAGCCCCAACACCACGAAAACggg  >  2:901728/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACGATAACGGG  >  NC_000913/3389047‑3389325

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