Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,916,815 A→G 8.1% G187G (GGT→GGC)  atpD ← F1 sector of membrane‑bound ATP synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,916,8150AG 8.1% 134.8 / 2.9 51G187G (GGT→GGC) atpDF1 sector of membrane‑bound ATP synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base G (2/2);  total (22/29)
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

GACCTTCGTCACGGAATTTCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGCCCTTAGCGAACG  >  NC_000913/3916670‑3916960
                                                                                                                                                 |                                                                                                                                                 
gACCTTCGTCACGGAATTTCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTAccc                                                                                                                                                >  1:16604/1‑149 (MQ=255)
    ttCGTCACGGAATTTCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGGCGGTCATTTCGTGGTAGAAGTCGTTACCCTCAc                                                                                                                                            <  2:479271/149‑1 (MQ=255)
       gTCACGGAATTTCTCAGCCCTGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCGCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGGCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAg                                                                                                                                         <  2:521998/149‑1 (MQ=255)
        tCACGGAATTTCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGt                                                                                                                                        >  2:970856/1‑149 (MQ=255)
            gtaaTTGCTCATCCATGGTCAGACCGGTCAGAGCAACGGGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGTATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGt                                                                                                                                    <  1:943699/147‑1 (MQ=255)
             gAATTGCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGtt                                                                                                                                   <  2:914958/149‑1 (MQ=255)
                 ttCTAAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGCCAATACACCTGGGATACTTTGTCGATAACGTTGGAGTCGGTCAAATCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCAcc                                                                                                                               >  2:875831/1‑149 (MQ=255)
                    tcAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATCTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTAc                                                                                                                            <  2:913447/149‑1 (MQ=255)
                     cAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCAATTCGTGGTAGAAGTCGTTGCCCTCACGAGTACGTTCACCTACg                                                                                                                           >  1:1008779/1‑149 (MQ=255)
                      aGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTTGAAGTCGTTCCCCTCACGAGTACGTTCACCTACGc                                                                                                                          >  1:893250/1‑149 (MQ=255)
                        ccATGGTCTGACCGGTCAGAGCAAGTGGCAGACGGTTTCCCGGCGTCTCGTTCATCTGGCCATACTCCTGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTAGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGccc                                                                                                                        <  1:875831/149‑1 (MQ=255)
                              tCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTGCCCTCACGAGTACGTTCACCTACGCCCGCAAac                                                                                                                  >  1:26126/1‑149 (MQ=255)
                               cAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAaca                                                                                                                 <  2:392767/149‑1 (MQ=255)
                                aGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAacac                                                                                                                >  1:15733/1‑149 (MQ=255)
                                  aCCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACag                                                                                                              >  1:886432/1‑149 (MQ=255)
                                   caggTCAGAGCATCGCGCAGAGGGTTTCCGGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTAGTCGATAACGTTGGAGTATGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACaga                                                                                                             <  2:270150/147‑1 (MQ=255)
                                            gCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCTTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGa                                                                                                    <  2:989633/149‑1 (MQ=255)
                                            gCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGa                                                                                                    <  1:294514/149‑1 (MQ=255)
                                             cAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAg                                                                                                   <  2:791494/149‑1 (MQ=255)
                                              aaCGCGCAGAGGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGt                                                                                                  <  1:767441/149‑1 (MQ=255)
                                               aCGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCCCCTACGCCCGCAAACACAGAGTAACCGGAGTg                                                                                                 >  1:609412/1‑149 (MQ=255)
                                                  cgcAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTGCCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTc                                                                                              <  2:26126/149‑1 (MQ=255)
                                                  cgcAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTc                                                                                              <  1:118558/149‑1 (MQ=255)
                                                     aGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGAt                                                                                           >  2:153138/1‑149 (MQ=255)
                                                           tttCCCGTCGGCTCGTTAATCTGGCCATACACCAGGGATACTTTGTCGATAACCTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGAt                                                                                     <  2:217097/149‑1 (MQ=255)
                                                              cccGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATCCCGGTGGAGTAGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTACCCGGAGTGCTCGATCGCGATGtt                                                                                  <  1:250683/149‑1 (MQ=255)
                                                                cggcggCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCAAGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTAc                                                                                >  2:642577/1‑149 (MQ=255)
                                                                cggcggCGCGTTCATCTGGCCATACACCCGGGATACTTTGTCGATTACGTTCGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTAc                                                                                <  1:863226/149‑1 (MQ=255)
                                                                       tCGTTCTTCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAg                                                                         >  1:351977/1‑149 (MQ=255)
                                                                         gTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGGCATTTCGTGGTAGAAGTCGTAACCCTCACTAGTACGTTCACCTACGCCCGCAACCACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCt                                                                       <  2:803966/149‑1 (MQ=255)
                                                                               cTGGACATACACCAGGGATACTTGGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCcatca                                                                 <  1:851321/149‑1 (MQ=255)
                                                                                 ggCCATACACCAGGGATACTTTGTCGAAAACGTTGGAGTCGGGCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATg                                                               <  2:31061/149‑1 (MQ=255)
                                                                                  cacATACACAAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCACACACAGAGTAACCGGAGTTCTCGATCGCGATGTTACGAATGCTCTCAAGCATGt                                                              >  2:112510/3‑149 (MQ=255)
                                                                                    cATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATGTCGTGGTAGCAGTCGTTACCCTCACGAGTACGTTCCCCTAAGCCCGCAAACTCAGAGTACGCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATAATGttt                                                            >  1:212017/1‑149 (MQ=255)
                                                                                         acCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTTCCCTCACGAGTACGTTCACCTACGCCCGCAATCACAGAGTTACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCCTGGTTACGGt                                                       >  2:622104/1‑149 (MQ=255)
                                                                                         acCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGt                                                       >  2:508601/1‑149 (MQ=255)
                                                                                              ggATACTTTGTCGCTAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACACTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTAc                                                  <  2:604850/149‑1 (MQ=255)
                                                                                              ggATACTTTGTCGATAACGTTGGAGTCGGGCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTAc                                                  <  1:377631/149‑1 (MQ=255)
                                                                                               gATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATAATGTTTACGGTTTTAcc                                                 <  2:834078/149‑1 (MQ=255)
                                                                                                 tagtttGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTa                                               <  1:373644/146‑1 (MQ=255)
                                                                                                 tACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTa                                               <  2:774897/149‑1 (MQ=255)
                                                                                                        tCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGc                                        >  1:516501/1‑149 (MQ=255)
                                                                                                          gATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGcac                                      >  1:1993/1‑149 (MQ=255)
                                                                                                               cGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGATCGATAGCGTTGTTACGAATGAGCTCCATCAGGTTTACGGTTTTACCTACACCCGCTCCACCg                                 >  2:933203/1‑149 (MQ=255)
                                                                                                                   ggAGGCGGGCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACa                             <  2:167676/149‑1 (MQ=255)
                                                                                                                                  gTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTAACCCCGCACCACCGAACAGACCAACTTTACCGc              >  1:569006/1‑149 (MQ=255)
                                                                                                                                   tGGTATAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGTGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGcc             <  2:138872/149‑1 (MQ=255)
                                                                                                                                   tGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGcc             >  2:467823/1‑149 (MQ=255)
                                                                                                                                             cGTTGCCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGCCCTTAGCGAAc   <  2:1008779/149‑1 (MQ=255)
                                                                                                                                              gTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCTCCACCGAACAGACGAACTTTACCGCCATTAGCGAACg  >  1:735936/1‑149 (MQ=255)
                                                                                                                                              gTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGCCCTTAGCGAACg  <  2:16604/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GACCTTCGTCACGGAATTTCTCAGCCATGGTCAGACCGGTCAGAGCAACGCGCAGACGGTTTCCCGGCGGCTCGTTCATCTGGCCATACACCAGGGATACTTTGTCGATAACGTTGGAGTCGGTCATTTCGTGGTAGAAGTCGTTACCCTCACGAGTACGTTCACCTACGCCCGCAAACACAGAGTAACCGGAGTGCTCGATCGCGATGTTACGAATGAGCTCCATCATGTTTACGGTTTTACCTACACCCGCACCACCGAACAGACCAACTTTACCGCCCTTAGCGAACG  >  NC_000913/3916670‑3916960

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