Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,497,170 A→T 14.8% M153L (ATG→TTG)  ydcL → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,497,1700AT14.8% 62.7 / 6.7 27M153L (ATG→TTG) ydcLlipoprotein
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base T (2/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-02

ACAGCATCGTCTGGAACCCAATCACTTATTATCCGGTACCGAAACCGTCGACCCAGGTAGGGCAGAAAGTTCTGGATAAAATTTTGAACTATACCAACACCGAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCA  >  CP009273/1496910‑1497438
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                            
cccgCATCGTCTGGAAACAAACCACTTATTGTCCGGTACCTGAACCGTCGACCCAGGTCAGTAAGAAAGTTCTAGATAAAATTTTGAACTATACCAAGACCGAAATGAAAGAAGTGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCa                                                                                                                                                                                                                                                                 <  2:185146/271‑1 (MQ=255)
                   aaTCACTTATTATCCGGTACCGAAACCGTCGACCCAGGTAGGGCAGAAAGTTCTGGATAAAATTTTGAACTATACCAACACCGAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCATGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGATGCGGGGACGCAAATCGCTACAGACCACCGTACCATGGATACTCTcc                                                                                                                                                                                                                                             >  1:648513/1‑275 (MQ=255)
                                                 gACCCAGGTAGGGCAGAAAGTTCTGGATAAAATTTTGAACTATACCAACACCGAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCCAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGATACAGGCAACCTCACCATGGATACTCGCCTCCATTTCGACGGTGNGCTGATTGATGc                                                                                                                                                                                                                 >  2:531533/1‑273 (MQ=255)
                                                     cAGGTAGGGCAGAAAGTTCTGGATAAAATTTTGAACTATACCAACACCGAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGNTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCATACCATGGATACTCGTCTCTATTTTGAAGGTGAGCTGCTTGATg                                                                                                                                                                                                                  <  1:531542/268‑1 (MQ=255)
                                                                                                     gAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGACGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGGTACTCGCCTCTTTTTTGAAGGTGCGCAGATTAAGGCAGAGACTAATAAACAGGTTAACAAAGGGGTTCGACAGGGCGAAGGTGAAgccgt                                                                                                                                                           >  2:202081/1‑271 (MQ=255)
                                                                                                                                       tGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAATTGGCTACAGGCCACCGTACCATGGATACTc                                                                                                                                                                                                                                                >  1:621467/1‑156 (MQ=255)
                                                                                                                                           tACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAATTGGCTACAGGCCACCGTACCATGGATACTCGCct                                                                                                                                                                                                                                            <  2:621453/156‑1 (MQ=255)
                                                                                                                                                                                          gtgtAGATACCAGCAAAGTAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGCGAGCTCATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGt                                                                       <  1:658353/274‑1 (MQ=255)
                                                                                                                                                                                          gtgtAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTCTTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGACAATGGCTTACGAAAATATAAAAAAAGTTATTGATGACATGCCGACCGATGCCACCATGTTTGACGt                                                                       >  2:658339/1‑274 (MQ=255)
                                                                                                                                                                                               gATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTaa                                                                                                                                                                 >  2:199217/1‑179 (MQ=255)
                                                                                                                                                                                                   ccAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGAc                                                                                                                                                             <  1:199219/179‑1 (MQ=255)
                                                                                                                                                                                                    cAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAAAAAGTACTTGATGACATGGCGACCGATGCCACCATATTTGACGTTAACAAAAag                                                             >  1:348490/1‑274 (MQ=255)
                                                                                                                                                                                                           gaagGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAAt                                                                                                                                                                                                                                                                             <  1:273511/59‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAATTGGCTACAGGCCACCGTACCATg                                                                                                                                                                                                                                                       >  2:728897/1‑59 (MQ=255)
                                                                                                                                                                                                                                 tGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATg                                                                                                                                                                                                                                                       >  2:731127/1‑59 (MQ=255)
                                                                                                                                                                                                                                     tCCTGTTGCATTAGTGGTTGCGGGGACGCAATTGGCTACAGGCCACCGTACCATGGATa                                                                                                                                                                                                                                                   <  1:728911/59‑1 (MQ=255)
                                                                                                                                                                                                                                     tCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATa                                                                                                                                                                                                                                                   <  1:731141/59‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAg                                                                                                                                                                                                               >  2:96835/1‑89 (MQ=255)
                                                                                                                                                                                                                                              aTTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTAt                                                                                                           >  1:593623/1‑186 (MQ=255)
                                                                                                                                                                                                                                               ttAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGAc                                                                                                                                                                                                           <  1:96836/89‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACtaataa                                                                                                                                                                                                     >  2:714163/1‑92 (MQ=255)
                                                                                                                                                                                                                                                  gTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATtgat                                                                                                       <  2:593610/186‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCg                                                                                                                                                                                                 <  1:714177/92‑1 (MQ=255)
                                                                                                                                                                                                                                                             gACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAgg                                                                                                                                                                           >  1:76424/1‑107 (MQ=255)
                                                                                                                                                                                                                                                             gACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAgg                                                                                                                                                                           >  1:588230/1‑107 (MQ=255)
                                                                                                                                                                                                                                                              caacgaaTAGCTCCAGGCCAGCCGACCATGGATACTGGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa  <  2:44742/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGa                                                                                                                                                                       <  2:588217/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTAGTTTGAAGGTGATCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGa                                                                                                                                                                       <  2:76423/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                            
ACAGCATCGTCTGGAACCCAATCACTTATTATCCGGTACCGAAACCGTCGACCCAGGTAGGGCAGAAAGTTCTGGATAAAATTTTGAACTATACCAACACCGAAATGAAAGAAGCGATAGCGCAGCGTAAACCACTGGTTACCACCGCTGGGCCGCGTAGTCTGATTTTCCGTGGGGCCATTACCGGTGTAGATACCAGCAAAGAAGGGCTGCAATTCTATGAAGTGGTTCCTGTTGCATTAGTGGTTGCGGGGACGCAAATGGCTACAGGCCACCGTACCATGGATACTCGCCTCTATTTTGAAGGTGAGCTGATTGATGCAGCGACTAATAAACCGGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCA  >  CP009273/1496910‑1497438

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: