Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,122 A→G 9.7% intergenic (+76/‑70) hisL → / → hisG his operon leader peptide/ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,1220AG 9.7% 89.7 / 3.2 42intergenic (+76/‑70)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (14/22);  new base G (2/2);  total (18/24)
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.41e-01

ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAA  >  NC_000913/2089984‑2090260
                                                                                                                                          |                                                                                                                                          
acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCCGGCGATGTGTGCTGGGAGACATTCAGCTCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGGAGATCAc                                                                                                                                        >  2:46292/1‑143 (MQ=255)
  agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGCTGCATGAACGCATGAGTCAGCCTCCGGCAGCTCAccg                                                                                                                                      >  1:839179/1‑142 (MQ=255)
        ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                <  2:497502/143‑1 (MQ=255)
         tttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                               <  2:42952/143‑1 (MQ=255)
          ttATGACACGCGTTCTATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                                                                                                              <  2:139332/143‑1 (MQ=255)
               acacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttt                                                                                                                         >  1:545194/1‑143 (MQ=255)
                    cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                                                                                                    <  2:367921/143‑1 (MQ=255)
                      ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGACCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgc                                                                                                                  <  1:579628/143‑1 (MQ=255)
                       tCAATCTAAACACCCCCATCATCACCATCAGCCTGACTAGTCTGTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                                                                                                 <  2:764075/143‑1 (MQ=255)
                          aTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCgg                                                                                                              <  1:574749/143‑1 (MQ=255)
                           tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGACGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAACCCCCGAAAGATCACCTTCCGGGGGCTTTATTATTAcgcgat                                                                                                             >  1:744253/1‑141 (MQ=255)
                                  accatcATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAAc                                                                                                      <  1:53302/143‑1 (MQ=255)
                                             caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTTGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                           <  2:418510/143‑1 (MQ=255)
                                                catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAGCATATAAATAATCGTTTGGGGGCAACCCAACCGGGGGCTTTTTTATTGCGCGGTCGATAACGGTTCAGACAGGtt                                                                                        >  2:188241/1‑143 (MQ=255)
                                                catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGCCAggat                                                                                        >  2:19051/1‑141 (MQ=255)
                                                catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGtt                                                                                        >  1:322828/1‑143 (MQ=255)
                                                  gcagccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATGCAGAGCTGCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTa                                                                                      <  1:556057/139‑1 (MQ=255)
                                                      ccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                  <  2:147422/143‑1 (MQ=255)
                                                        tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGGAAGCCCCCGGAAGATCACCTTCCGGGGGCTGTTTTATTGCGCGGTTGATGACGGTTCTGTGAGGttttcagggg                                                                                >  2:667109/1‑136 (MQ=255)
                                                        tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                                >  1:429801/1‑143 (MQ=255)
                                                             aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATaa                                                                           <  2:465090/143‑1 (MQ=255)
                                                              gTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATGGCGCGGGTGATCACGGTTCAGACGGGTTTGAGGAGGAATagc                                                                          >  1:779094/1‑141 (MQ=255)
                                                               tCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACa                                                                         <  1:147570/143‑1 (MQ=255)
                                                                 tttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTATAAAGAGGAaaaaaaaa                                                                       >  1:213687/1‑143 (MQ=255)
                                                                 ttgcAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaa                                                                       <  2:657435/140‑1 (MQ=255)
                                                                  ttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                      <  1:771167/143‑1 (MQ=255)
                                                                     aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga                                                                   >  1:18294/1‑143 (MQ=255)
                                                                                tcTGGAAGACATTCAGATCTTCCCGTGGTGCAACAGCCTGCGTTTTGTCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACt                                                        <  2:867460/142‑1 (MQ=255)
                                                                                    gAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGAGTCTCTTGTATTGCGCGGTTGATAACGGTTCAGACAGGTTCAACGAGCAATAACAAAATGACAGACACAACTCgct                                                    >  2:47451/1‑141 (MQ=255)
                                                                                        aCATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCCTTTTCATTGCGCGGTTCATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                                <  1:227842/143‑1 (MQ=255)
                                                                                           ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                             >  2:566434/1‑143 (MQ=255)
                                                                                            tCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATa                                            <  1:770581/143‑1 (MQ=255)
                                                                                                   ttCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGccacac                                     >  1:196533/1‑138 (MQ=255)
                                                                                                           gtgCTTGAATGCATGAGAAAGCACCCGGAAGATCAACTTCGGGGGGCGTTGTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAGTGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATcc                             <  1:565634/143‑1 (MQ=255)
                                                                                                            tgCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCg                            <  2:70611/143‑1 (MQ=255)
                                                                                                             gcttGAACGCATGAGAAAGCGCCCGGAAGATCACCTTCCGGGGGCTCTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAACGAGGAATAACAAAATGACAGTCAACCCTCGTTTACGCATAGCTACGCAGAAATCCgg                           <  2:616844/140‑1 (MQ=255)
                                                                                                                                cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCa        <  1:240806/143‑1 (MQ=255)
                                                                                                                                 ccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAGAGCTATGCAGAAATCCGGCCGATTAAGTGATGACTCAc       >  1:470384/1‑143 (MQ=255)
                                                                                                                                  cccTGAAGATCACCTCCCGGGGGCGTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAAGAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg      <  1:46292/143‑1 (MQ=255)
                                                                                                                                  cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAACGAGGAAGAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGCCCGCGTAAGTGATGACTCAcg      >  2:2019/1‑143 (MQ=255)
                                                                                                                                  cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGACTAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATACCGCCGTTTAAGTGATGACTCAca      >  1:872791/1‑142 (MQ=255)
                                                                                                                                      gAAGATCACCTTCCGGGGGCGTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGaa  <  1:681250/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAA  >  NC_000913/2089984‑2090260

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: