Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,542 C→T 7.5% intergenic (‑155/+12) otsB ← / ← araH trehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,5420CT 7.5% 162.0 / 2.8 54intergenic (‑155/+12)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (28/26)
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.00e-01

CTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATT  >  NC_000913/1982397‑1982689
                                                                                                                                                 |                                                                                                                                                   
ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAg                                                                                                                                                   >  3:74710/1‑148 (MQ=255)
       tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAg                                                                                                                                            >  5:109328/1‑148 (MQ=255)
         ttttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACTGGGGTCGCCATATGCAACGAAGAACATCAAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCAt                                                                                                                                         >  2:17031/1‑149 (MQ=255)
           ttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCAt                                                                                                                                         <  8:4202/147‑1 (MQ=255)
                 ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGt                                                                                                                                 >  6:99558/1‑149 (MQ=255)
                   gtgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGc                                                                                                                               >  1:235481/1‑149 (MQ=255)
                   gtgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGc                                                                                                                               >  6:166143/1‑149 (MQ=255)
                     gtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGt                                                                                                                             >  4:218108/1‑149 (MQ=255)
                           cAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGAGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCt                                                                                                                       >  4:129035/1‑149 (MQ=255)
                            aGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGTAGAAAAAAGCATCAGACAGTGCGTTTCGCtt                                                                                                                      >  3:133156/1‑149 (MQ=255)
                               aCGCGCTTTACAGACATGTAAGTATAGACTGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                                                                   <  1:3505/149‑1 (MQ=255)
                               aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                                                                   <  1:199603/149‑1 (MQ=255)
                               aCGCGCTTTACAGACAGGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                                                                   <  1:158584/149‑1 (MQ=255)
                                 gcgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                                                                 >  8:34434/1‑149 (MQ=255)
                                  cgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCt                                                                                                                >  8:151870/1‑149 (MQ=255)
                                   gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                                                               <  1:168840/149‑1 (MQ=255)
                                   gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                                                               <  2:235481/149‑1 (MQ=255)
                                       tACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGt                                                                                                             <  6:218333/147‑1 (MQ=255)
                                       tACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCACTTTTTGCTTGt                                                                                                             <  1:76558/147‑1 (MQ=255)
                                         cAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc                                                                                                          <  3:129035/148‑1 (MQ=255)
                                            aCATGTAAGTATAGACAGTGTTACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGATTGTAACGGTc                                                                                                      >  6:45405/1‑149 (MQ=255)
                                            aCATGTAAGTATAGACAGTGTGACTCTGGTAGCAATTTGAAACTAATAAAATCAAACCGTTATGGCGGTTGGCTATGGTGGGAAAACACGCTAAATTGTTGCAGAAAAAAGCATCAGCGAGTGCGTTTCGCTTTTTGCTTGTAACGGTc                                                                                                      >  2:12913/1‑149 (MQ=255)
                                                gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCg                                                                                                                         <  7:139304/126‑1 (MQ=255)
                                                gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCg                                                                                                                         >  8:139304/1‑126 (MQ=255)
                                                    gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAg                                                                                                  <  5:170738/145‑1 (MQ=255)
                                                    gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAg                                                                                                  >  6:170738/1‑145 (MQ=255)
                                                       taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCg                                                                                           >  1:97187/1‑149 (MQ=255)
                                                             gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCca                                                                                     <  5:99558/149‑1 (MQ=255)
                                                                tgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTGTTAACGGTCGAAGATCACCGCTGCcag                                                                                    >  6:26139/1‑147 (MQ=255)
                                                                   cGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGTAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGa                                                                               <  4:133156/149‑1 (MQ=255)
                                                                            ccATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGAGTTTCGCTTATTGCTTGTAACGGTCGAAGATCACCCATGCCAGCAGTATTAAGCcgc                                                                      >  3:156771/1‑149 (MQ=255)
                                                                                   cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCa                                                                <  2:97187/148‑1 (MQ=255)
                                                                                    aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCTCTTTTTGATTGTATCGGTCTAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                                                              <  5:45405/149‑1 (MQ=255)
                                                                                        aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt                                                          <  4:74710/149‑1 (MQ=255)
                                                                                             cAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAACACGCTAAATTCTTGCAGAAAAACGCATCAGACAGTGCGTTTCGCTTTTTGGTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgagc                                                      >  2:72183/1‑148 (MQ=255)
                                                                                             cAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACCATAAATTGTTGAAGAAAAAAGCATCAGACAGTTCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACTTACTtcgc                                                      >  2:118392/1‑148 (MQ=255)
                                                                                               gCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                   <  3:218108/149‑1 (MQ=255)
                                                                                               gCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                   <  6:109328/149‑1 (MQ=255)
                                                                                                ccAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa                                                  <  7:193100/149‑1 (MQ=255)
                                                                                                  aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGATGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg                                                 >  2:218325/1‑148 (MQ=255)
                                                                                                       ttATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGCTCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaaa                                           >  2:117725/1‑149 (MQ=255)
                                                                                                          tGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAatat                                        >  3:94662/1‑149 (MQ=255)
                                                                                                            gCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTa                                      <  6:24317/149‑1 (MQ=255)
                                                                                                             cGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTaa                                     <  1:17031/149‑1 (MQ=255)
                                                                                                                ttGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCa                                  <  5:166143/149‑1 (MQ=255)
                                                                                                                             gAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGtt                     <  7:34434/149‑1 (MQ=255)
                                                                                                                                aaaaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGGTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAAGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTc                  <  2:244961/149‑1 (MQ=255)
                                                                                                                                       tAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa                                                  <  2:40844/110‑1 (MQ=255)
                                                                                                                                       tAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa                                                  >  1:40844/1‑110 (MQ=255)
                                                                                                                                       tAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTg           >  2:57536/1‑149 (MQ=255)
                                                                                                                                        aatttgttgTGGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCT‑‑TAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                                                               >  7:46082/4‑94 (MQ=255)
                                                                                                                                        aatttgttgTGGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCT‑‑TAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                                                               <  8:46082/91‑1 (MQ=255)
                                                                                                                                              ttgCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTaaaa    >  2:108153/1‑149 (MQ=255)
                                                                                                                                              ttgCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGAGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTaaaa    >  8:3933/1‑149 (MQ=255)
                                                                                                                                                 cAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAAtt  >  2:4066/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATT  >  NC_000913/1982397‑1982689

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