Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,676,071 G→A 8.8% A97V (GCG→GTG)  pntB ← pyridine nucleotide transhydrogenase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,676,0710GA 8.8% 130.6 / 2.3 45A97V (GCG→GTG) pntBpyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (24/17);  new base A (2/2);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTATCCGG  >  NC_000913/1675928‑1676207
                                                                                                                                               |                                                                                                                                        
tGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCAATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAg                                                                                                                                      <  2:308769/148‑1 (MQ=255)
    gAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCa                                                                                                                                 <  7:261848/149‑1 (MQ=255)
                             gTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATATCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGTCACCAg                                                                                                        >  2:80806/1‑149 (MQ=255)
                                 ccGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCTTGATGCAGATAGCTGTTAAAGCCATCAAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTc                                                                                                     <  1:267515/148‑1 (MQ=255)
                                 ccGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCt                                                                                                    <  8:149537/149‑1 (MQ=255)
                                  cGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTg                                                                                                   >  6:42285/1‑149 (MQ=255)
                                             aGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCACCAGACCCACGAAGCTACGCAGGATCGccacca                                                                                                         >  7:75636/1‑132 (MQ=255)
                                             aGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCACCAGACCCACGAAGCTACGCAGGATCGccacca                                                                                                         <  8:75636/132‑1 (MQ=255)
                                              gATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCACCAGACCCACGAAGCTATGCAg                                                                                                                    <  4:246674/120‑1 (MQ=255)
                                              gATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCACCAGACCCACGAAGCTATGCAg                                                                                                                    >  3:246674/1‑120 (MQ=255)
                                                tACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCAtt                                                                                     <  6:158281/149‑1 (MQ=255)
                                                   cGAGGAACACTTCCGTAAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTc                                                                                   >  8:274738/1‑148 (MQ=255)
                                                      ggAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACt                                                                               >  2:39738/1‑149 (MQ=255)
                                                        aaCACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACttt                                                                             >  4:29513/1‑149 (MQ=255)
                                                         acacTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACACACGAAGCTATGCATGATCGCCACCAGTTCTGGCATTTCGTTCATTTCAACAttc                                                                            >  6:79270/1‑149 (MQ=255)
                                                             ttCCGTCAGGTGCATATTGACCAGAATCGGTGCCATTCCCGCGGCATGATGCAGATAGCTGGTAAAGCCAACCAGCACTGCCGCCAGACCCACGCAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCg                                                                        >  2:198356/1‑149 (MQ=255)
                                                             ttCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCCCGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCACCTTTCTTCg                                                                        >  4:164645/1‑149 (MQ=255)
                                                               acGTCAGGTGAATATTGACCATAATCTGTGCCATTCGCGCGTCATGATGCAGATAGCTGTTAAAGCCAACGAGCAGTGCCGCCAGACCCACGAAGCTATGCAGGATCGCGACCAGTTCTGGCATTTCGGTCATTTCTACTTTCTTCGc                                                                       <  6:153948/147‑1 (MQ=255)
                                                                cGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGTCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCa                                                                     <  1:80806/149‑1 (MQ=255)
                                                                 gTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGc                                                                       >  8:105895/1‑146 (MQ=255)
                                                                       tGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGa                                                               >  7:90203/1‑148 (MQ=255)
                                                                           tatTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATAcc                                                           >  3:168956/1‑148 (MQ=255)
                                                                                         ggTGCCATTCCCGCGTCATGCTGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACGCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCAGTTCGGTCATTTCAACTTTCTTAGCCGGCCGGATACCAATGGCCccacccag                                            >  6:79955/1‑146 (MQ=255)
                                                                                               aTTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAgg                                                                                                                   <  8:8225/72‑1 (MQ=255)
                                                                                               aTTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAgg                                                                                                                   >  7:8225/1‑72 (MQ=255)
                                                                                                     gcgTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGc                                                                       <  5:222204/110‑1 (MQ=255)
                                                                                                     gcgTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGc                                                                       >  6:222204/1‑110 (MQ=255)
                                                                                                         cATTATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCATACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAg                            <  5:79270/149‑1 (MQ=255)
                                                                                                         cATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGc                               >  2:303685/1‑146 (MQ=255)
                                                                                                          atgatgCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAg                            <  4:168956/148‑1 (MQ=255)
                                                                                                          atgatgCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGa                           <  1:39738/149‑1 (MQ=255)
                                                                                                           tgatgCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGAt                          <  7:274738/149‑1 (MQ=255)
                                                                                                                cAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGAcc                                        <  2:330277/130‑1 (MQ=255)
                                                                                                                cAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGAcc                                        >  1:330277/1‑130 (MQ=255)
                                                                                                                cAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGc                     >  2:37912/1‑149 (MQ=255)
                                                                                                                        tGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTa             >  1:283672/1‑149 (MQ=255)
                                                                                                                            aaaGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTAccc          <  7:105895/148‑1 (MQ=255)
                                                                                                                              aGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCg         >  2:241669/1‑147 (MQ=255)
                                                                                                                               gCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACa                >  6:110409/1‑139 (MQ=255)
                                                                                                                               gCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACa                <  5:110409/139‑1 (MQ=255)
                                                                                                                               gCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTAt      >  6:251093/1‑149 (MQ=255)
                                                                                                                                  aaCCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTATCCg   >  4:127016/1‑149 (MQ=255)
                                                                                                                                  aaCCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTGCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTATCCg   <  5:251093/149‑1 (MQ=255)
                                                                                                                                   accaccAATGCCGCCAGACCAACGAATCTATGCAGGATCGCCACCAGATCTGGCATTTCGGTCATTTAAACATTCTTCGCCAGACCTATAACAATATCCCCACCAAACACCATCGCCATCAATATCCAGACAAAAAAACCCCAATCCg   >  4:80739/1‑148 (MQ=255)
                                                                                                                                   aCCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTATCCgg  >  8:89496/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCAGCCAACATTACCCGTATCCGG  >  NC_000913/1675928‑1676207

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