Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,096,668 Δ1 bp 7.1% coding (55/612 nt) hisI → phosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,6680A. 7.1% 183.2 / 6.8 56coding (55/612 nt)hisIphosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (18/34);  new base . (2/2);  total (20/36)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAA  >  NC_000913/2096520‑2096805
                                                                                                                                                    |                                                                                                                                         
aCGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACACCAACGGCGCGAACTGGACTGGGAAAAAACCGCCGGCCTTa                                                                                                                                           >  2:218399/1‑149 (MQ=255)
      cgcTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCgg                                                                                                                                     <  7:20127/149‑1 (MQ=255)
           gCAGCTTCCGTCTTCCACAAACAAATAATCAATATTGGTGAATTAAACGCGTCCCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGTCGGTGa                                                                                                                                  <  6:101797/147‑1 (MQ=255)
           gCAGCTTCCGTATTCCACAAACAAATAATCAATCTTGGTGAATTAAAAGCGGACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGa                                                                                                                                  <  4:234790/147‑1 (MQ=255)
           gCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                                                >  5:381325/1‑149 (MQ=255)
           gCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                                                <  6:38742/149‑1 (MQ=255)
           gCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGACAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                                                <  6:338772/149‑1 (MQ=255)
             aGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATtgt                                                                                                                              <  5:361554/149‑1 (MQ=255)
             aGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGGACCTCGCAACACCGGGCGGGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATtgt                                                                                                                              <  6:343018/149‑1 (MQ=255)
                  ccGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCCTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCa                                                                                                                           <  2:300607/147‑1 (MQ=255)
                       ttCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACAcgcg                                                                                                                    <  2:113008/149‑1 (MQ=255)
                          cacaAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGACCAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTa                                                                                                                 <  2:133764/149‑1 (MQ=255)
                          cacaAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGt                                                                                                                  <  1:11305/148‑1 (MQ=255)
                          cacaAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGGGCAACACGCGGTa                                                                                                                 <  1:312087/149‑1 (MQ=255)
                                      atCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTg                                                                                                     <  3:191546/149‑1 (MQ=255)
                                           tatTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTTGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGGGCAACACGCGGTATCCGGCGAAGTGCTAAt                                                                                                <  2:47033/149‑1 (MQ=255)
                                                  gAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGc                                                                                         <  2:236988/149‑1 (MQ=255)
                                                   aaTTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCt                                                                                        <  4:215829/149‑1 (MQ=255)
                                                     ttAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCta                                                                                       >  1:258021/1‑148 (MQ=255)
                                                       aaaaGAGTACCTGGCAACACAGGGCGTGGATATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAATCCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCtata                                                                                     <  5:154397/148‑1 (MQ=255)
                                                        aaaGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                                                >  2:166562/1‑104 (MQ=255)
                                                        aaaGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                                                <  1:166562/104‑1 (MQ=255)
                                                              tACCTGGCAACACAGGGCGTGGAGATCATGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCg                                                                             <  1:253563/149‑1 (MQ=255)
                                                              tACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCt                                                                                             >  5:193492/1‑133 (MQ=255)
                                                              tACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCt                                                                                             <  6:193492/133‑1 (MQ=255)
                                                                 cTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGaa                                                                          <  5:173804/149‑1 (MQ=255)
                                                                      aaCACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAAc                                                                                >  8:18533/1‑138 (MQ=255)
                                                                      aaCACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAAc                                                                                <  7:18533/138‑1 (MQ=255)
                                                                      aaCACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCtt                                                                     <  7:87374/149‑1 (MQ=255)
                                                                       acacAGGGCGTGGAGATCAGGATATGTTAACAGAACATCAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTa                                                                    <  1:377999/149‑1 (MQ=255)
                                                                         ccaGGGCGGGGACATCAGGAGATGATAACAGAACACCCACGTCGCCAACTGGACTGGGAAAAAACCGACCGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGa                                                                  <  6:234190/148‑1 (MQ=255)
                                                                                          ggATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCg                                                 >  5:299807/1‑149 (MQ=255)
                                                                                             tatGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCgg                                                                                                                   >  7:42174/1‑80 (MQ=255)
                                                                                             tatGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCgg                                                                                                                   <  8:42174/80‑1 (MQ=255)
                                                                                             tatGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATAGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGc                                               >  1:163497/1‑148 (MQ=255)
                                                                                               tGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCa                                              >  4:200049/1‑147 (MQ=255)
                                                                                                gTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGTTTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCaaa                                            <  7:38059/148‑1 (MQ=255)
                                                                                                gTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAg                                           >  5:103847/1‑149 (MQ=255)
                                                                                                   aaCAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTc                                         <  6:381325/148‑1 (MQ=255)
                                                                                                       gAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAAccc                                                          <  8:22583/127‑1 (MQ=255)
                                                                                                       gAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAAccc                                                          >  7:22583/1‑127 (MQ=255)
                                                                                                        aacaacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAACTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACcttc                                   >  3:6407/1‑149 (MQ=255)
                                                                                                           aacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACaaaa                                                             <  7:259054/120‑1 (MQ=255)
                                                                                                           aacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACaaaa                                                             >  8:259054/1‑120 (MQ=255)
                                                                                                            acaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTtcttct                               <  6:299807/149‑1 (MQ=255)
                                                                                                              aacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTcg                             >  2:192608/1‑149 (MQ=255)
                                                                                                                      gAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTaa                     <  5:289098/149‑1 (MQ=255)
                                                                                                                         ctggactggGAAAAAACAGAAGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACa                  <  2:258021/149‑1 (MQ=255)
                                                                                                                          tggactggGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGCCAAAACCCTCGAAAGCGGCAAAGTCCCCGTCTTCTCGCGCACTAAACAg                 >  4:71031/1‑149 (MQ=255)
                                                                                                                            gactggGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGc                >  3:105308/1‑148 (MQ=255)
                                                                                                                             actggGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGa              >  8:66259/1‑149 (MQ=255)
                                                                                                                             actggGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGa              >  4:315674/1‑149 (MQ=255)
                                                                                                                               tggGAAAAAACCGACGGACT‑‑TGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACtgt          >  6:240619/1‑149 (MQ=255)
                                                                                                                               tggGAAAAAACCGACGGACT‑‑TGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACtgt          >  6:235201/1‑149 (MQ=255)
                                                                                                                                   aataaaCCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGACCCCGGAAGCCCTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCAc                        <  6:173359/130‑1 (MQ=255)
                                                                                                                                       aaCCGACGGACT‑‑TGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCaaa  <  5:240619/149‑1 (MQ=255)
                                                                                                                                       aaCCGACGGACT‑‑TGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCaaa  <  5:235201/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
ACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAA  >  NC_000913/2096520‑2096805

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