Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 547,929 A→G 13.4% K191K (AAA→AAG ftsH → cell division protein FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730547,9290AG13.4% 99.4 / 6.7 44K191K (AAA→AAGftsHcell division protein FtsH
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (3/3);  total (23/21)
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

TAAATCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGATGT  >  CP000730/547800‑548068
                                                                                                                                 |                                                                                                                                           
aaaaTCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTCTCTCTGATGTAGCGGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATc                                                                                                                                  <  2:106484/140‑1 (MQ=255)
        aaGCAAAAATGTACGATAATAATAGACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTc                                                                                                                          >  2:179122/1‑141 (MQ=255)
        aaGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTCATTTCTTGAAAGATAATAAAAAATTCGAAGAGATGGGATCTAGGATTc                                                                                                                          >  1:49017/1‑141 (MQ=255)
          gCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCt                                                                                                                        <  1:134403/141‑1 (MQ=255)
           cAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTa                                                                                                                       >  1:23293/1‑141 (MQ=255)
              aaaTGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGCTGATTTCATGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAg                                                                                                                    >  1:282240/1‑141 (MQ=255)
                aTGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGgt                                                                                                                  >  1:317149/1‑141 (MQ=255)
                     cGATAATAATAAACGGCGGGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCtt                                                                                                             <  1:100245/141‑1 (MQ=255)
                       ataataataAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTat                                                                                                           >  1:326344/1‑140 (MQ=255)
                           taataaACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACttgt                                                                                                       <  2:55915/141‑1 (MQ=255)
                                 aCGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAACAATTCACAGAAATGGGATCTAGGTTTCCTAAAGGTGCCCTACttcgtggccg                                                                                                 <  2:378899/141‑9 (MQ=255)
                                 aCGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGGCTTACTTGTTGGAcc                                                                                                 >  1:39868/1‑141 (MQ=255)
                                        gtTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTAGGACCTCCAggg                                                                                          >  2:55879/1‑140 (MQ=255)
                                           cGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACt                                                                                       >  1:316277/1‑141 (MQ=255)
                                           cGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACt                                                                                       >  1:143045/1‑141 (MQ=255)
                                           cGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTGGATTTCTTGAAAGATAATAAAAAATTCAAGGAAAAGGGATTTAGGATTCCTAAAGGAGTCTTACTCGTTggcgcgccaaggaat                                                                                       >  1:185787/1‑128 (MQ=255)
                                                 tctGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCAAGGTACTGGTaaa                                                                                 >  1:144692/1‑141 (MQ=255)
                                                   tGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAAc                                                                               <  2:62359/141‑1 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:158912/1‑141 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:169961/1‑141 (MQ=255)
                                                          gCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTATTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACtt                                                                        >  1:170053/1‑141 (MQ=255)
                                                             ggggCAGATGCAGAAAAACAAGAATGAATTGAAATTGTTGATTGCGTGAAAGATAACAAAAAATGCAAAGAGATGGGATCTAGGATTCCTAAAGGTGGCTTACTTGTTGGACCTCCAGGTACTGGTAAAAAGTTACTTGCt                                                                     >  2:90693/1‑141 (MQ=255)
                                                               ggCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAGTGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTag                                                                   <  1:10587/141‑1 (MQ=255)
                                                               ggCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTag                                                                   <  1:10585/141‑1 (MQ=255)
                                                                  aGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGc                                                                <  1:301898/141‑1 (MQ=255)
                                                                           aaaaCAAGAATTAATTGAAATTGTTACTGGTTTTGAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGATGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATGACTTGCTAGAGCGGTTGGAgg                                                       >  2:175489/1‑141 (MQ=255)
                                                                             aaCAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTg                                                     >  2:233847/1‑141 (MQ=255)
                                                                                           gAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCa                                       <  1:253596/141‑1 (MQ=255)
                                                                                           gAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATGACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCa                                       <  2:157486/141‑1 (MQ=255)
                                                                                                            gAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTgg                      >  1:383277/1‑141 (MQ=255)
                                                                                                              aaGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGtt                    >  2:369853/1‑141 (MQ=255)
                                                                                                                   aataaAACATTTAAGGAAAGGGGATCTAGGAGTCCTAAAGTTGTCCTACTAGTGGCTCCTCCAGGTACTTTTAACATTTTGCTTGCTAGAGCGGTTGCAGGTGAAGCTCGCGCACCATTCTTCTCTATTAGTGGTTCATAc               <  2:185787/141‑1 (MQ=255)
                                                                                                                   aataaAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGAc               <  2:143045/141‑1 (MQ=255)
                                                                                                                   aataaAAAATTCAAGGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGAc               <  2:316277/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACt              <  2:169961/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACt              <  2:170053/141‑1 (MQ=255)
                                                                                                                    ataaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCACTCCATTCTTCTCTATTAGTGGTTCAGACt              <  2:158912/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCGGGTGACGCTGGCGCACCATTCTTCTCTATTGGTGGTTCAggagt             >  2:111892/1‑137 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:139266/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:139241/141‑1 (MQ=255)
                                                                                                                     taaAAAATTCAAAGAAATGGGATCGAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACtt             <  2:49017/141‑1 (MQ=255)
                                                                                                                      aaaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTGTTCTCTATTAGTGGTTCAGACttt            >  2:35277/1‑141 (MQ=255)
                                                                                                                       aaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTg           <  1:20681/141‑1 (MQ=255)
                                                                                                                       aaaaaTTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAAGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACttgg           >  2:96800/1‑139 (MQ=255)
                                                                                                                                aaGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGAtgt  >  1:67746/1‑141 (MQ=255)
                                                                                                                                 |                                                                                                                                           
TAAATCTAAAGCAAAAATGTACGATAATAATAAACGTCGTGTTCGTTTCTCTGATGTAGCAGGGGCAGATGAAGAAAAACAAGAATTAATTGAAATTGTTGATTTCTTGAAAGATAATAAAAAATTCAAAGAAATGGGATCTAGGATTCCTAAAGGTGTCTTACTTGTTGGACCTCCAGGTACTGGTAAAACATTACTTGCTAGAGCGGTTGCAGGTGAAGCTGGCGCACCATTCTTCTCTATTAGTGGTTCAGACTTTGTAGAGATGT  >  CP000730/547800‑548068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: