Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,391,841 T→C T401A (ACC→GCC)  nuoM ← NADH:quinone oxidoreductase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,8410TC98.1% 180.8 / ‑4.2 52T401A (ACC→GCC) nuoMNADH:quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (26/25);  total (26/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 = 1.89e-01

AGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCATAAAGCTGACCACAA  >  NC_000913/2391713‑2391968
                                                                                                                                |                                                                                                                               
aGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACATAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATcc                                                                                                                      >  2:17807/1‑140 (MQ=255)
aGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATcc                                                                                                                      >  2:894791/1‑140 (MQ=255)
  cgcgATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCa                                                                                                                    >  2:1099953/1‑140 (MQ=255)
      aTGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCg                                                                                                                >  2:582875/1‑140 (MQ=255)
      aTGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCg                                                                                                                >  1:447894/1‑140 (MQ=255)
            cATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCa                                                                                                          <  1:894791/140‑1 (MQ=255)
               cGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATcc                                                                                                                      <  1:680739/125‑1 (MQ=255)
               cGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATcc                                                                                                                      >  2:680739/1‑125 (MQ=255)
                 ccAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCa                                                                                                     >  2:695835/1‑140 (MQ=255)
                     cGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGa                                                                                                 >  2:612342/1‑140 (MQ=255)
                       aaTAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAAc                                                                                               >  2:585487/1‑140 (MQ=255)
                         tAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAg                                                                                             >  2:104663/1‑140 (MQ=255)
                            aCAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGa                                                                                          <  2:547223/140‑1 (MQ=255)
                            aCAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGa                                                                                          <  1:721652/140‑1 (MQ=255)
                                aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGt                                                                                      >  1:1104407/1‑140 (MQ=255)
                                    aaaGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTg                                                                                  >  2:129008/1‑140 (MQ=255)
                                        aCCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAg                                                                              <  2:545212/140‑1 (MQ=255)
                                        aCCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAg                                                                              <  2:447894/140‑1 (MQ=255)
                                          cAGCCCAAAGGGAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGcc                                                                            >  2:1022365/1‑140 (MQ=255)
                                             cccAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCAtt                                                                         >  2:139872/1‑140 (MQ=255)
                                                  aGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATc                                                                    <  1:695835/140‑1 (MQ=255)
                                                    gTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCtt                                                                  <  2:60008/140‑1 (MQ=255)
                                                    gTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCtt                                                                  <  1:1206828/140‑1 (MQ=255)
                                                      agagATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGc                                                                <  2:1013142/140‑1 (MQ=255)
                                                       gagaTAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCt                                                               >  2:205808/1‑140 (MQ=255)
                                                            aaCGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCcaca                                                          <  1:104663/140‑1 (MQ=255)
                                                            aaCGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCcaca                                                          >  2:797290/1‑140 (MQ=255)
                                                            aaCGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCcaca                                                          >  2:857119/1‑140 (MQ=255)
                                                              cGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGa                                                        >  2:1235811/1‑140 (MQ=255)
                                                               ggTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTGCCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGAc                                                       <  2:1104407/140‑1 (MQ=255)
                                                                   aTCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGcc                                                   >  2:945219/1‑140 (MQ=255)
                                                                    tCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGccc                                                  <  1:1040959/140‑1 (MQ=255)
                                                                    tCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGccc                                                  >  2:1040959/1‑140 (MQ=255)
                                                                            aCAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATgcg                                          <  1:585487/140‑1 (MQ=255)
                                                                                  tGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTc                                    <  1:205808/140‑1 (MQ=255)
                                                                                   ggAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcg                                   <  1:950487/140‑1 (MQ=255)
                                                                                    gAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgc                                  >  2:297002/1‑140 (MQ=255)
                                                                                     aaGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgcg                                 >  2:601016/1‑140 (MQ=255)
                                                                                     aaGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgcg                                 >  2:492813/1‑140 (MQ=255)
                                                                                        cTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGt                              <  1:612342/140‑1 (MQ=255)
                                                                                          gCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgc                                  >  2:159249/1‑134 (MQ=255)
                                                                                          gCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgc                                  <  1:159249/134‑1 (MQ=255)
                                                                                          gCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCgg                                >  1:681050/1‑136 (MQ=255)
                                                                                          gCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCgg                                <  2:681050/136‑1 (MQ=255)
                                                                                           ccGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgc                                  <  1:907798/133‑1 (MQ=255)
                                                                                           ccGAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTcgc                                  >  2:907798/1‑133 (MQ=255)
                                                                                             gAACAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGa                         <  1:601016/140‑1 (MQ=255)
                                                                                                cAGAATCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATAc                      <  1:1167037/140‑1 (MQ=255)
                                                                                                     tCATAAATTCGCCGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCa                 <  1:139872/140‑1 (MQ=255)
                                                                                                                ccGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATAc                      <  1:1112097/124‑1 (MQ=255)
                                                                                                                ccGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATAc                      >  2:1112097/1‑124 (MQ=255)
                                                                                                                ccGACGAAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCATAAAGCTGAcc      <  1:1235811/140‑1 (MQ=255)
                                                                                                                  gacgaAGTTACCGGCGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCATAAAGCTGACcac    <  1:857119/140‑1 (MQ=255)
                                                                                                                    cgaAGTTACCGGCGCCAGTCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGATCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCATAAAGCTGACCACaa  <  2:924538/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                               
AGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCAAAGAACAGCGACAGTGCTGGCAGCCATTTCATCTTGCTCCACAGACCGCCCATCATGCGCATGTCGCGGGTATGGATACGTTCATAAAGCTGACCACAA  >  NC_000913/2391713‑2391968

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