Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,957 G→A 11.5% intergenic (‑44/+29) 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,9570GA11.5% 121.8 / 7.9 54intergenic (‑44/+29)ridA/pyrIenamine/imine deaminase, redox regulated chaperone/aspartate carbamoyltransferase, PyrI subunit
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base A (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 = 6.18e-01

ATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATG  >  NC_000913/4470816‑4471099
                                                                                                                                             |                                                                                                                                              
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)
                                                                                                                                             |                                                                                                                                              
ATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATG  >  NC_000913/4470816‑4471099

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: