Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,752 C→A 12.1% G165C (GGC→TGC)  modF ← fused molybdate transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273788,7520CA12.1% 71.5 / 2.4 33G165C (GGC→TGC) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base C (18/11);  new base A (2/2);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGC  >  CP009273/788613‑788894
                                                                                                                                           |                                                                                                                                              
cAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGa                                                                                                                                           >  2:273357/1‑145 (MQ=255)
 aaTCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGaa                                                                                                                                          <  2:361492/145‑1 (MQ=255)
    ccGCCAGCAGGCCAGCAAACTGGACAACCTCCGGGATCTGATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACGGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACgg                                                                                                                                       <  1:98159/145‑1 (MQ=255)
       ccagcacgccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCtc                                                                                                                                    <  2:322148/145‑1 (MQ=255)
       ccagcacgccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCAATCGAACGGCtc                                                                                                                                    >  1:16912/1‑145 (MQ=255)
       ccagcacgccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCAATCGAACGGCtc                                                                                                                                    >  1:74603/1‑145 (MQ=255)
              gccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTCCCAGCGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGa                                                                                                                             >  1:278652/1‑145 (MQ=255)
               ccagcaAACTGGACAAACTCCGGGATCTCATCGGAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGaa                                                                                                                            <  2:207546/145‑1 (MQ=255)
                cagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACAAGGCGAGTCGCTCAGCCAGCTGCTGACGTAAGGCAACATCCAGGCCATCGAACGGCTCATCAAAAAt                                                                                                                           >  2:137641/1‑145 (MQ=255)
                 agccAACTGGACCAACGCCGGCATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGGAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATc                                                                                                                          <  1:330872/145‑1 (MQ=255)
                    aaaCTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCAGCAAGAATcaac                                                                                                                       >  2:278364/1‑145 (MQ=255)
                     aaCTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATcaaca                                                                                                                      >  2:24741/1‑145 (MQ=255)
                                           gaTCGAAGCGATTGAGCACCAGTACCAGGGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGACAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGACTCCGACATCAgcga                                                                                                >  1:390944/2‑144 (MQ=255)
                                            aTCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATGAAGAATCAACAAGTCAGGCTCCGACATCAGCGcc                                                                                               >  1:134564/1‑145 (MQ=255)
                                                              cAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCAATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                             <  2:74603/145‑1 (MQ=255)
                                                              cAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCAATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                             <  2:16912/145‑1 (MQ=255)
                                                                ggaCCAGGATAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTcg                                                                           <  1:48685/143‑1 (MQ=255)
                                                                 tACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTcgc                                                                          >  1:51777/1‑145 (MQ=255)
                                                                            caCCGGACGGAGGTAACGAGCCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGt                                                               <  1:79333/144‑1 (MQ=255)
                                                                               cGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGa                                                            >  1:333654/1‑145 (MQ=255)
                                                                                      aTGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGCACGGCTCATCAGGAATCAACAAGTCAGGCTCCGACATCAGCGCCGGACACAGACGGGTTTTGCGCGTCTCGCCAGTGGAAAGGTAt                                                     >  1:140928/1‑145 (MQ=255)
                                                                                           aCGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTaaa                                                >  1:197663/1‑145 (MQ=255)
                                                                                             gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGc                                              >  2:403436/1‑145 (MQ=255)
                                                                                             gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACACCAGCGCCTGACACAGCAGGGCTTTTCGCGTCCCGCCAGCGGAAAGGGATGTAAAGc                                              >  1:265143/1‑145 (MQ=255)
                                                                                                gCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCCAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCctc                                           >  2:39725/1‑145 (MQ=255)
                                                                                                gCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtc                                           <  2:149917/145‑1 (MQ=255)
                                                                                                gCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtc                                           >  1:109543/1‑145 (MQ=255)
                                                                                                      cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCga                                     <  2:333654/145‑1 (MQ=255)
                                                                                                                          tGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAAc                 <  2:278652/145‑1 (MQ=255)
                                                                                                                              gCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgct             >  2:78064/1‑145 (MQ=255)
                                                                                                                              gCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATGTAAAGCGTCGGTCGAGGAGGGCGGTAATCCCGAAcggct             >  2:408781/1‑145 (MQ=255)
                                                                                                                                  cacccAGGCCAACGAACGGCACATCAACAATCAACAAGTCAGGCTCCGCCATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTTGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTgcgc         <  1:402910/142‑1 (MQ=255)
                                                                                                                                         gCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGGGCGTGACAATGCAGGGAACATAGTGAACAGTGTGCGTAATGCTATTTACAGCGTCGACTGAGGAGGGCGGTATTACCGGACTGCTGCGCCAGCtgc  >  1:246831/1‑145 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGC  >  CP009273/788613‑788894

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