Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,910,100 C→T 12.1% D196N (GAC→AAC)  mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,910,1000CT12.1% 111.8 / 5.8 52D196N (GAC→AAC) mazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (3/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

ATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCT  >  NC_000913/2909964‑2910240
                                                                                                                                        |                                                                                                                                            
aTTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCAt                                                                                                                                        >  1:1033603/1‑143 (MQ=255)
      agcccgcgcgcgGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTtcctc                                                                                                                                  <  2:775082/137‑1 (MQ=255)
             gcgcgGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTGTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGt                                                                                                                           <  2:295894/143‑1 (MQ=255)
                cgGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTGCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTa                                                                                                                        <  1:561781/143‑1 (MQ=255)
                  gcagcaCTACGCTCCACTTCGCGAAAACGACGTTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTGAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGc                                                                                                                      <  2:229796/143‑1 (MQ=255)
                  gCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTGTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGc                                                                                                                      <  2:580035/143‑1 (MQ=255)
                       aaTACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCATACCTATTTGGCGAGCCAGATTACAGGTGGCAAACAGCAAACAACCCATTGCTTCCTCCAGTTTAGCCTGGt                                                                                                                 >  1:488075/1‑143 (MQ=255)
                        aTAGGCTCCTCTTCTCGAAAGCGACGCTCGAATTTTTCGTTCGCTTTTTGGAATGCGATTTCTGTTTTCGTCCTTAAATGGCGTGCCAGATTAACCGCGGCAAACAGCAGGTCTCCCATTTCCTCCTCCAGTTTAGCCTGGTc                                                                                                                <  2:406532/143‑1 (MQ=255)
                         tACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACGGCAGGTCCCCCATTTCCTCCTCCAGTTCAGCCTGGTCg                                                                                                               >  2:637614/1‑143 (MQ=255)
                              tCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGacaac                                                                                                          >  1:440639/1‑143 (MQ=255)
                                cACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAg                                                                                                        <  2:669533/143‑1 (MQ=255)
                                   ttCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAGTGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCCTTTCCTCCTCCCGTTCACCATGGTCGACGACGGCCt                                                                                                     >  2:702693/1‑143 (MQ=255)
                                       cgAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGCGCGACAACAGCCTGCcg                                                                                                 >  2:433390/1‑143 (MQ=255)
                                             cgacgCTCGAATTTTGCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTc                                                                                           <  2:924927/143‑1 (MQ=255)
                                                   tCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACAt                                                                                     <  1:349882/143‑1 (MQ=255)
                                                         gtttCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTTCCTAGATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtc                                                                               <  2:881015/142‑1 (MQ=255)
                                                               ttcgCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGTGAGCCAGATTAACCGTGGCAAACAGCAGGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACAACACCTCGTCGAt                                                                         >  2:355046/1‑143 (MQ=255)
                                                               ttcgCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                                         >  2:355047/1‑143 (MQ=255)
                                                               ttcgCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                                         >  2:355040/1‑143 (MQ=255)
                                                                   cTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTTGCAAACAGCCGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACACCAGCCCGCCGCGCACCGTCCATCCCCTCGTCGCtccgt                                                                     >  2:930409/1‑140 (MQ=255)
                                                                          cAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATCTCCTCCTCCAGTTTCGCCTGGTCCACAACAGCCTGCCCCGCTTCGTCCATCACCTCCTCGACCTCTTCGCAGa                                                              >  2:666955/1‑143 (MQ=255)
                                                                             tGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCATCGATCTATTCGAAGACtt                                                           >  1:804440/1‑143 (MQ=255)
                                                                                gATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTAt                                                        <  1:376398/143‑1 (MQ=255)
                                                                                gATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTAt                                                        <  1:645738/143‑1 (MQ=255)
                                                                                gATTGCTGCTTGCGTCCCTGAATGGCGAGCCAGATCAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTAt                                                        <  1:681846/143‑1 (MQ=255)
                                                                                aaTTCCTGCGTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCCTAGACTTTAt                                                        <  2:1053047/142‑1 (MQ=255)
                                                                                            cGTCCCGAAAGGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGAc                                            <  2:702789/143‑1 (MQ=255)
                                                                                                ccTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGACTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTGCGTACACGTAATCGACTACCGGACCAAg                                        >  2:412391/1‑143 (MQ=255)
                                                                                                 cTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGACTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCACCTACCGAAcccacc                                       >  2:321777/1‑139 (MQ=255)
                                                                                                 cTAAATGGCGAGCCAGATTAACCGGGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGCTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGc                                       <  1:39383/143‑1 (MQ=255)
                                                                                                  tAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCCAGcg                                      >  2:125923/1‑143 (MQ=255)
                                                                                                  tAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGcg                                      >  1:898718/1‑143 (MQ=255)
                                                                                                      tGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGcgtcgt                                  >  1:991354/1‑143 (MQ=255)
                                                                                                         cGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTGGATCTCTTCGTAGACTTTATCGAGTACCGGACGAAGCGTCGTccg                               >  1:630285/1‑142 (MQ=255)
                                                                                                               aGATTAATCTTGGCAAACAGCAGGTGCCTCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGaa                         <  1:501862/143‑1 (MQ=255)
                                                                                                                gATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAg                        <  2:787414/143‑1 (MQ=255)
                                                                                                                     acggTGGCTATCAGCAGGTCCCCCATTTCCTCCTCGAGTTTTGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGTTCCTATTGTCGTCCAATGGAAGCCAAc                   <  2:154278/140‑1 (MQ=255)
                                                                                                                          cccAAATAGCAGGCCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTgg              <  2:1033603/141‑1 (MQ=255)
                                                                                                                            cAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCACCGTTGgcg            >  1:878095/1‑143 (MQ=255)
                                                                                                                             aaaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCACTCGAAGCCAACCTTgaccc           >  1:979636/1‑139 (MQ=255)
                                                                                                                             aaaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCACCGTTGgcgc           >  1:1032792/1‑143 (MQ=255)
                                                                                                                             aaaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTATATCGACTACCGGACCAAGCGGCGTCCAATCGAAGACAACGCTGacgc           >  2:587886/1‑143 (MQ=255)
                                                                                                                                cagtagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAAc        <  1:220361/143‑1 (MQ=255)
                                                                                                                                 agcagGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTTGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACg       <  1:355040/143‑1 (MQ=255)
                                                                                                                                 agcagGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACg       <  1:355047/143‑1 (MQ=255)
                                                                                                                                 agcagGTTCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACg       <  1:355046/143‑1 (MQ=255)
                                                                                                                                 agcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACg       >  1:201585/1‑143 (MQ=255)
                                                                                                                                  gcagGTCCCCCATGTCCGCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGt      <  2:32813/143‑1 (MQ=255)
                                                                                                                                    agGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGttt    <  1:250280/143‑1 (MQ=255)
                                                                                                                                     ggTCCCCCAGTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCATCGCTCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTc   <  2:1054532/143‑1 (MQ=255)
                                                                                                                                      gTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAAGCGAAGCCAACGTTGGCGCACCGCTTcc  >  1:437586/1‑142 (MQ=255)
                                                                                                                                      gTCCCCCATTGCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCt  <  2:48200/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
ATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTCGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCT  >  NC_000913/2909964‑2910240

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