Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,686 C→G intergenic (‑136/+149) araD ← / ← araA L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,6860CG96.9% 86.1 / NA 32intergenic (‑136/+149)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (12/19);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATT  >  NC_000913/66449‑66922
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
gcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCt                                                                                                                                                                                                                                 <  2:27589/251‑1 (MQ=255)
               gACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGcc                                                                                                                                                                                                                    >  1:176297/1‑249 (MQ=255)
                    ggCCCCATGTGAGCGTAACCAGGTTATGTTTTGGGTGCGCCAGGTTGGCCTCGAATACCTGGCGTTTGAGAGCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGACTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAGGACGCGCCAGCGTc                                                                                                                                                                                                              <  2:275383/249‑1 (MQ=21)
                                gCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTg                                                                                                                                                                                                 <  1:240136/251‑1 (MQ=255)
                                          gttatgtTTTGGCAGCGCCAGGTTGGCCTCTAATCCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAt                                                                                                                                                                                       <  1:79595/247‑1 (MQ=255)
                                               gtTTTGGCAGCGCCAGGTTGCCCTCTAATGCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGGGGCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCt                                                                                                                                                                                  <  2:62055/251‑1 (MQ=21)
                                                  ttGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCttt                                                                                                                                                                                <  1:317813/250‑1 (MQ=255)
                                                       aGCGCCAGGTTGGCCTCTCATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGCTAAGAAGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTAt                                                                                                                                                                          >  2:128417/1‑251 (MQ=21)
                                                              ggTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCt                                                                                                                                                                   <  1:306265/251‑1 (MQ=255)
                                                                    cctCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACaaa                                                                                                                                                              >  2:261765/3‑250 (MQ=255)
                                                                                ggCGTTTGAGATCTGCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTa                                                                                                                                                 <  1:116413/251‑1 (MQ=255)
                                                                                                   cATGTTGACTCCTTCGTACCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCATCGTGTCGGCCTGATAAGCCGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGCGTAGGCCTGATAAGCCTCGCTCGCGTCGCATCAGCCGTTGCATGCCGGCTTCGCTTTGCTTATCCGTCATACAAATTCCCTGCGTGTAGGccttatcagcagcgtcc                                                                                                                              >  2:308732/1‑236 (MQ=12)
                                                                                                      gTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCg                                                                                                                           >  1:286185/1‑251 (MQ=255)
                                                                                                           cTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGc                                                                                                                       >  1:2163/1‑250 (MQ=255)
                                                                                                             ccTTCGTGGCGGATGCGCTGTGATTAGCGGGCTAACAAAATCGCAGCGAGAAGGCCTGCGCAGACTCTCCAGCGTCGAATCAGGCGGTGAATGGCGGATGCGCTTTACGTATCCGGCCTTCAAGATCGGAACGGGTAGGCTTGCTAATACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGAGAAGACGCTGGAGCGTCGCATc                                                                                                                    <  2:313906/251‑1 (MQ=9)
                                                                                                               ttCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATc                                                                                                                    >  2:43541/1‑249 (MQ=255)
                                                                                                                 cGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGc                                                                                                                <  1:266164/251‑1 (MQ=255)
                                                                                                                      cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTAc                                                                                                           >  2:242754/1‑251 (MQ=255)
                                                                                                                        gATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                         <  2:28231/251‑2 (MQ=255)
                                                                                                                                          ggCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACtt                                                                                        >  2:241561/1‑250 (MQ=255)
                                                                                                                                           gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAg                                                                                      <  1:248199/251‑1 (MQ=255)
                                                                                                                                                                      gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGt                                                             <  1:99404/249‑1 (MQ=18)
                                                                                                                                                                       aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGGTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTcc                                                          >  1:240240/1‑251 (MQ=18)
                                                                                                                                                                         aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcag                                                        <  1:138992/251‑1 (MQ=18)
                                                                                                                                                                               cgcCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATGCACTTCGTTCcagcgcagc                                                  <  2:181840/251‑1 (MQ=18)
                                                                                                                                                                                             cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGGCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGCAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                    >  2:245600/1‑251 (MQ=18)
                                                                                                                                                                                                    tGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAACCGCTGGccgg                             >  2:138804/1‑248 (MQ=18)
                                                                                                                                                                                                    tGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCat                              <  2:170008/250‑2 (MQ=18)
                                                                                                                                                                                                         gCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgtcgagt                        <  2:51700/251‑7 (MQ=18)
                                                                                                                                                                                                                          tATCCGGCATACAAAATCGGCACGCGTAGGCCTGCTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATCCGTACTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGCTAACAATCAcc        >  2:116385/1‑250 (MQ=18)
                                                                                                                                                                                                                           aTCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGt      <  1:241561/251‑1 (MQ=21)
                                                                                                                                                                                                                               ggCCTACAAAATCGGAGCGCGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCTACTACTTAGCGACGAAATCCGTAATACACTTCGGTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGTGGGTATCGTTATCAATCACCTTGAtt  <  2:282992/251‑1 (MQ=18)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
GCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATT  >  NC_000913/66449‑66922

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