Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,933,885 G→T L30I (CTT→ATT)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,8850GT90.0% 89.5 / ‑2.4 42L30I (CTT→ATT) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (2/0);  major base T (21/15);  minor base A (1/3);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGT  >  CP009273/3933751‑3934007
                                                                                                                                      |                                                                                                                          
aCCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTGTTAataata                                                                                                                        >  1:2251919/1‑139 (MQ=255)
aCCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAataata                                                                                                                        <  1:2107386/139‑1 (MQ=255)
        gATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATc                                                                                                                <  1:2878104/139‑1 (MQ=255)
         aTCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCa                                                                                                               >  1:345268/1‑139 (MQ=255)
          tCGGCTAGCCATCCGCTTACCGGAATAAGCATCGACACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGAAAGGCGGGATTAATAATAGTGGCATCaa                                                                                                              >  1:1058520/1‑139 (MQ=255)
              cTAGCCCTCCGCTTACCGGAATAAGCAGCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTg                                                                                                           <  1:3230269/138‑1 (MQ=255)
               tAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGc                                                                                                          <  1:921201/138‑1 (MQ=255)
                   cATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCa                                                                                                     <  1:766311/139‑1 (MQ=255)
                           tACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTg                                                                                                                     >  1:657762/1‑115 (MQ=255)
                           tACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAAAATAGTGGGATCAAGGGCCTGCATGaagaag                                                                                             >  1:4001621/1‑139 (MQ=255)
                            aCCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGaagaag                                                                                             >  1:997636/1‑138 (MQ=255)
                                 aaTAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATc                                                                                       >  1:1966399/1‑139 (MQ=255)
                                 aaTAAGCATCGCCACCGACAGCGGATAACTGATGATGTCTGATTGCATCGCGGGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCGAGAAGAAGGCCaac                                                                                       >  1:3013953/1‑137 (MQ=255)
                                        aTCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATtaataa                                                                                                                          >  1:583617/1‑97 (MQ=255)
                                           gCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAaaaaaa                                                                                                                        >  1:2573891/1‑94 (MQ=255)
                                           gCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATcc                                                                             >  1:417158/1‑139 (MQ=255)
                                                cGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGc                                                                        >  1:1457387/1‑139 (MQ=255)
                                                  tCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCaa                                                                      >  1:3158032/1‑139 (MQ=255)
                                                                atgatgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAAACCGCCATACTgc                                                        >  1:2885728/1‑139 (MQ=255)
                                                                atgatgGCTGATAGAATCGCGAGAGTAGAACGATTCAGGCTATGGGCGATTGCGGGTAAGGCGGTATTAATAATAGGGGCGTCGAGTGCGTGCCTGAAAAAGGCCATCGCCGCGATCCACAGGAAACCCGCCCtgctgc                                                        >  1:4076456/1‑139 (MQ=255)
                                                                 tgatgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATGGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCGCGGCAAACCCGCCATACTgcg                                                       <  1:1503087/139‑1 (MQ=255)
                                                                 tgaggGCTGATTGCATCGGGAGAGGAGAACGATTAAGGCTATGAGCGTTTGGGGGTAGGGGGGTATTAAAAATGGGGGCATCAAGTCCCTGCAGGAAGAAGGCCATCGCCGCGATCCCCGGAAACCCCGCCATACTgcg                                                       <  1:3169477/139‑1 (MQ=255)
                                                                          aTTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCttttt                                              >  1:1665301/1‑139 (MQ=255)
                                                                                tCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAataata                                                                                                                        >  1:1720802/1‑59 (MQ=255)
                                                                                 cgcgGGAGGAAAACGATTAGGGAATGGCGGGATGGGGGGAAAGGGGGTATAAATAATAGTGGCAACAAGTGCCTGCATGAAGAGGTCCACCGCCGCGACCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTc                                       <  1:3626644/139‑1 (MQ=255)
                                                                                   cgAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAataata                                                                                                                        >  1:7304/1‑56 (MQ=255)
                                                                                    gagaGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAg                                                                                                                       >  1:2833597/1‑56 (MQ=255)
                                                                                     agagGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCAAACTGCGCTTCTTTTTATCGCCCATTc                                   >  1:3751219/1‑139 (MQ=255)
                                                                                            aaCGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATc                                                                                                                <  1:754663/55‑1 (MQ=255)
                                                                                              cGATTAAGGCTAGGAGCTATTGCGGGTAAGGGGGTATTAATAATAGGGGCCTCAAGTGCCTGCATGAAGAAGCCCATCGCCGCGACCCACGGCAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTg                          <  1:1831402/139‑1 (MQ=255)
                                                                                                aTTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGc                                                                        <  1:160877/91‑1 (MQ=255)
                                                                                                 ttAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGGCCTGTTa                       >  1:3411712/1‑139 (MQ=255)
                                                                                                  tAAGGCTATGGGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGCCCATCGCCGCGATCCCCGGCAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTAt                      <  1:986172/139‑1 (MQ=255)
                                                                                                   aaGGCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATc                     >  1:895610/1‑139 (MQ=255)
                                                                                                      gCTATGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGc                                                                                                      >  1:2926843/1‑55 (MQ=255)
                                                                                                       cTAGGAGCGATTGCGGGAAAGGGGTTAATAAAAAGGGGGGAACCAAGTGCCTGCATGAAGAGGCCCTTCGCCGCGACCCACGGAAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGt                 <  1:3628361/139‑1 (MQ=255)
                                                                                                        tATGGGCGTTGGGGGGAAGGGGGGTATTAAAAATAGTGCCATCAAGTCCCGGCATGAAGAGGCCTATCGCCGCGATCCACGGAAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGt                 <  1:2518073/138‑1 (MQ=255)
                                                                                                         aTGAGCGATTGCGGGTAAGGCGGTTTTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTa               >  1:494196/1‑139 (MQ=255)
                                                                                                          tGAGCGATTGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCa                                                              >  1:2333966/1‑91 (MQ=255)
                                                                                                           gAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATc             >  1:2466343/1‑139 (MQ=255)
                                                                                                                 ttGCGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCACTGTCCTGTTATCGGGTTATCACTTa        >  1:4236783/1‑138 (MQ=255)
                                                                                                                    ggggTAGGGGGTTATAAATAAGAGGGGAATAAGGTGCCTGCATGAAGAGGGCCACCGCCGCGACCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCa     <  1:3417259/137‑1 (MQ=255)
                                                                                                                    cGGGTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCa                                  >  1:2145569/1‑109 (MQ=255)
                                                                                                                      ggTAAGGCGGTATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGATTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGt  >  1:287999/1‑139 (MQ=255)
                                                                                                                                tATTAATAATGGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGAAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTc                             <  1:1279436/102‑1 (MQ=255)
                                                                                                                                tATTAATAATATTGGCATAAAGTGCCTCCATGAAGAAGGCCATCGCCGCGATCCACGGCAACCCCGCCATACTGCGCTTCTTTTTATCGCTc                                       <  1:2169331/92‑1 (MQ=255)
                                                                                                                                tATTAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATcgccg                                                                                   <  1:3341424/48‑1 (MQ=255)
                                                                                                                                  ttAATAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGAAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGt  <  1:790201/127‑1 (MQ=255)
                                                                                                                                      |                                                                                                                          
ACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGT  >  CP009273/3933751‑3934007

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