Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,353 A→G 72.2% intergenic (‑86/+135) hisJ ← / ← argT histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,420,3530AG72.2% 54.1 / 21.3 40intergenic (‑86/+135)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (2/9);  new base G (20/8);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  CP009273/2420220‑2420498
                                                                                                                                     |                                                                                                                                                 
gcTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGtt                                                                                                                                    >  2:375657/1‑149 (MQ=255)
  tGGAGAAGGCCAGAACCAGAGGGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGccc                                                                                                                                        <  1:63149/143‑1 (MQ=255)
  gggggaAGGCCACAACCAGAGCGAGCGATCGCACCAGTTTTTTCGTAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTGACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAg                                                                                                                                  <  2:149447/145‑1 (MQ=255)
    gagaAGGCCAGAACCAGCAGGAGCGATAGCACCAGTTTTTTTATGTTGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCACGGCAAGATATCCCGCCCCGTTAGAc                                                                                                                                <  2:216994/149‑1 (MQ=255)
          gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTTTTTTTTTTTTCAACTGACAGTCGTTTTATAGGGAATATTTACAGGCGTGTGACACTTATCGTGCCTTCAAATTGCCCATCAAGGCAAGAGATCCCGCCCCGtaggac                                                                                                                                >  1:226152/1‑138 (MQ=255)
          gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCTTTATTTTTTCAACTGACAGTCCTTTTATAGGGAATCTTTACAGGTTTGTTACACTTATCGTGCCATACAATTGACGAACACGGCAAGATATCCCGCCACCTTAGAc                                                                                                                                >  1:112997/1‑143 (MQ=255)
          gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATACCGCCCCGTTAGACGtatta                                                                                                                          >  2:262451/1‑149 (MQ=255)
          gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTAGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGAc                                                                                                                                >  1:288987/1‑143 (MQ=255)
          gCCAGAACCAGAGAGAGCGATAGCACCAGGTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCCTAAAATTGACCATCAAGGCAAGATATCCCGcccccgtaggc                                                                                                                                >  1:143276/1‑136 (MQ=255)
              gAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGtat                                                                                                                            <  1:368273/143‑1 (MQ=255)
                       agaggcattgCACCAGTTTTTTCATAAGGTCCTCACCTGACAGCATTTTTATAGGGAATATTTACAGGTTTGTGATACTTATCGTGCCATAAAATTGACCATCAAGTCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCa                                                                                                                   <  1:36117/134‑1 (MQ=255)
                       agagGGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTCTAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATCTCCCGCCCCGTTAGACGTATTATTTTTCa                                                                                                                   <  1:11226/143‑1 (MQ=255)
                       agagCGATAGCACCATTTTTTTCCTAAGGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTTACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCa                                                                                                                   <  1:32276/143‑1 (MQ=255)
                       agagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTaaaa                                                                                                             >  2:390315/1‑149 (MQ=255)
                                 cACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGCCCATCAAGGTAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGc                                                                                                         <  1:355664/143‑1 (MQ=255)
                                  aCCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCt                                                                                                        >  1:118444/1‑143 (MQ=255)
                                              atgaTGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATg                                                                                      <  2:175003/146‑1 (MQ=255)
                                              aTAATGTCCTCAACTGTCAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATg                                                                                      <  2:247710/149‑1 (MQ=255)
                                                       taaaCTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGcac                                                                                   <  1:180232/141‑1 (MQ=255)
                                                            tGACAGTCTTTTTATAGGGAATATTTCCACGAAGACTATAGTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTCGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACgc                                                                              >  1:175606/1‑142 (MQ=255)
                                                                     ttttATAGGGAATATTTACAGGTTTGTTAAACAACAAATGTCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaac                                                                     >  1:273692/1‑143 (MQ=255)
                                                                     ttttATAGGGAATATTTACAGGGTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGAAGTATTATTTTTCATTAAAAATGCTGAAGGATGCATTTCAATGCACGACGATAGTGcaac                                                                     >  1:96632/1‑143 (MQ=255)
                                                                      tttATAGGGAATTTTGCCAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTg                                                              <  2:354204/149‑1 (MQ=255)
                                                                               aaTATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCg                                                           >  1:198876/1‑143 (MQ=255)
                                                                                  atTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatg                                                  >  2:355720/1‑149 (MQ=255)
                                                                                  atTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGAGCACCatg                                                  >  2:355726/1‑149 (MQ=255)
                                                                                   tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACTACGATAGTGCAACAAACCTGGCGTGCa                                                       >  1:373734/1‑143 (MQ=255)
                                                                                          gTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCGTGATGATcct                                          >  2:88851/1‑147 (MQ=255)
                                                                                              gTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTCAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCCACAAACCTGACGCGCACCATGATGATCATTCCa                                      >  2:406848/1‑149 (MQ=255)
                                                                                                 acacTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCa                                   <  2:49826/149‑1 (MQ=255)
                                                                                                  cacTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAg                                  <  2:340736/149‑1 (MQ=255)
                                                                                                  aactttctgTGCCATAAAATTGACCGTAAAGGTAGTATATCCGGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAg                                  <  2:205529/141‑1 (MQ=255)
                                                                                                      tacagTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATc                                    <  1:230714/139‑1 (MQ=255)
                                                                                                              cATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAATCCATCAGGTCCAg                            >  1:391886/1‑143 (MQ=255)
                                                                                                                     ttGACCATCAAGGCAAAATATCCCGCCCCGTTCGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACACCCCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCa                     <  1:133443/143‑1 (MQ=255)
                                                                                                                          cATCAAGGCAAGATATCCCGCCCCGTTAGACGTGTTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGAc                <  1:194716/143‑1 (MQ=255)
                                                                                                                          cATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCa          >  2:124143/1‑149 (MQ=255)
                                                                                                                            tCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCCCCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCGATCAGGTACAGCTTCCCAGCGACGTATCAGt        >  2:66651/1‑149 (MQ=255)
                                                                                                                                 gCAAGATATCCCGCCCCGTTAGACGTATTATCTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGCTACAGCTTCCCCGCGACGTATCAg         >  1:360384/1‑143 (MQ=255)
                                                                                                                                  cAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTAAccat  >  2:156507/1‑147 (MQ=255)
                                                                                                                                     |                                                                                                                                                 
GCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  CP009273/2420220‑2420498

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