Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,958 G→T 11.2% intergenic (‑45/+28) ridA ← / ← pyrI enamine/imine deaminase, redox regulated chaperone/aspartate carbamoyltransferase, PyrI subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9580GT11.2% 120.2 / 9.7 54intergenic (‑45/+28)ridA/pyrIenamine/imine deaminase, redox regulated chaperone/aspartate carbamoyltransferase, PyrI subunit
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base T (2/4);  total (22/32)
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.30e-01

ATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGG  >  NC_000913/4470816‑4471100
                                                                                                                                              |                                                                                                                                              
aTCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATCATAgg                                                                                                                                                >  2:677661/1‑143 (MQ=255)
 tCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCAGTAAAAAATATGGCACATTATAggg                                                                                                                                               >  1:613620/1‑143 (MQ=255)
    gcccGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                                                                            <  2:137594/141‑1 (MQ=255)
     aCCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTAAGAAGTCTCGCTCATTAGAGGGAGcg                                                                                                                                           >  2:344195/1‑142 (MQ=255)
       cGGAGGTGATGGTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCTGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACACCGGTAAAAACTCTGGCTCATTATACGTAGGCAg                                                                                                                                         >  2:477178/1‑143 (MQ=255)
        ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                                                                        <  1:401960/143‑1 (MQ=255)
         gAGGTTATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                                                                       <  1:227206/143‑1 (MQ=255)
             tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                                                                   <  1:621660/143‑1 (MQ=255)
               atgatCATATTGCCCAGTGCTACGCCCTGTACGTAAGGACCGAGAGTTGGCGTGGCATTTTCCGTCGGTATAGGTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGTTCCCTATAATGAGCCAGACttttta                                                                                                                                 <  2:559495/143‑3 (MQ=255)
               aggctCATATTGCCCGGATCACCGCCTTGTACGTAAGGCCCGCTAGCTGCGGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttttta                                                                                                                                 <  2:370539/139‑3 (MQ=25)
                 gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTCTAGGGAGCCAGACTTttatta                                                                                                                               >  2:427560/1‑143 (MQ=255)
                 gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                                                               <  1:701521/143‑1 (MQ=255)
                     gtatTGCCCACGTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGCGCATTTTCCGTCGCGATAGTTTCGCTCATGATTTCTCCTTTATTACCGCGGTAATTAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                                                                           <  2:202085/142‑1 (MQ=255)
                                 tCAGCGCCCTGTACGTAGGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCcattaatg                                                                                                               <  2:880088/143‑2 (MQ=255)
                                  cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatta                                                                                                              >  1:129597/1‑142 (MQ=255)
                                   atcGCCCTGTACGTTAGGTCCGTTAGCTGCCGGTGCATTTTCCGGCGTGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                             <  2:16260/141‑1 (MQ=255)
                                   aaCGCCCTGTGCGTAATGACCGTTAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                             <  2:235816/143‑1 (MQ=255)
                                     cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                           <  2:7999/143‑1 (MQ=255)
                                       cccTGTACGTAAGGACCGATAGCCGCCGGTCCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                         <  2:46714/143‑1 (MQ=255)
                                        ccTGTACGTCAGGACCAATAGCTCCCCGTGCCCTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                                                                        <  1:902137/143‑1 (MQ=255)
                                            tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCGTTATAGGGAGCCAGACCTTTATTACTAAACGGAATCAACTGGCCAGaacca                                                                                                    >  2:679305/1‑143 (MQ=255)
                                               gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                                                                 >  2:308595/1‑143 (MQ=255)
                                                    gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                            <  1:317433/143‑1 (MQ=255)
                                                        gagtgCTGTCGGCGCAGTTTCCGGCGCGATCGTGTTGCTGTTGGGTGCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                        <  1:298945/139‑1 (MQ=255)
                                                            gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                                    >  1:748121/1‑143 (MQ=255)
                                                              cgCCGGTGCATTTGCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                                                                  <  1:786983/142‑1 (MQ=255)
                                                                      cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCGAGCACCACATTATGGGAAAAATCTTTTTccaag                                                                          >  1:722148/1‑139 (MQ=255)
                                                                        ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACATTTATCACCAACCGCAATTACTTGGCCAGCACCACATGATGGGAACACTATTATACACAgca                                                                        >  2:225025/1‑141 (MQ=255)
                                                                              gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGCGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCAGTTTCACAGTATTTACa                                                                  >  2:20429/1‑143 (MQ=255)
                                                                               tCGCGATAGGTTTGCTCATGATTTGTCCTTTATTACGAAGCTGAAAAGTCTGGCTCATTATAGGGAGCCGGACGCTTAAGACCAACCGCAATCAATTTGCCACCACCACATGATGGGAAAACTCTTTTTCACAGTATTGGCAt                                                                 >  2:239494/1‑143 (MQ=255)
                                                                                 gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGTCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                                                               <  1:308595/143‑1 (MQ=255)
                                                                                       gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                         <  1:849950/143‑1 (MQ=255)
                                                                                        ttttGCTCATGATTTCTCCTTTATTACCGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                        <  1:859295/143‑1 (MQ=255)
                                                                                         tttGCTCATGATTTCTCCTTTATTACAGCGTTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTCCCCACTATAATTAGTTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                       <  2:871844/143‑1 (MQ=255)
                                                                                          ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCATTTTCACAGTATTTGCATTTGAGCGCCat                                                      >  2:639470/1‑143 (MQ=255)
                                                                                               cttGATTTCGCCTTCCTTACAGCGTTAGAAAGTCTGGCTCATGATAGGGAGCCTCACTTTCATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                 <  2:280997/141‑1 (MQ=255)
                                                                                                    tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                            >  2:578142/1‑143 (MQ=255)
                                                                                                        gccgtttttaCAGCGGGAAAAGGCCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                        <  1:499254/136‑1 (MQ=255)
                                                                                                        gccTTTATTACAGCGGTAAAAAGGCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAATCTCGTGTTCTTAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                        <  1:179021/142‑1 (MQ=255)
                                                                                                         ccTTTATTACATCTGTAAAAAGTCTGCTGCATTCTAGGGAGCAAGGCTTTAGCCTCCATGCGGAATTCGTTCGCAAGCACCACATTATGGGAATACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                       <  2:918897/143‑1 (MQ=255)
                                                                                                           tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                                     <  2:209768/143‑1 (MQ=255)
                                                                                                             tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGTATTTGCGCGCGATAGCCTTGGCGCGGTTTCgg                                   >  1:585057/1‑142 (MQ=255)
                                                                                                                   aGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAACACTCTTTTTCCCAGTATTTGCATTTGAGCGCGCTCTCCTTGGCGCCTTTTCGCACGGCa                             >  2:96522/1‑143 (MQ=255)
                                                                                                                         aaaaaGTCTGGATCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTg                       <  1:221997/143‑1 (MQ=255)
                                                                                                                             aGTCTGGCTCATTATAGGGAGCCAGACTTTTTTTATCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATg                   >  2:334758/1‑143 (MQ=255)
                                                                                                                              gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCGTTGGCGCGTCTTCGCACGTCAACACTGGATGa                  >  1:802069/1‑143 (MQ=255)
                                                                                                                              gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGCTTTCGCACGGCAAAGCTGGATGa                  >  1:942628/1‑143 (MQ=255)
                                                                                                                                cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCCAAGCTGGCTGaaa                >  2:485391/1‑143 (MQ=255)
                                                                                                                                 tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTTTTTTTCACAGTATTTGCGTTTGAGCGCGCTATCATTGGCGCGGTTTCGCACGGCAAAGCTGGGTGAAAc               >  1:362515/1‑143 (MQ=255)
                                                                                                                                 tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                           >  1:708011/1‑131 (MQ=255)
                                                                                                                                 tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                           <  2:708011/131‑1 (MQ=255)
                                                                                                                                  ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc              <  1:862894/143‑1 (MQ=255)
                                                                                                                                   gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg             <  1:907164/143‑1 (MQ=255)
                                                                                                                                     tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATAAAACCGGt           >  2:911572/1‑143 (MQ=255)
                                                                                                                                      cATTATAGGGAGCCAGACTGTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCGTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGtt          <  1:6363/143‑1 (MQ=255)
                                                                                                                                       aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc         <  1:400153/143‑1 (MQ=255)
                                                                                                                                             gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGCGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg   <  1:751789/143‑1 (MQ=255)
                                                                                                                                              ggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg  >  1:759046/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGG  >  NC_000913/4470816‑4471100

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