Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,022,234 C→A 12.1% W205L (TGG→TTG)  yedI ← DUF808 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,022,2340CA12.1% 104.9 / 2.1 33W205L (TGG→TTG) yedIDUF808 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (16/13);  new base A (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TCCCACCGCCGACGAGGAACATCGCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTCAAACTTCAGCGGATCCTGCGC  >  CP009273/2022109‑2022451
                                                                                                                             |                                                                                                                                                                                                                         
tCCCACCGCCGACGAGGAACATCGCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGtcatca                                                                                                                                                                                                          >  1:46139/1‑143 (MQ=255)
    aCCGCCGACGAGGAACATCGCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCt                                                                                                                                                                                                      <  2:46137/143‑1 (MQ=255)
                     tCGCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTa                                                                                                                                                                                                    >  2:347338/1‑128 (MQ=255)
                      cGCCAGCGTGGCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACTTGATTAAGCAACGGCGCTTCGGCCACAATcccc                                                                                              >  1:489647/1‑229 (MQ=255)
                         cAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCACGCAATGCCCCAGGCCATCAATCTTAACGa                                                                                                                                                                                                <  1:347350/128‑1 (MQ=255)
                          aGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                                          <  2:489634/229‑1 (MQ=255)
                                ccGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATa                                                                                                                                                                                                                       >  1:282733/1‑98 (MQ=255)
                                    cAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATAcccc                                                                                                                                                                                                                   <  2:282725/98‑1 (MQ=255)
                                                   tcatcaGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCAAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAg                                                                                                                                                        >  1:177433/1‑142 (MQ=255)
                                                       caGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCAAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGcgcg                                                                                                                                                    <  2:177427/142‑1 (MQ=255)
                                                         gCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAAt                                                                                                                                                                                                                        >  1:39867/1‑72 (MQ=255)
                                                             gggCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATAccc                                                                                                                                                                                                                    <  2:39866/72‑1 (MQ=255)
                                                              ggCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCAAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGtgcg                                                                                                                                                    >  2:374806/1‑135 (MQ=255)
                                                               gcgcGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAAcccc                                                                                                                                                                                         >  2:469546/1‑97 (MQ=255)
                                                                  cgATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCAAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGTGCGATGc                                                                                                                                                <  1:374819/135‑1 (MQ=255)
                                                                   gATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTAcc                                                                                                                                                                                     <  1:469559/97‑1 (MQ=255)
                                                                    aTAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGNTTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATATCCGCAGACAATATAAAATCGGTACGAATCGCCCCATTAATTTTGTCCTTTTCAAACTACAGCGGacactg     >  1:423043/1‑267 (MQ=255)
                                                                    aTAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACgg                                                                                                                 >  1:256881/1‑164 (MQ=255)
                                                                    aTAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTCCGGCCAAAATCCCCATCGTGATGGAGACGATTCCCGCAAACAATATAAAATCGGTACCANTCCCCCCTTTTATTTGTTCNTTTTTTAACTTNaacctttcccttcc  >  2:250019/1‑262 (MQ=255)
                                                                        tCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCt                                                                                                             <  2:256873/164‑1 (MQ=255)
                                                                            tAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCgccagcc                                                                                                                                                                                                                           <  2:331033/50‑1 (MQ=255)
                                                                                          aTGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAAt                                                                                                                                                                                                                        >  2:111660/1‑39 (MQ=255)
                                                                                          aTGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAAt                                                                                                                                                                                                                        >  1:168270/1‑39 (MQ=255)
                                                                                          aTGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTa                                                                                                                                                                                                    >  1:21568/1‑59 (MQ=255)
                                                                                          aTGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCt                                                                                                                                                                                        >  1:306934/1‑71 (MQ=255)
                                                                                              cTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATAccc                                                                                                                                                                                                                    <  2:168264/39‑1 (MQ=255)
                                                                                              cTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATAccc                                                                                                                                                                                                                    <  1:111666/39‑1 (MQ=255)
                                                                                              cTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGa                                                                                                                                                                                                <  2:21567/59‑1 (MQ=255)
                                                                                              cTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAATCCCTACCa                                                                                                                                                                                    <  2:306925/71‑1 (MQ=255)
                                                                                                gCATCAGCGCGCTGGATTTTTCCGCCAGGCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAagcacc                                                                                                                                     >  1:372289/1‑116 (MQ=255)
                                                                                                       cgcgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGcg                                                                                  >  1:24038/1‑160 (MQ=255)
                                                                                                           cTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGcgacga                                                                              <  2:24037/160‑1 (MQ=255)
                                                                                                            tGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGaa                                                >  1:453973/1‑189 (MQ=255)
                                                                                                                tttttCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGc                                            <  2:453960/189‑1 (MQ=255)
                                                                                                                             |                                                                                                                                                                                                                         
TCCCACCGCCGACGAGGAACATCGCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTCAAACTTCAGCGGATCCTGCGC  >  CP009273/2022109‑2022451

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