Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,506,454 C→T 28.5% intergenic (‑269/‑67) mntH ← / → nupC manganese/divalent cation transporter/nucleoside (except guanosine) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,506,4540CT28.5% 50.0 / 19.4 35intergenic (‑269/‑67)mntH/nupCmanganese/divalent cation transporter/nucleoside (except guanosine) transporter
Reads supporting (aligned to +/- strand):  ref base C (9/16);  new base T (3/7);  total (12/23)
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

GAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGC  >  CP009273/2506320‑2506595
                                                                                                                                      |                                                                                                                                             
gAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCa                                                                                                                                 >  1:425871/1‑149 (MQ=255)
         ggttatgTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATa                                                                                                                        >  1:118402/1‑149 (MQ=255)
          gttatgTTTAATGTTAGCACCTTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACAtcg                                                                                                                       <  2:135464/149‑3 (MQ=255)
            tatgTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGaa                                                                                                                     <  2:25840/149‑1 (MQ=255)
            tatgTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGaa                                                                                                                     <  1:115983/149‑1 (MQ=255)
            tatgTTTAATGTTAGCACATCTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGaa                                                                                                                     <  2:365494/149‑1 (MQ=255)
             atgCTTAATGTAAGCACATTTACATAACTTTCAGCTTCCATACACAATCTAGCAGAAACGTATGACAGATCACTATTTTGGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAAt                                                                                                                    <  1:130255/149‑1 (MQ=255)
                tttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAAtgtg                                                                                                                 >  1:444687/1‑149 (MQ=255)
                         aGCACATTTACATCACTTTCAGCTTCCATACACAACACAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGCCCCATGCCGCCATTATACTCTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGaact                                                                                                        <  2:157076/149‑3 (MQ=255)
                           cacaTTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAAttt                                                                                                      <  2:142079/149‑1 (MQ=255)
                            acaTTTACATAACTTTCAGCTTCGATACACAACAAAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTa                                                                                                     >  1:257524/1‑149 (MQ=255)
                               tttACATAACTTTCAGCTTCCACACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCATATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAAcc                                                                                                  <  2:125259/149‑1 (MQ=255)
                                     tAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCTTGTTCCCTCGAACGTCCACGGAACTTTAACCTGcccc                                                                                            >  2:229835/1‑147 (MQ=255)
                                            cAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTg                                                                                     <  1:422856/149‑1 (MQ=255)
                                            cAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTg                                                                                     <  2:286032/149‑1 (MQ=255)
                                              gCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGa                                                                                   <  2:260530/149‑1 (MQ=255)
                                                   cATACAAAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTGGCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCaaa                                                                              <  1:225072/149‑1 (MQ=255)
                                                     tACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCGTTGTAGTGGGTGTCAGATATCGTTGTCCTTAACCATGTTGCATAGGATGTGCACGGAAATTTAACCTGCCTCATGTTGGGAGCAAata                                                                            >  2:515/1‑149 (MQ=255)
                                                        acaacaTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTCGAGCAAATATgg                                                                         >  1:244833/1‑149 (MQ=255)
                                                          aacaTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGAc                                                                       <  1:324124/149‑1 (MQ=255)
                                                                   gAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCtt                                                              <  2:118402/149‑1 (MQ=255)
                                                                            gACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTa                                                     >  1:436739/1‑149 (MQ=255)
                                                                              cAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACt                                                   <  1:29357/149‑1 (MQ=255)
                                                                               aGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATTTTCACGTAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTg                                                  <  2:425871/149‑1 (MQ=255)
                                                                               aGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTGACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTg                                                  <  2:261869/149‑1 (MQ=255)
                                                                                      tATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTg                                           >  1:291615/1‑149 (MQ=255)
                                                                                                gCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGc                                 >  1:133465/1‑149 (MQ=255)
                                                                                                         atgACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGc                        <  1:393639/147‑1 (MQ=255)
                                                                                                          ggACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCa                       <  2:369319/149‑1 (MQ=255)
                                                                                                              gTCATCAGAGTGTGTGTCAGATCTTGTCTTCCTTAACCATTTTACATAGCATGTTCACGGAAATTTACCCCGCCTCATTTTTGGAGCAACTATGGACCGCGTTCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCCctgc                   <  2:178199/149‑1 (MQ=255)
                                                                                                              atCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCActgc                   <  2:380565/148‑1 (MQ=255)
                                                                                                                 aTTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGACTCATCTTTGGAGCAACTCTGGACCGCGTCCTTCATTTTGTAATGGCACTTGCCGTTGTTGCGATTCTCGAACTGCggg                >  1:96682/1‑149 (MQ=255)
                                                                                                                       gtgtgtgtCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAagca          <  2:18640/149‑1 (MQ=255)
                                                                                                                          tgtgtCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCg       <  2:436739/149‑1 (MQ=255)
                                                                                                                               cAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTCCGATTCTCGCACTGCTGGTAAGCAGCGACCGc  >  2:59458/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                             
GAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGC  >  CP009273/2506320‑2506595

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: