Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,836,655 A→G 6.4% Y112H (TAC→CAC)  hisB ← imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,836,6550AG 6.4% 169.5 / 3.0 64Y112H (TAC→CAC) hisBimidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

TGTATTACCTCCACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGACAATGCCGATATC  >  CP000730/2836518‑2836790
                                                                                                                                         |                                                                                                                                       
tgtATTACCTCCACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATgg                                                                                                                                      <  2:543727/141‑1 (MQ=255)
 gtATTACCTCCACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATggg                                                                                                                                     >  1:9395/1‑141 (MQ=255)
        cTCCACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCAc                                                                                                                              <  1:142399/141‑1 (MQ=255)
           cACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTCTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTta                                                                                                                           <  1:18948/141‑1 (MQ=255)
            aCGAATTAAATCAATAGGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTTTGGGCGACAACTtat                                                                                                                          >  1:340460/1‑141 (MQ=255)
                    aaTCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACaa                                                                                                                  <  2:70470/141‑1 (MQ=255)
                       cAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAGCTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGa                                                                                                               <  1:516540/141‑1 (MQ=255)
                         atatGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGaca                                                                                                             <  2:224895/141‑1 (MQ=255)
                             gCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTg                                                                                                         >  1:448444/1‑141 (MQ=255)
                             gCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTg                                                                                                         >  1:239727/1‑141 (MQ=255)
                                    tAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGt                                                                                                  <  2:524160/141‑1 (MQ=255)
                                      aTCTTGCATTGATTACGCCCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGttt                                                                                                <  1:17607/141‑1 (MQ=255)
                                        cTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCAACTTATATCCACAACGACACGTGCTAATGTTTCa                                                                                              >  2:559527/1‑141 (MQ=255)
                                          tGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAGCTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCCAATggatgaaa                                                                                            <  2:41682/141‑8 (MQ=255)
                                                  ttACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCAACAACGACACGTGCTAATGTTTCATCCATTGGaa                                                                                    >  2:168313/1‑141 (MQ=255)
                                                   tACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGaag                                                                                   >  1:486512/1‑140 (MQ=255)
                                                   tACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAAt                                                                                   >  1:202784/1‑141 (MQ=255)
                                                     cGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACGACGACACGTGCTAATGTTTCATCAATTTGAATGt                                                                                 >  2:18081/1‑141 (MQ=255)
                                                      gACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTa                                                                                >  1:423057/1‑141 (MQ=255)
                                                       aCCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTGTGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTAc                                                                               >  1:379036/1‑141 (MQ=255)
                                                       aCCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTGTGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTAc                                                                               >  1:378872/1‑141 (MQ=255)
                                                       aCCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTAc                                                                               <  1:536103/141‑1 (MQ=255)
                                                          gCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATc                                                                            <  2:377831/141‑1 (MQ=255)
                                                               aaaaaaTTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTcc                                                                       <  2:257815/141‑1 (MQ=255)
                                                                  aaaTTCTTCTCCTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATa                                                                    <  1:239754/141‑1 (MQ=255)
                                                                   aaTTCTTCTACTTACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATaa                                                                   <  1:142373/141‑1 (MQ=255)
                                                                   aaTTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCCTTACTTCATGATGCATTGAATGATAGGAATGGGCGCCCACTTATTTGCCGAACTAGACGCGCTAATGTTTCAGCCATTGGCATGGAACTCGGTACATaa                                                                   >  1:270376/1‑141 (MQ=255)
                                                                        ttctACTAACTCCGGATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGCCACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATaacgaac                                                              <  1:187730/141‑1 (MQ=255)
                                                                          ctactaACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAacgaacga                                                            >  2:23975/1‑141 (MQ=255)
                                                                            actaACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGaaa                                                          >  1:34273/1‑141 (MQ=255)
                                                                                 cTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCtt                                                     >  2:358150/1‑141 (MQ=255)
                                                                                  tCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCttt                                                    >  1:549909/1‑141 (MQ=255)
                                                                                    cGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCtattt                                                  >  1:554692/1‑141 (MQ=255)
                                                                                       aTCAAACGTACCAACTTTTTCTTTACTTAATGATTCATTGAATGATACGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTAtt                                               >  2:155705/1‑140 (MQ=255)
                                                                                        tCAAACGTACCAACTTTTTCTTTTCTTAATGATGCATTGAATGATAGGTATTGGCGCCAACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCt                                              <  2:377375/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTTCTTAATGATTCATTGAATGATAGGTATGGGCGCCCACTTATATCTACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCTTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  1:164454/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTGTGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  2:378872/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTGTGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  2:379036/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  2:68556/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  1:546661/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  1:449453/141‑1 (MQ=255)
                                                                                            aCGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTaa                                          <  1:212194/141‑1 (MQ=255)
                                                                                             cGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCCCCCCTTCGAATGGACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAAt                                         <  1:251168/141‑1 (MQ=255)
                                                                                                 ccAACTTTTTCTTTGCTTAATGATGCATTGAATGATAGGTATGGCCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCAtt                                     >  1:312153/1‑141 (MQ=255)
                                                                                                 ccAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCAtt                                     >  1:184117/1‑141 (MQ=255)
                                                                                                   aaCTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTc                                   >  1:91585/1‑141 (MQ=255)
                                                                                                    aCTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCa                                  >  1:511393/1‑141 (MQ=255)
                                                                                                     cTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCAAATGATATCCACCACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGGAATGCTTTTTATCTTTAATCAtttaca                                 >  2:20265/1‑138 (MQ=255)
                                                                                                     cTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGACGCCCACTTATATACACAACAACACATGCTAATGTTTCATCCATTGGAATGTACATCGTTGCATAACAAACGAAATGCTTTTTATCTTTAATAATTTCaa                                 >  2:20261/1‑141 (MQ=255)
                                                                                                       tttcccttacttacAGATGCATTGAATGATAGGTATGGGCGCCCACCTATATCCCCAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGt                               <  2:352650/135‑1 (MQ=255)
                                                                                                        ttccTTTCCTTAATGATGCATCGCATGATAGGCATGGGCGCCCACTCATATCCACAACGACACGTGCTAAAGTTTCGTTCGTTATATTGTTTATCTTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTa                              <  2:261553/138‑1 (MQ=255)
                                                                                                           cTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACACCGTTCCATAACGAACGAAATGGTTTTTATCTTTAATCATTTCACGTAACa                           >  1:246133/1‑141 (MQ=255)
                                                                                                            tttACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACaa                          >  1:445594/1‑141 (MQ=255)
                                                                                                            tttACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACaa                          <  2:448444/141‑1 (MQ=255)
                                                                                                            ctgaCTTCATGGTGCATTGAATTATAGGTATTGTCGCCGCCTTATATCCACAACGACGCGTGCCAATGTTTCATCCTTTGGAATGTACTCCGTTCCATAACGAACGAAATGCTTTTTACCTTTAATCATTTCAAGTAACaa                          <  2:464539/138‑1 (MQ=255)
                                                                                                             ttACTTAATGATGCGTTGACTGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAATTAACAAt                         >  1:314645/1‑141 (MQ=255)
                                                                                                             ttACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAAt                         >  1:256069/1‑141 (MQ=255)
                                                                                                               aCTTAATGATGCATTGAATGATAGGTATGGGCTCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATTGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTg                       <  1:382638/141‑1 (MQ=255)
                                                                                                                      gatgCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATg                <  2:9395/141‑1 (MQ=255)
                                                                                                                       atgCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGa               <  2:428179/141‑1 (MQ=255)
                                                                                                                       atgCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGa               <  1:168313/141‑1 (MQ=255)
                                                                                                                        tgCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGACCAAACTGCTCTGTATCGTTAATCATTTCCAGTAACAATTGGCCAATGAc              >  1:105588/1‑141 (MQ=255)
                                                                                                                         gCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGACa             <  1:18081/141‑1 (MQ=255)
                                                                                                                                    tAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGACAATGCCGATATc  <  2:530474/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TGTATTACCTCCACGAATTAAATCAATATGCGTTGTTAATCTTGCATTGATTACGACCGCTCTAAAAAATTCTTCTACTAACTCCGTATCAAACGTACCAACTTTTTCTTTACTTAATGATGCATTGAATGATAGGTATGGGCGCCCACTTATATCCACAACGACACGTGCTAATGTTTCATCCATTGGAATGTACATCGTTCCATAACGAACGAAATGCTTTTTATCTTTAATCATTTCAAGTAACAATTGGCCAATGACAATGCCGATATC  >  CP000730/2836518‑2836790

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: