Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 12.9% 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,4370TC12.9% 108.0 / 10.0 55N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base C (4/3);  total (32/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 = 9.69e-02

GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAG  >  NC_000913/4472298‑4472578
                                                                                                                                            |                                                                                                                                             
gCTTTCAGCTGCGCCGCTGTCGCCAGCACCAGATGAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAg                                                                                                                                            <  1:115288/143‑1 (MQ=255)
gCTTGCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAg                                                                                                                                            <  1:522740/143‑1 (MQ=255)
   ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTGATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGGTACCTTTTATCTACTGACGCGTGAGCAAGAAAAAGCGCCTCGATTGAgggg                                                                                                                                         >  2:992569/1‑143 (MQ=255)
     cAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                       >  1:370566/1‑143 (MQ=255)
       gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttt                                                                                                                                     <  1:849519/143‑3 (MQ=255)
            gccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttttttt                                                                                                                                <  1:788386/143‑8 (MQ=255)
                   tCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                                                                         >  1:987596/1‑143 (MQ=255)
                     gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGACTGAGGGGCTGGGAATGGGGGATCaa                                                                                                                       >  1:831225/1‑143 (MQ=255)
                      ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAAc                                                                                                                      >  1:524729/1‑143 (MQ=255)
                      ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTATTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCaca                                                                                                                      >  2:435620/1‑141 (MQ=255)
                       cAGCACCAGATTAAGGTCATCGCGAGTAAGGTCGTTTATGGCAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                     <  2:693305/143‑1 (MQ=255)
                       aagtaCCAGAGTAAAGTCAGCGAGATCAAGGTGTTTAATCGAAATGATATGTTGGTAATATAGCGGATGAGCGCTTTTTTATCTCCAGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                     <  2:280747/139‑1 (MQ=255)
                        aGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                                                                                    <  1:1009401/143‑1 (MQ=255)
                           aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAAGCGCCTCGCTTGAGGGGCGGAGGATGGGCGGTCACGGGga                                                                                                                  >  2:618848/1‑143 (MQ=255)
                           aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                 >  1:218042/1‑143 (MQ=255)
                                 gtAAGGTCATCGCTACTAAGGTCGTTTATGGAAATGATATGTTGCTGATATAGCGGATTAGCCATCTTTTATCTCCGGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                           <  2:1086336/142‑1 (MQ=255)
                                    aGGTCATCGCGACTAAGGTCGTTTATGGAACTGATATCTTTCTGATATAGCGGATTAGCCATCTTTTATCTCATGACCCCTTTGCAAAAAATGTCCCGTCGATGGAGTGGCCTGGCATGGGTGATCGACGGGCACAGAAACgg                                                                                                        >  2:543009/1‑143 (MQ=255)
                                       tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGCTTGAGGGGCTGGGAATGGGAGATCAACGGGAAGAAAAACGGCAg                                                                                                     >  2:320352/1‑143 (MQ=255)
                                       tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAGCGGGAAGAAAGGCGGCAg                                                                                                     >  1:77692/1‑143 (MQ=255)
                                       tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGcgg                                                                                                     >  1:613434/1‑141 (MQ=255)
                                         accgcgGCTAAGTGCGTTTATGTAAATGATATGTGTCTGAGATAGCGGATTAGCCATGTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                                                                   <  2:936503/141‑1 (MQ=255)
                                          tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                                                                  >  1:46290/1‑143 (MQ=255)
                                          tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                                                                  >  1:325721/1‑143 (MQ=255)
                                             cgACTAAGGTCGTTTATGGAAATGATATGTGTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                                                                               <  1:345133/143‑1 (MQ=255)
                                                 tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCGAGCGtgt                                                                                           >  2:156664/1‑141 (MQ=255)
                                                  aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                                                                          >  1:221763/1‑143 (MQ=255)
                                                       cGTTTCTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTAGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                     <  1:473174/143‑1 (MQ=255)
                                                       cGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                     <  2:370566/143‑1 (MQ=255)
                                                           tATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                                                                 >  2:190037/1‑143 (MQ=255)
                                                                   tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCCCGATTGCGGGGCTCGGCATGGGTGATCAACGGGACGAAAAACGGCAGGCCAGCGTCTTCTTTCAGCCGCGGTAAg                                                                         >  2:324042/1‑143 (MQ=255)
                                                                       atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaa                                                                     >  1:822016/1‑143 (MQ=255)
                                                                          tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                                                                  >  2:562997/1‑143 (MQ=255)
                                                                            tccgATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                <  2:427404/140‑1 (MQ=255)
                                                                                       ggATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGGCGAACCCACta                                                     >  1:323502/1‑142 (MQ=255)
                                                                                          ttacgcATCTTTTCTTCACTGTCGCCCGGCCAGAAAAGAGCCCCTCGAGTGAGGGCCTGGGAATGGGTGATCAACGAGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                  <  2:160043/138‑1 (MQ=255)
                                                                                              ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                              >  1:500254/1‑143 (MQ=255)
                                                                                                aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACa                                            >  2:423262/1‑143 (MQ=255)
                                                                                                       atcagcTGACCCATAGGCAAAAAAAAGCCCCCCGAGTGAGGGGCAGGGAATGGGTGATCAACGGGAAGAAAAACGGCCGTCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                     <  1:156664/138‑1 (MQ=255)
                                                                                                         ctcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                   >  1:740525/1‑143 (MQ=255)
                                                                                                            cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTACCCGcca                                >  1:796602/1‑141 (MQ=255)
                                                                                                             tGACGCCTGGGCAAAAAAACGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                               <  1:608669/143‑1 (MQ=255)
                                                                                                             tGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                               <  2:987596/143‑1 (MQ=255)
                                                                                                                       gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGcc                     >  1:292479/1‑142 (MQ=255)
                                                                                                                          aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCGAATTGTCCGGCAtt                   >  1:133875/1‑143 (MQ=255)
                                                                                                                           aaaaaaaGCCCCTCGATTTAGGCGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTGTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                  <  1:708883/143‑1 (MQ=255)
                                                                                                                             aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATtata                <  1:648632/143‑1 (MQ=255)
                                                                                                                                caGCCCCTCCATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATAAATCATCCCGGCAAATTGTCCGGCATTATACTc             >  2:144569/2‑143 (MQ=255)
                                                                                                                                caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc             >  2:133697/2‑143 (MQ=255)
                                                                                                                                caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc             >  1:385831/2‑143 (MQ=255)
                                                                                                                                caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc             >  2:831513/2‑143 (MQ=255)
                                                                                                                                  gCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCAt           <  1:656086/143‑1 (MQ=255)
                                                                                                                                     cgtCGATTGAGGGGCTGTGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt        <  1:722778/141‑1 (MQ=255)
                                                                                                                                     ccTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt        <  1:133697/143‑1 (MQ=255)
                                                                                                                                       tCGTTTGAGTGGCCGGGAACGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTATTTTCAGACGCGGTAAGACAAGATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa      <  1:1101100/143‑1 (MQ=255)
                                                                                                                                         gATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCTGCATTATACTCATCGTCAGa    >  1:75236/1‑143 (MQ=255)
                                                                                                                                           ttGAGGGGCTGGGAATGGGTGATCAACGGGAATAAAAACGGCAGGCCAGCGCCTTTCTTTATTCGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg  >  2:235426/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAG  >  NC_000913/4472298‑4472578

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