Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,391,574 A→T 16.0% A56A (GCT→GCA moaE ← molybdopterin synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,391,5740AT16.0% 60.5 / 7.2 25A56A (GCT→GCAmoaEmolybdopterin synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (9/12);  new base T (2/2);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

ACGATGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTGTTCTGTTTGTAT  >  CP000730/2391439‑2391711
                                                                                                                                       |                                                                                                                                         
aCGATGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTGGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTGCAATTTGTGCCAATTTCTTTTCAGCCAt                                                                                                                                      >  2:192890/1‑141 (MQ=255)
 cGATGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCAtt                                                                                                                                     >  1:50719/1‑141 (MQ=255)
 cGATGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCAtt                                                                                                                                     >  1:186618/1‑141 (MQ=255)
    tGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTTGCCATTGGaa                                                                                                                                 >  2:58757/1‑141 (MQ=255)
    tGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTTGCCATTGGaa                                                                                                                                 >  2:51741/1‑141 (MQ=255)
       ggTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAata                                                                                                                               <  2:76510/141‑1 (MQ=255)
       ggTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAata                                                                                                                               <  2:122740/141‑1 (MQ=255)
         tGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAAtata                                                                                                                             <  2:178871/141‑1 (MQ=255)
                          aTACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAat                                                                                                            <  2:262293/141‑1 (MQ=255)
                                    atCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGtttt                                                                                                  >  1:213095/1‑141 (MQ=255)
                                               gTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAg                                                                                       <  1:39806/141‑1 (MQ=255)
                                                  aTGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTcc                                                                                    >  1:54567/1‑141 (MQ=255)
                                                                       tACCCGTTATCGTTCCAGGCCCTTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTTAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGAGCACATGACCGGTaa                                                               >  2:44484/1‑141 (MQ=255)
                                                                            gTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAaaca                                                          >  1:71044/1‑141 (MQ=255)
                                                                              tATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAaacaac                                                        <  1:86331/141‑1 (MQ=255)
                                                                              tATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAaacaac                                                        <  2:217469/141‑1 (MQ=255)
                                                                                         gCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCtt                                             <  2:221340/141‑1 (MQ=255)
                                                                                                ttCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATGCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTc                                      <  1:230794/141‑1 (MQ=255)
                                                                                                                 cAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCAc                     <  2:186618/141‑1 (MQ=255)
                                                                                                                 cAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCAc                     <  2:50719/141‑1 (MQ=255)
                                                                                                                       gtgCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATAtt               <  1:29374/141‑1 (MQ=255)
                                                                                                                       gtgCCAATTTCTTTTTGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTg              <  1:51741/141‑1 (MQ=255)
                                                                                                                       gtgCCAATTTCTTTTTGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTg              <  1:58757/141‑1 (MQ=255)
                                                                                                                          ccAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTGtt            >  1:107591/1‑141 (MQ=255)
                                                                                                                                    ttCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTGTTCTGTTTGTAt  >  1:184877/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                         
ACGATGCGGTGAAGAAACCGCAATTAATACAGCGATATCTGAAATTTGTAATGGCCCTATTCTATGAACAATACTCGTTATCGTTCCAGGCCATTTTTCATTTATTTCATCTCCAATTTGTGCCAATTTCTTTTCAGCCATTGGAATATACGCTTCATATTCTAAATATTCCGTTTTGACGCCTTTAGTCCATTCGCGAACATGACCGGTAAAAACAACTACTGCACCTTGATATTCATTTATAGTGAATTCACGATATTGTTCTGTTTGTAT  >  CP000730/2391439‑2391711

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