Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 9.3% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC 9.3% 120.8 / 5.0 54N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (31/17);  new base C (3/2);  total (34/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 = 7.92e-02

GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAA  >  NC_000913/4472298‑4472577
                                                                                                                                           |                                                                                                                                             
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)
                                                                                                                                         aTTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGaa  >  2:650846/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAA  >  NC_000913/4472298‑4472577

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: