Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,434 G→A 11.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,4340GA11.6% 111.8 / 8.6 54R34* (CGA→TGA) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base A (3/3);  total (31/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 = 4.73e-02

GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGA  >  NC_000913/4472298‑4472576
                                                                                                                                         |                                                                                                                                              
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:325721/1‑143 (MQ=255)
                                          tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                                                                >  1:46290/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)
                                                                                                                                         |                                                                                                                                              
GCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGA  >  NC_000913/4472298‑4472576

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: