Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 780,262 A→G K222R (AAA→AGA)  cpoB → cell division coordinator CpoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913780,2620AG95.9% 156.8 / ‑3.3 49K222R (AAA→AGA) cpoBcell division coordinator CpoB
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (23/24);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTACCTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAAAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATGTAATGC  >  NC_000913/780126‑780392
                                                                                                                                        |                                                                                                                                  
aCTTACCTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGt                                                                                                                                 >  2:137849/1‑140 (MQ=255)
   tACCTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCgg                                                                                                                              <  1:854755/140‑1 (MQ=255)
     ccTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCg                                                                                                                            >  2:640371/1‑140 (MQ=255)
     ccTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTCCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCg                                                                                                                            >  1:708597/1‑140 (MQ=255)
      cTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGt                                                                                                                           >  1:181736/1‑140 (MQ=255)
        gCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGtca                                                                                                                         <  2:181736/140‑1 (MQ=255)
        gCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGtca                                                                                                                         <  2:49650/140‑1 (MQ=255)
        gCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGAAGCGATGTTTAGAGTCGGCGtca                                                                                                                         <  1:362707/140‑1 (MQ=255)
         ccAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGtcat                                                                                                                        >  2:370009/1‑140 (MQ=255)
                   attattGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAgg                                                                                                                <  2:812251/138‑1 (MQ=255)
                   attattGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAgg                                                                                                                >  1:812251/1‑138 (MQ=255)
                               gTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACAcc                                                                                                  >  2:749278/1‑140 (MQ=255)
                               gTCAGTTAAACTACAACAAGGGGAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACAcc                                                                                                  >  2:569977/1‑140 (MQ=255)
                                            caacaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGaa                                                                                       >  1:690697/1‑138 (MQ=255)
                                            caacaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGaa                                                                                       <  2:690697/138‑1 (MQ=255)
                                            caacaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAg                                                                                     <  1:749278/140‑1 (MQ=255)
                                            caacaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAg                                                                                     >  1:382474/1‑140 (MQ=255)
                                              acaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGcc                                                                                   >  1:417457/1‑140 (MQ=255)
                                              acaaGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGcc                                                                                   >  2:248167/1‑140 (MQ=255)
                                                   ggTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGtata                                                                              <  1:248167/140‑1 (MQ=255)
                                                    gTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATAc                                                                             >  2:800829/1‑140 (MQ=255)
                                                            gatgatGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGt                                                                     <  1:800829/140‑1 (MQ=255)
                                                                atgCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATc                                                                 <  1:370009/140‑1 (MQ=255)
                                                                      cgTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTaaa                                                           <  1:186994/140‑1 (MQ=255)
                                                                      cgTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTaaa                                                           <  2:821992/140‑1 (MQ=255)
                                                                        tACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATa                                                         <  2:14710/140‑1 (MQ=255)
                                                                         tcTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATAc                                                        <  2:496254/139‑1 (MQ=255)
                                                                           tATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATAccc                                                      >  2:178293/1‑140 (MQ=255)
                                                                                     cGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGAt                                            >  1:119628/1‑140 (MQ=255)
                                                                                     cGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGAt                                            >  1:93487/1‑140 (MQ=255)
                                                                                       gTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATgg                                          >  2:637129/1‑140 (MQ=255)
                                                                                              aaaaCTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTaaa                                   <  2:93487/140‑1 (MQ=255)
                                                                                              aaaaCTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTaaa                                   <  1:637129/140‑1 (MQ=255)
                                                                                               aaaCTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAAc                                  <  1:137849/140‑1 (MQ=255)
                                                                                                  cTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAgg                               >  1:785052/1‑140 (MQ=255)
                                                                                                          aGTCACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACaaaaa                       >  2:849584/1‑140 (MQ=255)
                                                                                                             cACCAAAGGCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGt                    >  2:432570/1‑140 (MQ=255)
                                                                                                                    ggCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcg             <  1:413487/140‑1 (MQ=255)
                                                                                                                    ggCTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcg             <  1:178293/140‑1 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcgcg           <  1:849584/140‑1 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcgcg           <  2:708428/140‑1 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcgcg           <  1:227170/140‑1 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcgcg           >  2:297978/1‑140 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcg             >  1:207739/1‑138 (MQ=255)
                                                                                                                      cTGCAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAAcg             <  2:207739/138‑1 (MQ=255)
                                                                                                                         cAGATGCGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATg        <  1:297978/140‑1 (MQ=255)
                                                                                                                               cGATGTTTAGAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATGTAATGc  >  2:323147/1‑140 (MQ=255)
                                                                                                                               cGATGTTTAAAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATg        >  1:395058/1‑134 (MQ=255)
                                                                                                                               cGATGTTTAAAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATg        <  2:395058/134‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
ACTTACCTGCCAAACGCCAATTATTGGCTGGGTCAGTTAAACTACAACAAGGGTAAAAAAGATGATGCGGCGTACTATTTTGCTTCGGTAGTGAAAAACTATCCGAAGTCACCAAAGGCTGCAGATGCGATGTTTAAAGTCGGCGTCATCATGCAGGACAAAGGTGACACCGCAAAAGCGAAAGCCGTATACCAGCAGGTTATCAGTAAATACCCTGGTACTGATGGCGCTAAACAGGCACAAAAACGTCTGAACGCGATGTAATGC  >  NC_000913/780126‑780392

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