Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,113 G→T H226Q (CAC→CAA sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,1130GT96.5% 154.3 / NA 57H226Q (CAC→CAAsthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (26/29);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATA  >  NC_000913/4159975‑4160251
                                                                                                                                          |                                                                                                                                          
ggTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGAt                                                                                                                                        <  2:967944/143‑1 (MQ=255)
    cGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGga                                                                                                                                    <  2:892931/143‑1 (MQ=255)
         gTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCCCTGTTCCCGAATTGATAGgagagag                                                                                                                               >  1:251519/1‑143 (MQ=255)
                tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTCGTGACGAATCACTACGCCCCTGTTCCAGAAGTGATAGGAGAGAGAATCTGa                                                                                                                        <  2:603410/143‑1 (MQ=255)
                     caggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGGCAtct                                                                                                                   >  1:996416/1‑143 (MQ=255)
                           gTCAGCTTTCAGTTTTTTACCCGACTTAAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGAAATCTATTGaa                                                                                                             >  2:466974/1‑142 (MQ=255)
                           gTCAGCTTTCACTTTCTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGAt                                                                                                             <  1:830743/143‑1 (MQ=255)
                              aGCTTTCAGTTTTTTACCCGACTGCAGATGCATGAGCACACCAGCGTCACAGCCTTCGATCTTCTCGTACGCTTCGTTGGGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCga                                                                                                          <  1:472665/143‑1 (MQ=255)
                               gtgggacgttgttTACCCGTCTGCAGTTGCATGATCGCACCGTCGGCACAGCCTTCGAGCTTCTCGTACTCTTGGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACAGCTCTTGATCtag                                                                                                         <  2:231550/136‑1 (MQ=255)
                                     ggTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGc                                                                                                   <  2:961848/142‑1 (MQ=255)
                                     aGTTTTTTACCCGACTTCAGATCCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGc                                                                                                   >  1:608140/1‑143 (MQ=255)
                                        ttttttCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCcag                                                                                                <  2:77841/143‑1 (MQ=255)
                                         ttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCCCTGTTCCAGAATTGATAGGAGCGAGAATCTGACACCTCTTGCTCGAGCACCGCcagc                                                                                               >  1:443049/1‑143 (MQ=255)
                                               ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATCGCAGAGAGAATCTGGCTTCCCTTGATCGAGAAATGCGAGCAGGCGa                                                                                         >  1:98248/1‑143 (MQ=255)
                                               ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                         >  2:16693/1‑143 (MQ=255)
                                                cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGAt                                                                                        <  1:749475/143‑1 (MQ=255)
                                                 gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTAATAGGAGAGAGAATCTGACAGATCTTGATCGAGAAATGCCAGCAGGCGATc                                                                                       >  1:845955/1‑143 (MQ=255)
                                                    ttCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgcg                                                                                    >  1:55240/1‑143 (MQ=255)
                                                     tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                   <  1:69635/143‑1 (MQ=255)
                                                     tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                   <  1:701864/143‑1 (MQ=255)
                                                        gATGCGTGATCACACCATCGTCACAGCCTTCGATCTTCTCGTAATCTTTGGTGTGACGAGTCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtg                                                                                <  1:43380/143‑1 (MQ=255)
                                                         aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGGCACAGTTACAGAATTCATAGGAGAGAGAATCTGACACCTCTTGATCGAGAAATGCCAGCAGGCGGTCGCGGgtgt                                                                               >  2:624118/1‑143 (MQ=255)
                                                             aTGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGGTGAt                                                                           >  1:254727/1‑143 (MQ=255)
                                                                       gctCGTCCGAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATGTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                 <  1:417097/141‑1 (MQ=255)
                                                                         tCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                               <  1:910936/143‑1 (MQ=255)
                                                                         tCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                               <  1:277036/143‑1 (MQ=255)
                                                                          cgccccaGCCCTCGATCTTCTCGCTCTCTTCGTTGTGACGAATCACGACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTATTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                              <  1:144325/138‑1 (MQ=255)
                                                                          cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCCGGCGATCGCGGGTGTTGGTCAGATCCACtgtt                                                              >  1:617107/1‑143 (MQ=255)
                                                                           gtggcaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                             <  1:78309/139‑1 (MQ=255)
                                                                             cacaGCCGTCGGTGTGCTCGTACTCTTCGTTGGGACGAATCACGACGCCACTGTTCCAGAATGGATAGGCGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCTCTTTTACa                                                           <  2:1012805/143‑1 (MQ=255)
                                                                               caGCCTGCGATCTTCTCGGTCTCTTCGGCGGGGCTGATCCCTACGCCACTGCTCCAGAATTGATAGGAGATAGAATCTGACAGCTCTGGAGCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                         <  2:1041835/143‑1 (MQ=255)
                                                                                   cTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGTGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                     <  2:352112/143‑1 (MQ=255)
                                                                                    ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCGGAATTGATAGGAGAGGGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCGGAGCGACTGTTACAGCCATAc                                                    >  1:445547/1‑143 (MQ=255)
                                                                                    ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                    >  2:1137859/1‑143 (MQ=255)
                                                                                    ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGGTGATCAGATCCCCTTTTACATCCAtgc                                                    >  1:563106/1‑141 (MQ=255)
                                                                                    ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGACAATCTGACATCTCTTGATCGAGGAATGCCAGCAGGTGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                    >  2:580793/1‑143 (MQ=255)
                                                                                     tCGATCTTCTTGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAcc                                                   >  1:559746/1‑143 (MQ=255)
                                                                                             ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCCGATCCACTTTTACATCCATACCGCGGAAgg                                           >  1:204165/1‑142 (MQ=255)
                                                                                               cGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGACCCACTTTTACATCCATACCGCGGAAGCTc                                         >  1:446228/1‑143 (MQ=255)
                                                                                                    tctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                    <  1:91690/143‑1 (MQ=255)
                                                                                                      ttCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                  >  2:193426/1‑143 (MQ=255)
                                                                                                      ttCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCAACTTTTACATCCATACAGCGGAAGATCGACGCat                                  >  2:573323/1‑143 (MQ=255)
                                                                                                       tCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGAGTGTAGATCAGATACACTAATACAACAATACCGCGGAAGATCGACGCata                                 >  1:1054296/1‑143 (MQ=255)
                                                                                                       tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGCGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                 <  2:520133/143‑1 (MQ=255)
                                                                                                             tgACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATAGAGAAATGCAAGCAGGCGATCGCGGGTGTTGATCAGAACAACTTTTGCATCCATACCGCGGAAGATCGAAGCATAttacca                           >  2:340281/1‑139 (MQ=255)
                                                                                                               gcGAATCACTACGCCCCTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                         <  2:996416/142‑1 (MQ=255)
                                                                                                               aCGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                         <  2:579585/143‑1 (MQ=255)
                                                                                                               aCGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                         <  1:812338/143‑1 (MQ=255)
                                                                                                                cGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                        <  1:1103938/143‑1 (MQ=255)
                                                                                                                  aaTCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCATGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa                      <  2:743012/143‑1 (MQ=255)
                                                                                                                    gcACTACGTCACTGTTCCAGAATTGATAGGACAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGCTGATCAGATCCACTTTTCCCTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                    <  1:231834/142‑1 (MQ=255)
                                                                                                                     cACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCa                   <  1:16693/143‑1 (MQ=255)
                                                                                                                        tACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGCTCACTc                >  2:1039690/1‑143 (MQ=255)
                                                                                                                            ccACTGTTTCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGc            <  2:1045983/143‑1 (MQ=255)
                                                                                                                                tGTTCCATAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCa        <  1:1046771/143‑1 (MQ=255)
                                                                                                                                      aGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATGCCGCGGAAGATCGAGGCATATTCACAGCGGATCACTCCAGCACCATAGATa  >  1:533543/1‑143 (MQ=255)
                                                                                                                                      aGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGaca  >  1:994433/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATA  >  NC_000913/4159975‑4160251

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: