Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,122 A→G 9.2% 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.2% 135.6 / 5.1 56intergenic (+76/‑70)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (2/3);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACG  >  NC_000913/2089984‑2090256
                                                                                                                                          |                                                                                                                                      
acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCCTCCTGACTAGTATTTCAGGCGCTGTGTGCTGGAAGACATTCAGACCTTACCGTGCTGCATGAACGCCTGAGAGAGCGCCCAGACGATCAc                                                                                                                                    >  2:101495/1‑143 (MQ=255)
acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAGCGCATGAGAAAGCCCCCGGAAGATCAc                                                                                                                                    >  1:208177/1‑143 (MQ=255)
       cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCCCCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgg                                                                                                                             <  1:21704/143‑1 (MQ=255)
       cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCAGCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                             <  2:210477/143‑1 (MQ=255)
       cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATATTCTTGTCTAGTCGCTCAGGGCTTGAGCGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                             <  2:470449/143‑1 (MQ=255)
        ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGCAGACATTCAGATCTTCCAGTGGTGCATGTACGCATGTGAAAGCCCCCGGAAGATCTCCTTCCggg                                                                                                                            >  1:329489/1‑143 (MQ=255)
         tttATGACCCGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                           <  1:206203/143‑1 (MQ=255)
         tttATGACACGCGTTCAATTTAAACACCACCATCATCTCCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCGGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                           <  2:600820/143‑1 (MQ=255)
          ttATGACACGACATTAATTTAAACAAGCCCATCATCGCGCTCATCCTGGCTAGTCTTCCAGGCGATGTGTGCTGGAAGACATTCAGATGTTCCAGTGGCGCACGAACTCATGAGACAGCCCCCGGAAGATCACCTTCCggggg                                                                                                                          <  2:1039371/143‑1 (MQ=255)
                cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                                                                                                    <  1:908123/143‑1 (MQ=255)
                cacGCGTTCAATTCAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                                                                                                    <  1:298373/143‑1 (MQ=255)
                 acGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttttt                                                                                                                   >  1:813070/1‑143 (MQ=255)
                 acGCGTTCAATTCAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttttt                                                                                                                   <  1:560168/143‑1 (MQ=255)
                   gcgTTCAATTTATACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCTGAAGTTCACCTTCCGGGGGCTTTTTTAt                                                                                                                 <  2:90720/143‑1 (MQ=255)
                      ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgc                                                                                                              <  1:1036000/143‑1 (MQ=255)
                          aTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCgg                                                                                                          >  1:412055/1‑143 (MQ=255)
                           tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGCGGACCTTCAGATCTTCCAGCGGTGAATCAACGCAAAAGAAAGCCCCCGGACGATCAACATCCGCGGGCTCTTTCATTACGCGGt                                                                                                         >  2:16534/1‑143 (MQ=255)
                           tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGt                                                                                                         >  1:65074/1‑143 (MQ=255)
                                  accaccATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAAc                                                                                                  >  2:97510/1‑143 (MQ=255)
                                   ccgccatcatcagcatcatcCGGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTGCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACg                                                                                                 <  2:424414/140‑1 (MQ=255)
                                   ccaccaTCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACg                                                                                                 <  1:610558/143‑1 (MQ=255)
                                             caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGGGCACGAACGCATGAGAAAGCCCCCGGAAGATCACCTGCCGGGGGCTGTTGTATTGAGCGGTTGATAACGGCTCAGACAg                                                                                       >  2:855252/1‑143 (MQ=255)
                                             caccatcatcCTGACGAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                       <  1:595560/143‑1 (MQ=255)
                                               ccatcatcCTGACTAGTCTTTCGGGCGATGTGTGGTGGAAGACATTCAGGTCTTCCAGTGGTGCATGAACGCATGAGAAGGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGt                                                                                     <  2:57992/143‑1 (MQ=255)
                                                        tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAAGGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGGTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                            >  1:793881/1‑143 (MQ=255)
                                                          aCTAATCTTCCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGaa                                                                          <  2:854978/143‑1 (MQ=255)
                                                                   tCAGGGGATGTGGGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAt                                                                 <  2:652884/143‑1 (MQ=255)
                                                                     cggCGATGTGTGCGGGGAGACATTCAGACCCTCCAGTGGTGCAAGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga                                                               <  1:101495/142‑1 (MQ=255)
                                                                     aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga                                                               >  2:372616/1‑143 (MQ=255)
                                                                     aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAGga                                                               >  1:254255/1‑143 (MQ=255)
                                                                      ggCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgac                                                              <  1:811875/143‑1 (MQ=255)
                                                                       gtgaggtgcgCGGGAAGACATTCAGGTGTCCCAGTGGTGAGTGAACGGATGCGAAACCCCCCGGAAGAGCACCTTCCGGGGGGTTGTTTATTGCGCGGTTGATAACGGTGCAGACAGGTTTAAAGAGGAATAACAAACTgaca                                                             <  2:485481/138‑1 (MQ=255)
                                                                       gCGAGGTGCGCGGGAAGAGATTCAGAGCGTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgaca                                                             <  1:8549/143‑1 (MQ=255)
                                                                             tgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGACCGCATGAGAAAGCCCCCGGAAGATCACCTTGCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCACACAGGTTTAAAGAGGAATAACAAAATGACAGacaac                                                       >  2:233343/1‑143 (MQ=255)
                                                                                   ggAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGt                                                 <  1:83383/143‑1 (MQ=255)
                                                                                    gAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt                                                >  1:1034937/1‑143 (MQ=255)
                                                                                     aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAACAACAAAATGACAGACAACACTCGatt                                               >  2:906621/1‑143 (MQ=255)
                                                                                       tacacTCAGATCTTCGAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTAc                                             <  2:683474/142‑1 (MQ=255)
                                                                                        acgtgcAGATCGTCCGGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCATCCGGGGGTTGTTTTATTGTGCGGGTGATAACGGTGCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                            <  2:261924/138‑1 (MQ=255)
                                                                                          cttCAGATCTTGCAGTGGTGCTTGCACGTATGAGAAAGCCCCCGGAAGATCACCTTCCGGGCGCTTTTTTATTGCGCGGTTGATAACGGTTCTGACAGGCTTACAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                          <  2:324226/142‑1 (MQ=255)
                                                                                          aTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                          >  1:500689/1‑143 (MQ=255)
                                                                                           ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                         >  1:49009/1‑143 (MQ=255)
                                                                                                  cTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                  <  2:244691/143‑1 (MQ=255)
                                                                                                  cTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                  <  1:97510/143‑1 (MQ=255)
                                                                                                  cTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                  <  1:908927/143‑1 (MQ=255)
                                                                                                      cAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTACAGACAGGTTTAAAGAGGAATAACACCATGACAGACAACACTCGTTTACGCATAGCTATGCAGa                              >  2:857155/1‑143 (MQ=255)
                                                                                                      cAGTGGTGCATGAAAGCATGAGAAAGTTTCTGGAAGATCATCTTTTGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGa                              <  1:786771/143‑1 (MQ=255)
                                                                                                           gtgCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAACATGACAGACCACACTCGTTTACGCATCGCTATGCGGAAATcc                         >  1:4978/1‑143 (MQ=255)
                                                                                                                  aaCGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCCTAGCTATGCAGAACACGGGCCGtt                  >  2:156377/1‑143 (MQ=255)
                                                                                                                       aTGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGt             >  2:533587/1‑143 (MQ=255)
                                                                                                                                cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCa    >  1:846053/1‑143 (MQ=255)
                                                                                                                                cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGGTTAAGTGCTGCCTcc    >  2:822953/1‑142 (MQ=255)
                                                                                                                                  cgcGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg  <  2:458662/141‑1 (MQ=255)
                                                                                                                                  cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg  >  2:970777/1‑143 (MQ=255)
                                                                                                                                  cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGCTTACAGACGAATAACACAATGACAGACACCACTCGTTTACGCATATCTATGCAGAAATCAGGCCGTTTACGACATGAcgcccg  >  2:619528/1‑138 (MQ=255)
                                                                                                                                  cccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGCATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGGTTACGTGCTGACTCAcc  >  2:773465/1‑142 (MQ=255)
                                                                                                                                          |                                                                                                                                      
ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACG  >  NC_000913/2089984‑2090256

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: ‑