Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,969,465 A→G 12.0% W19R (TGG→CGG)  ubiG ← ubiquinone biosynthesis O‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,969,4650AG12.0% 123.6 / 8.5 50W19R (TGG→CGG) ubiGubiquinone biosynthesis O‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (2/4);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

CGCAGGGCCATGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGTTTTT  >  NC_002947/1969317‑1969604
                                                                                                                                                    |                                                                                                                                           
cgcAGGGCCATGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCccacca                                                                                                                                             <  7:310039/149‑1 (MQ=255)
    gggCCATGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCgg                                                                                                                                         <  6:157021/149‑1 (MQ=255)
       ccATGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTgg                                                                                                                                      <  2:244806/149‑1 (MQ=255)
         aTGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGc                                                                                                                                    <  1:37857/149‑1 (MQ=255)
            gCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCGCTGGTgggacc                                                                                                                                  <  1:180203/148‑4 (MQ=255)
             cTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCGCTGGTgggaccgc                                                                                                                                <  7:131142/149‑6 (MQ=255)
             cTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCGCTGGTgggac                                                                                                                                   >  4:97869/1‑144 (MQ=255)
             cTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCGCTGGTgggac                                                                                                                                   <  3:97869/146‑3 (MQ=255)
             cTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAgc                                                                                                                                >  7:296842/1‑149 (MQ=255)
                 gCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGGTGACGCGCAGCGGGTTGATTTCGTGGATCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCt                                                                                                                             <  7:32471/148‑1 (MQ=255)
                       gATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGaa                                                                                                                        >  8:305746/1‑147 (MQ=255)
                       gATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGaa                                                                                                                        <  7:305746/147‑1 (MQ=255)
                         tGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTg                                                                                                                    >  5:23058/1‑149 (MQ=255)
                          gccgccACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTgg                                                                                                                   >  1:297801/1‑149 (MQ=255)
                            cgccACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCAAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGc                                                                                                                  >  8:215991/1‑148 (MQ=255)
                            agccACCGCAACCCACGTCCAGCAGTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCAAGTTGACGCGCAGAGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGAGGTCCCACCAGCGGTGGGCCAGAGCTTCGAACTTGGc                                                                                                                  >  2:211121/2‑148 (MQ=255)
                             gccACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGa                                                                                                                >  3:133813/1‑149 (MQ=255)
                                 ccGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTc                                                                                                            >  2:100733/1‑149 (MQ=255)
                                       cccACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGcgcg                                                                                                      >  3:60881/1‑149 (MQ=255)
                                       cccACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGcgcg                                                                                                      >  5:223353/1‑149 (MQ=255)
                                        ccACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCTGGCCCACCGCTGGTGGGACcgcgaaagc                                                                                                                          >  8:312549/1‑123 (MQ=38)
                                        ccACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCTGGCCCACCGCTGGTGGGACcgcgaaagc                                                                                                                          <  7:312549/128‑6 (MQ=38)
                                          aCGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTc                                                                                                   <  6:223353/149‑1 (MQ=255)
                                                aGCACTTTCTTGCCGGCCAGGCAGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGtt                                                                                             >  7:50578/1‑149 (MQ=255)
                                                 gCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTg                                                                                            <  6:287768/149‑1 (MQ=255)
                                                   aCTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCt                                                                                          <  4:60881/149‑1 (MQ=255)
                                                     tttCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCa                                                                                        >  3:262024/1‑149 (MQ=255)
                                                        cTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCAtt                                                                                      <  2:177063/148‑1 (MQ=255)
                                                                 cAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCct                                                                            >  3:8474/1‑149 (MQ=255)
                                                                  aGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCctc                                                                           >  2:285905/1‑149 (MQ=255)
                                                                           cgcTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTc                                                                  <  1:285905/149‑1 (MQ=255)
                                                                                 tcgATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCaaaa                                                            <  1:100733/149‑1 (MQ=255)
                                                                                 tcgATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCaaaa                                                            >  6:193987/1‑149 (MQ=255)
                                                                                   gATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAAtt                                                          >  5:109479/1‑149 (MQ=255)
                                                                                    aTCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTg                                                         >  3:215917/1‑149 (MQ=255)
                                                                                       cAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTgcgc                                                      <  4:215917/149‑1 (MQ=255)
                                                                                               gcgcAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGta                                               >  7:155553/1‑148 (MQ=255)
                                                                                               gcgcAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAACTGCGCCGGAGta                                               >  8:122793/1‑148 (MQ=255)
                                                                                                                 gTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCcgc                            <  2:297801/149‑1 (MQ=255)
                                                                                                                     aGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAgg                        <  8:296842/149‑1 (MQ=255)
                                                                                                                        ggCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTc                     <  4:8474/149‑1 (MQ=255)
                                                                                                                                 tCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATc            >  3:123676/1‑149 (MQ=255)
                                                                                                                                   gCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGc          <  2:41057/149‑1 (MQ=255)
                                                                                                                                   gCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGc          >  5:183394/1‑149 (MQ=255)
                                                                                                                                     tttCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTg                                 >  6:247862/1‑124 (MQ=255)
                                                                                                                                     tttCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTg                                 <  5:247862/124‑1 (MQ=255)
                                                                                                                                        cgcgGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGt      <  4:123676/148‑1 (MQ=255)
                                                                                                                                         gcgGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGtt     <  6:109479/148‑1 (MQ=255)
                                                                                                                                          cgGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGtttt   <  5:193987/149‑1 (MQ=255)
                                                                                                                                           ggTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGttttt  <  4:133813/149‑1 (MQ=255)
                                                                                                                                           ggTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGttttt  <  8:155553/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CGCAGGGCCATGGCTTCGCTGAGGATGCCGCCACCGCAACCCACGTCCAGCACTTTCTTGCCGGCCAGGCTGACGCGCTCGTCGATCCAGTTGACGCGCAGCGGGTTGATTTCGTGCAGCGGCTTGAACTCGCTTTCGCGGTCCCACCAGCGGTGGGCCAGCGCTTCGAACTTGGCGATTTCGGCGCGGTCGACGTTGCTCATTGAACAGTCCTCTGAAACTTCGCAAAATTGCGCCGGAGTATACCCGAGCGCTTGGCCGCGAGGGTCGCTATAATCGCCGGTTTTT  >  NC_002947/1969317‑1969604

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