Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,710,898 (T)6→7 coding (29/1812 nt) oatA ← possible acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,710,8921.T95.7% 63.5 / ‑3.4 23coding (35/1812 nt)oatApossible acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (16/6);  total (17/6)
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

AACCTGAAATAACGAAGAATGTATCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCACTGATTTTCATTAATTGATTTAACCCTTA  >  CP000730/2710755‑2711027
                                                                                                                                          |                                                                                                                                       
aaCCTGAAATAACGAAGAATGTATCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACAttt                                                                                                                                      <  2:319224/141‑1 (MQ=255)
              aagaaTGTATCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCGGGTAAGTATCGCGGGGAATACATTTTTTCTAAACGTTTa                                                                                                                        >  1:68309/1‑141 (MQ=255)
               agaaTGTATCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTGAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTa                                                                                                                        <  1:76137/141‑1 (MQ=255)
                       tCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCttt                                                                                                                <  1:113634/141‑1 (MQ=255)
                                         ccACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGAAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATaaa                                                                                              >  2:376592/1‑141 (MQ=255)
                                               aaCCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGcccc                                                                                        >  2:59604/1‑141 (MQ=255)
                                                          aTTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTtc                                                                             >  2:21079/1‑141 (MQ=255)
                                                            tCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTtctc                                                                           <  1:183190/141‑1 (MQ=255)
                                                                gTGATAAATAATGATTCCAATAACTGCGAATGCCCTCACTCCATCTACTCCAGGTAAGTATCGCGGGGAATACATTTTTTTCGAAACGTTTACAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTAtt                                                                       >  2:262957/1‑141 (MQ=255)
                                                                    taaataaTGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTAt                                                                                                           <  2:398177/101‑1 (MQ=255)
                                                                    taaataaTGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTAt                                                                                                           >  1:398177/1‑101 (MQ=255)
                                                                       ataatGATTCCTATAACTGCGAATGCCCTCAATCCATCTAGTCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGc                                                                >  1:419042/1‑141 (MQ=255)
                                                                             aTTCCTATAACTGCGAATGCCCTCAATTCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATca                                                          >  2:275419/1‑141 (MQ=255)
                                                                                    taACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATAtt                                                   >  1:69272/1‑141 (MQ=255)
                                                                                          cGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAATGTTTAAAGTCATTTGTATCCATGTTAATAAACGCCCAATTTATTTTTCTCTATTTTATAGAATTTCACAATATTTTTGaa                                             >  2:103013/1‑141 (MQ=255)
                                                                                                        ccATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCa                               <  2:381811/141‑1 (MQ=255)
                                                                                                          aTCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATGTATTTTTCTCTATTTTGTAGCATATCACAATATCTTTGAAAATAAAATATTGCACt                             <  2:429569/141‑1 (MQ=255)
                                                                                                           tCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAACTAAAATATTGCACTg                            >  2:418007/1‑141 (MQ=255)
                                                                                                            cTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGAATATCACAATATTTTTGAAAATAAAATATTGctttga                           >  1:200462/1‑136 (MQ=255)
                                                                                                             tAATCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCACTGAt                          >  1:381839/1‑141 (MQ=255)
                                                                                                              aaTCCAGGTAAGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTTTTTTTTTAATATTGCACTGAtt                         >  2:149851/1‑141 (MQ=255)
                                                                                                                       aaGTATCGCGGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCACTGATTTTCATTAAt                >  1:121015/1‑141 (MQ=255)
                                                                                                                              gcgGGGAATACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAAAGCCCCATTTATTATTTACTATTTGGTAGCATAGCACAATATTTTTGAATATAAAATATTGCACTGATTTTCATTAATTGATTTa         >  1:237122/1‑141 (MQ=255)
                                                                                                                                     aTACATTTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCACTGATTTTCATTAATTGATTTAACCCtta  >  1:171290/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                       
AACCTGAAATAACGAAGAATGTATCTACTCCTAAAAAGCCCCCACTTAACCATTGTGCATTCAAGTGATAAATAATGATTCCTATAACTGCGAATGCCCTCAATCCATCTAATCCAGGTAAGTATCGCGGGGAATACATTTTTTCTAAACGTTTAAAGTCTTTTGTATCCATGTTAATAAACGCCCCATTTATTTTTCTCTATTTTGTAGCATATCACAATATTTTTGAAAATAAAATATTGCACTGATTTTCATTAATTGATTTAACCCTTA  >  CP000730/2710755‑2711027

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