Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,243,153 A→G 8.1% S195S (AGT→AGC)  PP_3720 ← NAD(P)H quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,243,1530AG 8.1% 135.3 / 2.8 49S195S (AGT→AGC) PP_3720NAD(P)H quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GAGAAAGGTCACGTCAGCGGGTACAACGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTCGCGCCCGCCCAGG  >  NC_002947/4243022‑4243300
                                                                                                                                   |                                                                                                                                                    
gagaAAGGTCACGTCAGCGGGTACAACGCTTCATCGAACTGCGCCAGGCGCGCGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATg                                                                                                                                     >  7:155531/1‑149 (MQ=255)
                gCGGGTACAACGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAAcggc                                                                                                                     <  6:79948/149‑1 (MQ=255)
                   ggTACAACGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAAcggcgg                                                                                                                   >  2:178776/1‑148 (MQ=255)
                    gTACAACGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGgca                                                                                                                 >  2:49967/1‑149 (MQ=255)
                       cAACGCTTCATCGAACTGCGCCAGGCGCGGTAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGgcagc                                                                                                               <  3:41660/148‑1 (MQ=255)
                        aaCGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACggcggc                                                                                                                  <  1:59769/144‑1 (MQ=255)
                        aaCGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACggcggc                                                                                                                  >  2:59769/1‑144 (MQ=255)
                        aaCGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCAc                                                                                                             >  3:190763/1‑149 (MQ=255)
                         aCGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCAcc                                                                                                            >  3:15291/1‑149 (MQ=255)
                                aTCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAc‑c                                                                                                    <  3:232243/149‑1 (MQ=255)
                                 tCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAc‑cc                                                                                                   <  4:15291/149‑1 (MQ=255)
                                   gAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGa                                                                                                 <  1:291613/149‑1 (MQ=255)
                                    aaCTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGa                                                                                                 >  7:250829/1‑148 (MQ=255)
                                          gcCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGc                                                                                          <  1:9886/149‑1 (MQ=255)
                                           ccAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGcc                                                                                         <  4:190763/149‑1 (MQ=255)
                                              ggCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGc                                                                                      <  1:49967/149‑1 (MQ=255)
                                                           cAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGaa                                                                         >  5:312004/1‑149 (MQ=255)
                                                           cAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGa                                                                          >  5:284835/1‑148 (MQ=255)
                                                           cAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGAGCATCAGCACTGATCTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGATC‑CCGACATAGGCCAGCGTGCACCGCAGaa                                                                         >  7:239597/1‑149 (MQ=255)
                                                             gTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATc                                                                       >  3:75805/1‑149 (MQ=255)
                                                                       gtcCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCa                                                                            >  3:132992/1‑134 (MQ=255)
                                                                       gtcCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCa                                                                            <  4:132992/134‑1 (MQ=255)
                                                                          cGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAAcggcgg                                                                                                                   <  2:84671/93‑1 (MQ=255)
                                                                          cGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAAcggcgg                                                                                                                   >  1:84671/1‑93 (MQ=255)
                                                                           gagaGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTc                                                         >  7:174335/1‑149 (MQ=255)
                                                                            agagCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGTCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCt                                                        <  6:312004/149‑1 (MQ=255)
                                                                            agagCTGCTCCAGACGCTGCCGCTAGTGCTGCAGAAAGTCCTGCCGCGCATCAGCACGGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCGCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCt                                                        <  8:239597/149‑1 (MQ=255)
                                                                              agCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTc                                                      <  8:174335/149‑1 (MQ=255)
                                                                               gctgctCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAACGCCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTc                                                      <  1:70304/149‑1 (MQ=255)
                                                                                         aCGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGt                                           >  7:150264/1‑149 (MQ=255)
                                                                                              gCCGGTATTGCTGCAGAAACTCCTGCCGCGCATCAGCACTGATGGACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATc                                      <  2:218216/149‑1 (MQ=255)
                                                                                               ccGGTATTGCTGCAGAAAGTCCTGCCGCGCAGCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATc                                      <  4:75805/148‑1 (MQ=255)
                                                                                                 ggTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATcgcg                                   >  4:289837/1‑149 (MQ=255)
                                                                                                         tgcAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGc                            >  5:207599/1‑148 (MQ=255)
                                                                                                          gcAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGcc                           >  1:254951/1‑148 (MQ=255)
                                                                                                           cAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCaa                         <  6:134388/149‑1 (MQ=255)
                                                                                                             gAAAGTCCTGGCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACTAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCTTGCCCCGCAGAATCGGCCGGAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAAc                       <  8:155531/149‑1 (MQ=255)
                                                                                                             gAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAAc                       <  6:284835/149‑1 (MQ=255)
                                                                                                                aGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACAtg                    <  8:150264/149‑1 (MQ=255)
                                                                                                                   ccTGCCGCGCATCAGCACTGATGTACGGCACATGACACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATCTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGc                 >  1:33279/1‑149 (MQ=255)
                                                                                                                    gtGCCGTACATCAGTGCTGATGCGCGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGc                                                                             >  6:120363/2‑88 (MQ=37)
                                                                                                                    gtGCCGTACATCAGTGCTGATGCGCGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGc                                                                             <  5:120363/87‑1 (MQ=37)
                                                                                                                     tGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATcgc                                    >  2:316340/1‑128 (MQ=255)
                                                                                                                     tGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATcgc                                    <  1:316340/128‑1 (MQ=255)
                                                                                                                      gCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTcg              >  8:112209/1‑149 (MQ=255)
                                                                                                                       ccgtacATCAGTGCTGATGCGCGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGcc                                 <  4:252234/124‑1 (MQ=255)
                                                                                                                       ccgtacATCAGTGCTGATGCGCGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGcc                                 >  3:252234/6‑129 (MQ=255)
                                                                                                                       ccGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTc               <  3:289837/147‑1 (MQ=255)
                                                                                                                                  cACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTCGCGCCCGCCCAgg  <  1:199447/149‑1 (MQ=255)
                                                                                                                                   aCTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTCgcgc           >  6:65213/1‑139 (MQ=255)
                                                                                                                                   aCTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTCgcgc           <  5:65213/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                    
GAGAAAGGTCACGTCAGCGGGTACAACGCTTCATCGAACTGCGCCAGGCGCGGGAACACCAGTGGCTGGTCGTCCGAGAGCTGCTCCAGACGCTGCCGGTATTGCTGCAGAAAGTCCTGCCGCGCATCAGCACTGATGTACGGCACATGCCACGCCACGAACGGCGGCAGCACCTCGAAC‑CCGACATAGGCCAGCGTGCCCCGCAGAATCGGCCGCAGCATGTCTTCCAGCGGCCCGTGAATCGCGCCGTCGCCAAACATGTGCTCGCGCCCGCCCAGG  >  NC_002947/4243022‑4243300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑