Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 8.1% 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 8.1% 114.3 / 2.3 50N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (2/2);  total (27/23)
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.43e-01

TTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGC  >  NC_000913/4472296‑4472579
                                                                                                                                             |                                                                                                                                              
ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTg                                                                                                                                               >  2:766389/1‑143 (MQ=255)
ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTg                                                                                                                                               >  2:118848/1‑143 (MQ=255)
    tttCAGTTTCGCCGCTGTCGCCAGCCCCAGATTCAGGTCATCGCGACCAAAGTCGTTTATGGAAATGATATGTTTCAGATATAGCGGATGAGCCATCTTTTATCTCCTGACGCCTTGGCAAGAAAAAGCCCCTCGCTTGAggg                                                                                                                                           >  2:902282/1‑143 (MQ=255)
     ttCAGTTTCGCCGCTGTCGCCAGCACTAGATTAAGGTTATAGCGACTCAGGTCGTTTATGGACATGATGTGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCATGACGCCTGGGCAAAAAACAGCCCCTCGATCTAgggg                                                                                                                                          >  2:286551/1‑143 (MQ=255)
     ttCAGTTTCGCCGCTGTCGCCAGCACCAGCTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGAGTAGCCATCTTTTATCTCCTGACGCCTGGGCACAAAAACGCCCCGCGATTGAgggg                                                                                                                                          >  1:283483/1‑143 (MQ=255)
     ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAgggg                                                                                                                                          <  1:61117/143‑1 (MQ=255)
            tCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGACATCTTTTATATCCTGACGACTGGGCAAAAAAAAGCCCCTGGATTGCGGGGCTGGGaa                                                                                                                                   >  2:938418/1‑143 (MQ=255)
            tCGCCACTGTCGCCAGCACCAGATTGATGTCATCGCGACTAAGGTCGTTTATGGAAAGGATATGTTTCTGATATCGCGGATTAGCCATCTTTTATCTCCGGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                                                                   <  2:871250/143‑6 (MQ=255)
               ccgcTGTCGCCAGCACCAGATCACGCTCATCGCGCCTCAGGTCGTTTCTCGCAACGCTCTGTTTCTGCTACCGCGGACTAGCCATCCTTTACTCCCTGCCGCCTGGGCACAGAAGAGGCCCACTCTTGACGGCCTggatatgc                                                                                                                                >  2:147317/1‑137 (MQ=255)
                   tGTCCCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATTTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                            <  1:40927/143‑1 (MQ=255)
                   ggTCGCGAGCAGCAGGCTAGGGTCATCGCGACTAAGGTTGGTTGTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                            <  2:57612/142‑1 (MQ=255)
                         cAGCACCAGATTAAGGGCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                      <  1:203099/143‑1 (MQ=255)
                            cACCAGTTTCAGGCCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                                   <  1:835440/143‑1 (MQ=255)
                             aCCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGCGTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                  <  1:516291/143‑1 (MQ=255)
                               cAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAGCGGGTGCTCAACGCgaaga                                                                                                                >  2:767452/1‑143 (MQ=255)
                                           aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                                                                    <  2:393374/143‑1 (MQ=255)
                                           aTCGCGACGAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                                                                    <  2:48683/143‑1 (MQ=255)
                                            tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGAGTGCTATAAGGTTTGGAAAACGGCAGGcc                                                                                                   >  1:725728/1‑143 (MQ=255)
                                                       gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCACAAAAAGGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                        >  1:793648/1‑143 (MQ=255)
                                                       gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGCAGAAACACGGCAGGCCAGCGTCttttt                                                                                        >  1:701047/1‑143 (MQ=255)
                                                               tGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgc                                                                                <  1:369853/143‑1 (MQ=255)
                                                                ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCAAGCATCTTTTTTCAGAcgcg                                                                               >  2:700429/1‑143 (MQ=255)
                                                                  aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCCGt                                                                             >  2:820129/1‑143 (MQ=255)
                                                                     tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGCATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                          >  2:418475/1‑143 (MQ=255)
                                                                             ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                                                                  <  1:766389/143‑1 (MQ=255)
                                                                             ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                                                                  <  1:118848/143‑1 (MQ=255)
                                                                              tCTGATATAGCGGACTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                 <  1:244757/143‑1 (MQ=255)
                                                                               ccgATATGGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCGGGGAATGGGGGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                <  1:341231/141‑1 (MQ=255)
                                                                                  atgtaGCAGATTAGCCCTCTTCCGTGTCTTGACCCGGGGGCAAAAAAAAGCCCCTCGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                             <  2:862768/140‑1 (MQ=255)
                                                                                      aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                         <  1:358036/143‑1 (MQ=255)
                                                                                               gCCACCTTTCAACTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCAt                                                <  1:333558/143‑1 (MQ=255)
                                                                                                ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGCGGGGCCGGGGATGGGTGATCAACGGGAAGAAAAACGGCAGGCCCGCGTCTTTGTGCAGACGCGGTAAGACAAAATGTCAGACGCCCTGAACCATa                                               >  2:897465/1‑143 (MQ=255)
                                                                                                 cATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGGCGCGGTAGGACAAAATGTCGGACACACTGAAccctcc                                              >  1:170553/1‑139 (MQ=255)
                                                                                                  aTCTTTTATCTTCTGACGTTTTGGTATAAAAAAGCCCCTCGTTTGCGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACa                                             <  1:79645/143‑1 (MQ=255)
                                                                                                  aTCTTTTATCTCCTGACGCCTGGGAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAACAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAGGTCGAACACAGTGAACAATACa                                             >  1:596758/1‑143 (MQ=255)
                                                                                                  aTCTTTTATATCCTGACGCATGGGCAAAAAAAAGCCCCTCGATTGAGGTGCTATGCAGGGGCTAACAGCGGGCAGAAAAACGGCAGGGCAGCGTGTTTTCTTAGACGAGGTAAGACTAAATGTCGAACACACTGAACCATACa                                             >  2:221928/1‑143 (MQ=255)
                                                                                                     ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGCACACAGTGAACCGTACAtcc                                          >  1:140111/1‑143 (MQ=255)
                                                                                                         aTCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                      <  1:791581/143‑1 (MQ=255)
                                                                                                            tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGCTCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACACCCTCCCGGc                                   >  2:556203/1‑143 (MQ=255)
                                                                                                                   gCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCAGGCAAATTgc                            >  2:586216/1‑142 (MQ=255)
                                                                                                                   gCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTACCGGCAAATTTTc                           >  2:735325/1‑143 (MQ=255)
                                                                                                                     cTGGGCAAAAAAAAGCCCCTCGATTGGGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTcc                          <  1:829966/143‑1 (MQ=255)
                                                                                                                         gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAACTTGTCCGGCa                      <  1:604193/143‑1 (MQ=255)
                                                                                                                         gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCCTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATCGTCCGGCa                      <  1:257152/143‑1 (MQ=255)
                                                                                                                           aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCGAATTGTCCGGCAtt                    >  1:555282/1‑143 (MQ=255)
                                                                                                                            aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGGCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                   >  1:611557/1‑143 (MQ=255)
                                                                                                                               agaaGCCGCGCGAGTGAGGGGCTGTGAATGGGGGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATAc                <  2:610614/141‑1 (MQ=255)
                                                                                                                                ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGCTCAACGCGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCACATTGTCCGGCATTATACt               >  2:287094/3‑143 (MQ=255)
                                                                                                                                       cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACGTCCTCCCGGCAAATTGTCCGGCATGATACTCATCGTc        >  2:785753/1‑143 (MQ=255)
                                                                                                                                       cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCgtcgtc        >  1:739873/1‑143 (MQ=255)
                                                                                                                                             tGAGGGGCTGGGAGTGGGTGACCAACGGTAAGAAACACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAACTTGTCCGGCATTATACTCATCGTCAGAAgc  <  2:796591/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGC  >  NC_000913/4472296‑4472579

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: