Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,665,299 G→A 7.6% G145S (GGC→AGC)  nuoI → NADH‑quinone oxidoreductase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,665,2990GA 7.6% 149.8 / 2.3 53G145S (GGC→AGC) nuoINADH‑quinone oxidoreductase subunit I
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (26/27)
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.26e-01

CATTTTCTGCGGCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATGGAATTCGCTTTCTACTTC  >  NC_002947/4665169‑4665446
                                                                                                                                  |                                                                                                                                                   
cATTTTCTGCGGCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTAc                                                                                                                                    >  6:237228/1‑148 (MQ=255)
 aTTTTCTGCGGCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACaa                                                                                                                                  >  4:225070/1‑149 (MQ=255)
         cGGCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAATCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTa                                                                                                                            <  7:331160/147‑1 (MQ=255)
          ggCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCg                                                                                                                         <  1:342803/149‑1 (MQ=255)
            ccTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCgtg                                                                                                                       <  3:225070/149‑1 (MQ=255)
             cTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCgtgt                                                                                                                      <  4:143057/149‑1 (MQ=255)
                   gaagaaGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCgg                                                                                                                 <  6:296042/148‑1 (MQ=255)
                    aagaagCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCggg                                                                                                                >  6:235747/1‑148 (MQ=255)
                    aagaagCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGt                                                                                                               >  6:353063/1‑149 (MQ=255)
                       aagCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATg                                                                                                            <  5:353063/149‑1 (MQ=255)
                               ccGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGc                                                                                                    >  8:242853/1‑149 (MQ=255)
                                      ccGCTATCCAGATGAATACGGATTTCGAAATGGCCGAGTTCAAGCGTCATTACCTGTTGTACGAGAAAGAAGATCTGCTGATCTCCTGCAACTGCAAGAAACCTGACTACAACTTATACATTGTTGAGGGTATGGCAATAGATGGc                                                                                                >  2:126337/1‑146 (MQ=255)
                                           aTCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGATAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGaa                                                                                        <  1:212549/149‑1 (MQ=255)
                                           aTCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGaa                                                                                        <  8:43283/149‑1 (MQ=255)
                                           aTCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGaa                                                                                        <  4:42947/149‑1 (MQ=255)
                                               aGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGc                                                                                    >  5:50965/1‑149 (MQ=255)
                                                gCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAAcc                                                                                                                                            >  3:233070/1‑92 (MQ=255)
                                                gCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAAcc                                                                                                                                            <  4:233070/92‑1 (MQ=255)
                                                gCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGct                                                                                   >  4:175941/1‑149 (MQ=255)
                                                gCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGc                                                                                    >  3:237428/1‑148 (MQ=255)
                                                      tCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCGTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGca                                                                              >  6:369064/1‑148 (MQ=255)
                                                      tCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGca                                                                              >  7:8135/1‑148 (MQ=255)
                                                        cGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGcac                                                                             <  3:175941/147‑1 (MQ=255)
                                                        cGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGca                                                                              >  6:249004/1‑146 (MQ=255)
                                                        cGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCAGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGca                                                                              >  7:266496/1‑146 (MQ=255)
                                                        cGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCAGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGca                                                                              >  3:334202/1‑146 (MQ=255)
                                                               cGAAATGGCCGAGTTTAAGCTTCAGGACCTTGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGTCAAGAACCCTGACTAGAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGTGCTCTGCACAGATCGaa                                                                     <  8:271293/148‑1 (MQ=255)
                                                                gAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAagc                                                                   >  4:131469/1‑149 (MQ=255)
                                                                      gCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAagccgagc                                                              <  7:242853/148‑1 (MQ=255)
                                                                            ttCAAGCGTCAGGAACTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACGGAACGAAGCAGAGCCGAtaa                                                        >  2:261541/1‑146 (MQ=255)
                                                                               aaGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGt                                                    >  1:157536/1‑149 (MQ=255)
                                                                               aaGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGt                                                    >  1:152221/1‑149 (MQ=255)
                                                                                       ggACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCt                                            >  4:294218/1‑149 (MQ=255)
                                                                                       ggACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCt                                            <  6:50965/149‑1 (MQ=255)
                                                                                              gtgtACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTccc                                     <  5:237228/149‑1 (MQ=255)
                                                                                              gtgtACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTc                                       <  5:235747/147‑1 (MQ=255)
                                                                                              gtgtACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTc                                       <  6:58441/147‑1 (MQ=255)
                                                                                                   cGAGAAAGAAGATCTGCTGATCTCCGGCCCCAGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATaa                                 <  4:334202/148‑1 (MQ=255)
                                                                                                   cGAGAAAGAAGATCTGCTGATCTCCGGCCCCAGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATaa                                 <  8:266496/148‑1 (MQ=255)
                                                                                                    gagaAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAg                                >  8:36610/1‑148 (MQ=255)
                                                                                                     agaaagaaGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGa                              >  3:131176/1‑149 (MQ=255)
                                                                                                      gaaagaaGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGa                              <  4:131176/148‑1 (MQ=255)
                                                                                                         agaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGa                          <  3:131469/149‑1 (MQ=255)
                                                                                                          gaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATg                                                                                                            <  7:16611/66‑1 (MQ=255)
                                                                                                          gaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATg                                                                                                            >  8:16611/1‑66 (MQ=255)
                                                                                                          gaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAg                                <  8:318534/142‑1 (MQ=255)
                                                                                                          gaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAg                                >  7:318534/1‑142 (MQ=255)
                                                                                                          gaagaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGa                          >  4:85033/1‑148 (MQ=255)
                                                                                                            agaTCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAg                       <  3:85033/149‑1 (MQ=255)
                                                                                                                   ctgATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATg                    <  8:8135/145‑1 (MQ=255)
                                                                                                                                ccGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATGGAATTCGCTTTCTACtt   <  2:157536/149‑1 (MQ=255)
                                                                                                                                ccGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATGGAATTCGCTTTCTACtt   <  2:152221/149‑1 (MQ=255)
                                                                                                                                 cGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATGGAATTCGCTTTCTACTTc  >  4:98645/1‑149 (MQ=255)
                                                                                                                                  |                                                                                                                                                   
CATTTTCTGCGGCCTGTGTGAAGAAGCGTGCCCGACCACCGCGATCCAGCTGACTCCGGATTTCGAAATGGCCGAGTTCAAGCGTCAGGACCTGGTGTACGAGAAAGAAGATCTGCTGATCTCCGGCCCCGGCAAGAACCCTGACTACAACTTCTACCGTGTTGCGGGTATGGCAATCGCTGGCAAGCCGAAGGGCTCTGCACAGAACGAAGCCGAGCCGATCAACGTGAAGAGCTTGCTCCCATAAGGACAGAAAGATGGAATTCGCTTTCTACTTC  >  NC_002947/4665169‑4665446

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: ‑