Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,434 G→A 9.6% R34* (CGA→TGA)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4340GA 9.6% 126.2 / 5.8 52R34* (CGA→TGA) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base A (3/2);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

GGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGA  >  NC_000913/4472292‑4472576
                                                                                                                                              |                                                                                                                                               
gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCATAAAAAAGCCCCTCg                                                                                                                                                 >  1:680106/1‑143 (MQ=255)
      gCTTTCATTTCCGCCGCTGTCGCCAGCACCAGGTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCCGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAg                                                                                                                                           <  2:588878/143‑1 (MQ=255)
         ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgggg                                                                                                                                        >  1:761607/1‑143 (MQ=255)
           cAGTTTCGCCGCTGTGGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCt                                                                                                                                      <  2:610709/143‑1 (MQ=255)
              tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCCggg                                                                                                                                   >  1:292298/1‑143 (MQ=255)
              tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCCCCTGACGCCTGGGCAAGAAAAAGCCCCTCGATTGGGGGGCgggg                                                                                                                                   >  1:615599/1‑143 (MQ=255)
                tCGCCGCGGTCGCCAGCACCAGAGTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                                                                 <  2:569819/143‑6 (MQ=255)
                 cgccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAt                                                                                                                                >  1:676334/1‑143 (MQ=255)
                     gcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTCAGGTCGTTTATGCAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                            <  1:59309/143‑1 (MQ=255)
                        ggcGCCCGTACCAGATTAAGGGCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                         <  2:776614/141‑1 (MQ=255)
                              agaaacAGATTAAGGTGATCGCGACGAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCGTTTATCTCCGGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                                                                                   <  2:690760/138‑1 (MQ=255)
                                          aGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAACACgg                                                                                                       >  2:152856/1‑143 (MQ=255)
                                            gTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCACCTCGATTGAGGGGGTGGGAATGGGAGAGCAACGGGAAGAACAACGGCa                                                                                                     >  2:433150/1‑143 (MQ=255)
                                            gTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCAATCTTTTAGCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCACCGCAACGaaaaaggcca                                                                                                     >  2:262235/1‑138 (MQ=255)
                                                   ggatgaaGGGCGGTTATGGAAATGATATGTGTCTGATATAGTGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                                                                              <  2:352474/138‑1 (MQ=255)
                                                           gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGGGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                     >  1:567659/1‑143 (MQ=255)
                                                               tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGAGGCCTGGGCAAAAAAACGCCCCTCGATTGAGGGGCTGGGGATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAg                                                                                 >  1:693474/1‑143 (MQ=255)
                                                                  aTGGAAAGGATGTGGCTCTGATATAGGGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                                                                              <  2:209198/143‑1 (MQ=255)
                                                                     gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                                                           >  2:701327/1‑143 (MQ=255)
                                                                       aaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCCGCGTCTTTTTTCAGACGCGGTa                                                                         >  1:35033/1‑143 (MQ=255)
                                                                         ttatatGTCTCTGGTATAGCCGATTAGCCATCTTTTATCTCCTCACGCCTGCGCAAGAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGa                                                                      <  2:337471/141‑1 (MQ=255)
                                                                           atatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCATGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAGGCGGTAAGAc                                                                     >  2:127540/1‑143 (MQ=255)
                                                                            tatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACa                                                                    >  1:614294/1‑143 (MQ=255)
                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                                                            <  2:676334/143‑1 (MQ=255)
                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                                                            >  1:352504/1‑143 (MQ=255)
                                                                                     gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAGATGGCGa                                                           >  1:154160/1‑143 (MQ=255)
                                                                                        ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAaca                                                        >  2:434181/1‑143 (MQ=255)
                                                                                         taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGACCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac                                                       >  2:120324/1‑143 (MQ=255)
                                                                                           gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACGAAATGTCGAacacac                                                     >  1:511248/1‑143 (MQ=255)
                                                                                               aTTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa                                                 >  1:735579/1‑143 (MQ=255)
                                                                                                  aGCCATCTTTTATCTCCTGACGCCTGGGCAAACAAAGCCCCTCGATTGAGGGGCTGGGACTGCGTGTATAACGTTAATAGCACTATACTTCTTGCTTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCAt                                             <  1:421688/143‑1 (MQ=255)
                                                                                                      aTCTTTTATCGCCTGACTGCGGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAAGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAAATCATGTCGAACACACTGAACCATACa                                          <  2:800828/143‑1 (MQ=255)
                                                                                                          tttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcct                                      >  2:387216/1‑143 (MQ=255)
                                                                                                             aTCACCAGACGCATGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                   <  1:629964/143‑1 (MQ=255)
                                                                                                                  cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGGGGGGCTGGGGATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaa                              >  1:821046/1‑143 (MQ=255)
                                                                                                                    gACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAACCTCCCGGCAAAt                            >  2:829600/1‑143 (MQ=255)
                                                                                                                       gCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                         >  1:12428/1‑143 (MQ=255)
                                                                                                                        ccTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAATAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTc                        <  2:446820/143‑1 (MQ=255)
                                                                                                                          tGGGCAAAGACAAGCCCCTTGATGGAGGGACTGGGAATGGGTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCGTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCg                      <  2:367048/143‑1 (MQ=255)
                                                                                                                           gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGACCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCgg                     <  1:120324/143‑1 (MQ=255)
                                                                                                                                aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                <  2:352504/143‑1 (MQ=255)
                                                                                                                                 ttcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTACGGCGCGGTAAGACAAAATCTCGAACACACTGAACCAGCCATCCTCCCGGCAAGTTGTCCGGCCTtat               >  2:365522/6‑143 (MQ=255)
                                                                                                                                    ccaGCGCCTCAATCGAGGGGCTTGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCACCCGGCAAACTGTCGGGCATTATAc             >  2:757790/3‑143 (MQ=255)
                                                                                                                                    aaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACt            <  2:292298/143‑1 (MQ=255)
                                                                                                                                     caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGTGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc           >  1:240256/2‑143 (MQ=255)
                                                                                                                                     aaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAACCGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATAAATCATCCCGGCAACTTGTCCGGCACTCTACTc           >  2:624622/1‑143 (MQ=255)
                                                                                                                                          ccgcGATTGAGGGGCTGGGAGTGGGTGATCAACGGGAAGAAAAGCGGCAGGCCAGCGTCTCTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt      <  2:417221/140‑1 (MQ=255)
                                                                                                                                          ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATCCTCATCGt      >  1:171434/1‑143 (MQ=255)
                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc     >  1:808264/1‑143 (MQ=255)
                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTAATAGTc     >  1:149081/1‑143 (MQ=255)
                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGAGTCTTTTTTTAGACGCGGTAAGACAAGATGTCGAACCCACTGAACCATACATCCTCCCGGCACATTGTCCGGCATTAGATTCAtgttc     >  1:39163/1‑139 (MQ=255)
                                                                                                                                            tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCCGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACAAACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCACCGACa    >  2:96658/1‑143 (MQ=255)
                                                                                                                                            ccGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa    <  1:720394/142‑1 (MQ=255)
                                                                                                                                              gATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTCTACTCATCGTCAGa  >  2:551951/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGA  >  NC_000913/4472292‑4472576

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: