Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,676,624 T→C H427R (CAC→CGC)  pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,676,6240TC91.5% 215.7 / 6.1 72H427R (CAC→CGC) pntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (34/31);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TCTGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTG  >  NC_000913/1676490‑1676757
                                                                                                                                      |                                                                                                                                     
tCTGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGccc                                                                                                                                   >  2:1431626/1‑139 (MQ=255)
                     caaTAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAAc                                                                                                              <  1:210913/139‑1 (MQ=255)
                           tCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  2:264898/129‑1 (MQ=255)
                           tCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  1:264898/1‑129 (MQ=255)
                                       cgTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGGAGCGCCCAAGGAATTCTTCCGGCGCAACGCTTGCCATCCa                                                                                                  <  2:513977/133‑1 (MQ=255)
                                       cgTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGGAGCGCCCAAGGAATTCTTCCGGCGCAACGCTTGCCATCCa                                                                                                  >  1:513977/1‑133 (MQ=255)
                                       cgTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCaaa                                                                                            >  1:1562186/1‑139 (MQ=255)
                                           ggTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGa                                                                                        >  2:1503946/1‑139 (MQ=255)
                                                  gACATCAACGGTGTATGTAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGcc                                                                                               >  2:992358/1‑125 (MQ=255)
                                                  gACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGcc                                                                                               <  1:992358/125‑1 (MQ=255)
                                                          cGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGcc                                                                         <  1:1431626/139‑1 (MQ=255)
                                                            gtgtATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCAt                                                                       <  1:1584026/139‑1 (MQ=255)
                                                             tgtATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATc                                                                      <  1:19040/139‑1 (MQ=255)
                                                              gtATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  1:606849/1‑94 (MQ=255)
                                                              gtATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  2:606849/94‑1 (MQ=255)
                                                                 tGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcg                                                                  >  2:1690817/1‑139 (MQ=255)
                                                                      cgcgTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAt                                                             <  1:226841/139‑1 (MQ=255)
                                                                       gcgTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAtt                                                            >  1:908404/1‑139 (MQ=255)
                                                                           gCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACg                                                        >  1:1491766/1‑139 (MQ=255)
                                                                             gATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGcc                                                      <  2:1720769/139‑1 (MQ=255)
                                                                                aCATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACg                                                   <  1:1503946/139‑1 (MQ=255)
                                                                                    tCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCACCGCGTATTTACGCCACGGTGa                                               >  1:1025338/1‑139 (MQ=255)
                                                                                      cacaCCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGc                                             <  2:908404/139‑1 (MQ=255)
                                                                                         accacGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAgg                                          >  2:1721284/1‑139 (MQ=255)
                                                                                            acGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGc                                                                                >  2:295230/1‑98 (MQ=255)
                                                                                            acGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGc                                                                                <  1:295230/98‑1 (MQ=255)
                                                                                            acGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGa                                                                                        >  1:1433521/1‑90 (MQ=255)
                                                                                            acGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGa                                                                                        <  2:1433521/90‑1 (MQ=255)
                                                                                             cGTTATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATc                                                                                                    >  2:1116259/1‑77 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATc                                                                                                    <  1:1116259/77‑1 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATc                                                                                                    >  2:1527354/1‑77 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATc                                                                                                    <  1:1527354/77‑1 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCaaaaa                                                                                          >  1:700232/1‑87 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCaaaaa                                                                                          <  2:700232/87‑1 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACa                                      >  1:579550/1‑139 (MQ=255)
                                                                                             cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACa                                      >  1:444052/1‑139 (MQ=255)
                                                                                              gTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  1:1069536/1‑62 (MQ=255)
                                                                                              gTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  2:1069536/62‑1 (MQ=255)
                                                                                               taataaCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACAtt                                    <  2:579550/139‑1 (MQ=255)
                                                                                                     cggACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCATCGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCt                              <  2:1025338/137‑1 (MQ=255)
                                                                                                     ccGACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCt                              <  2:1491766/139‑1 (MQ=255)
                                                                                                       gACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  1:192270/53‑1 (MQ=255)
                                                                                                       gACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  2:192270/1‑53 (MQ=255)
                                                                                                       gACAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCa                            >  1:1227261/1‑139 (MQ=255)
                                                                                                        aCAACGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAg                           >  2:207727/1‑139 (MQ=255)
                                                                                                          aaCGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGtt                         >  2:1071411/1‑139 (MQ=255)
                                                                                                          aaCGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTcc                               <  1:1492064/133‑1 (MQ=255)
                                                                                                          aaCGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTcc                               >  2:1492064/1‑133 (MQ=255)
                                                                                                           aCGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTcc                               <  1:642462/132‑1 (MQ=255)
                                                                                                           aCGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTcc                               >  2:642462/1‑132 (MQ=255)
                                                                                                            cGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  2:1018738/48‑1 (MQ=255)
                                                                                                            cGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  1:1018738/1‑48 (MQ=255)
                                                                                                            cGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCa                                                                                                                >  1:1098292/1‑50 (MQ=255)
                                                                                                            cGCAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCa                                                                                                                <  2:1098292/50‑1 (MQ=255)
                                                                                                              cAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcgc                                                                 >  2:1410055/1‑95 (MQ=255)
                                                                                                              cAGGCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcgc                                                                 <  1:1410055/95‑1 (MQ=255)
                                                                                                                 gCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGcc                                                                         >  2:781403/1‑84 (MQ=255)
                                                                                                                 gCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGcc                                                                         <  1:781403/84‑1 (MQ=255)
                                                                                                                 gCCAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACtt                  >  1:1146427/1‑139 (MQ=255)
                                                                                                                  ccAGCGCGAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTc                 <  1:1690817/139‑1 (MQ=255)
                                                                                                                   cAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTcc                >  1:130802/1‑139 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGGGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  1:56878/38‑1 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  <  1:1289802/38‑1 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  2:56878/1‑38 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGCGCCCAAGGAATTCTTTCGgcg                                                                                                                  >  2:1289802/1‑38 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTg                                                <  2:1550400/104‑1 (MQ=255)
                                                                                                                      cgcgAAAACGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTg                                                >  1:1550400/1‑104 (MQ=255)
                                                                                                                           aaaCGGTGAAGCGCCGAAGGAATTCTTTCGGCGCAACGCTTGCCATCCa                                                                                                  <  1:951070/49‑1 (MQ=39)
                                                                                                                           aaaCGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCa                                                                                                  >  2:951070/1‑49 (MQ=255)
                                                                                                                           aaaCGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAt                                                             <  2:512005/86‑1 (MQ=255)
                                                                                                                           aaaCGGTGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAt                                                             >  1:512005/1‑86 (MQ=255)
                                                                                                                                 tGAAGCGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtg  >  1:739126/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TCTGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTG  >  NC_000913/1676490‑1676757

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