Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,776 A→G 7.3% N91D (AAT→GAT)  yaiV → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,7760AG 7.3% 159.7 / 3.2 55N91D (AAT→GAT) yaiVputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/22);  new base G (2/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

ACTTTCGTTATTTTAGAGGGAGTTATCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCATTGGGCATAATTCCTACG  >  NC_000913/394635‑394923
                                                                                                                                             |                                                                                                                                                   
aCTTTCGTTATTTTAGAGGGAGTTATCTCTATGCGTAGAGAAGAAAACGTACTTATCGTTATTACCAAGGATCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACAATACAAATTAATATAAGAAGGAAATTGTAc                                                                                                                                              >  1:406850/1‑149 (MQ=255)
  tttCGTTATTTTAGAGGGAGTTATCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACgg                                                                                                                                            >  1:34198/1‑149 (MQ=255)
      gTTATTTTAGAGGGAGTTATCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGata                                                                                                                                        >  3:384474/1‑149 (MQ=255)
              agagGGAGTTATCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTAc                                                                                                                                >  7:359771/1‑149 (MQ=255)
                  tgAGTTATCTCTATGCGTAGAGAAGAAAACGTACTTATCTGTATTACCCAGGCTCCTTATATTATGGGGATGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGc                                                                                                                            <  2:406850/148‑1 (MQ=255)
                        aTCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCCGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCcaaac                                                                                                                       >  1:76959/1‑148 (MQ=255)
                          ctctctGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCcaaacaaa                                                                                                                    >  4:40226/1‑149 (MQ=255)
                            ctctGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCaaacaaac                                                                                                                   >  2:408659/1‑148 (MQ=255)
                               tGCGTAGAGAAGAAAACGTACTTATCGGTATTACACAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCAt                                                                                                                >  4:58906/1‑148 (MQ=255)
                               tGCGTAGAGAAGAAAACGTACGTATCGGTATTACCCAGGCTCCTTATATTATGGGGGTGGCGGATGGTTTTATGAAAAAAGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCAt                                                                                                                >  2:64649/1‑148 (MQ=255)
                                 cGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTAc                                                                                                             <  2:34198/149‑1 (MQ=255)
                                  gTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACg                                                                                                            >  6:395351/1‑149 (MQ=255)
                                          gaaAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTg                                                                                                     >  4:128687/1‑148 (MQ=255)
                                          gaaAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTg                                                                                                     >  2:203478/1‑148 (MQ=255)
                                          gaaAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTg                                                                                                     >  7:53301/1‑148 (MQ=255)
                                                    ttATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCaa                                                                                          >  4:80408/1‑149 (MQ=255)
                                                      aTCGGTATTACCCAGTCTCCTTATATTATGGGGCTGGCTGATTGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATAATCTACCAGCCAAACAAACCATTAAGCTTATTGAACAAAATaaaat                                                                                        >  7:393617/1‑147 (MQ=255)
                                                      aTCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAAc                                                                                         <  1:408659/148‑1 (MQ=255)
                                                         ggTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTg                                                                                     <  3:58906/149‑1 (MQ=255)
                                                          gTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTgg                                                                                    >  2:141062/1‑149 (MQ=255)
                                                             ttACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCga                                                                                 <  3:40226/149‑1 (MQ=255)
                                                                cccAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCgaga                                                                               <  8:53301/148‑1 (MQ=255)
                                                                   aGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCt                                                                           >  3:145670/1‑149 (MQ=255)
                                                                         cTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTAc                                                                     <  1:141062/149‑1 (MQ=255)
                                                                           tataTTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAGATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTg                                                                   >  2:91791/1‑149 (MQ=255)
                                                                           tataTTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAGATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTg                                                                   >  2:86013/1‑149 (MQ=255)
                                                                                     ggctggctGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTgg                                                         >  1:26663/1‑149 (MQ=255)
                                                                                           ctgATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATTTCATCTACCAGCCAAACAAACCAGTACGCTTATTGAACAAAATCACCTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAAt                                                   <  5:349542/149‑1 (MQ=255)
                                                                                            tgATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCaaaa                                                    >  4:304861/1‑147 (MQ=255)
                                                                                             gATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACGAGCCAAACAAACCATTACGCGTTTTGAAGAAAATCAGGTCTGGAGAGGGGGTTTTTACTGGTTAGCGGGGGATAATAg                                                 >  2:313704/1‑149 (MQ=255)
                                                                                              aTGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATa                                                  <  3:304861/147‑1 (MQ=255)
                                                                                                 gTTTAATGAAAAACGAAATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAAtt                                             <  6:312051/149‑1 (MQ=255)
                                                                                                  tttAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGt                                                                 <  2:90593/128‑1 (MQ=255)
                                                                                                  tttAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGt                                                                 >  1:90593/1‑128 (MQ=255)
                                                                                                   ttAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCt                                           >  2:155694/1‑149 (MQ=255)
                                                                                                    tAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAGATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCt                                           <  1:86013/148‑1 (MQ=255)
                                                                                                    tAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAGATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCt                                           <  1:91791/148‑1 (MQ=255)
                                                                                                     aaTGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTg                                          <  1:203478/148‑1 (MQ=255)
                                                                                                      aTGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGa                                        <  4:400089/149‑1 (MQ=255)
                                                                                                       tGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGaa                                       <  3:128687/149‑1 (MQ=255)
                                                                                                        gAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTg                                                                   >  5:206231/1‑120 (MQ=255)
                                                                                                        gAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTg                                                                   <  6:206231/120‑1 (MQ=255)
                                                                                                             aCGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCAACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTAcgc                                 <  1:313704/149‑1 (MQ=255)
                                                                                                                   taCCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGt                            <  4:145670/148‑1 (MQ=255)
                                                                                                                   taCCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTg                           <  1:64649/149‑1 (MQ=255)
                                                                                                                    aCCATACAAATTAATATCAGAAGGAAATTGTACGCGATATCATCTACCACCCAAACAAAACATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTGAATGGTTAGCATGGCAAAATAGAATTCTGGAATTACGCGACGtc                           >  1:227558/1‑147 (MQ=255)
                                                                                                                     ccATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTg                           >  2:276292/1‑147 (MQ=255)
                                                                                                                       aTACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGc                       >  2:28549/1‑149 (MQ=255)
                                                                                                                        tACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCttt                                                                         >  6:84459/1‑98 (MQ=255)
                                                                                                                        tACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCttt                                                                         <  5:84459/98‑1 (MQ=255)
                                                                                                                            aaTTAATATCAGAAGGTAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACTAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCAtt                  >  5:195065/1‑149 (MQ=255)
                                                                                                                              ttAATATCAGAAGGTAATTGTACGGGCTATCATCTTCCAGCCAAACAAACCATTACTCTTATTGAACAAAATCAACTCTTGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCATTgg                <  8:393617/149‑1 (MQ=255)
                                                                                                                                          ggAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCATTGGGCATAATTCCTa    >  8:35810/1‑149 (MQ=255)
                                                                                                                                           gAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGGGCAGCTCATTGGGCATAATTCCTAc   <  2:26663/149‑1 (MQ=255)
                                                                                                                                           gAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCGTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCATTGGGCATAATTCCTAc   <  1:417934/149‑1 (MQ=255)
                                                                                                                                             aaTTGTACGGGATATCATCAACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGTCGCGACGATTTTTACTGGTTAGCCAGGCAAAATAGAAGTATGGAATAACGCGACGTGCAGCTCATTGGGCATAATTCCTACg  >  2:387282/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ACTTTCGTTATTTTAGAGGGAGTTATCTCTCTGCGTAGAGAAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGAACAAAATCAACTCTGGCGAGACGCTTTTTACTGGTTAGCCTGGCAAAATAGAATTCTGGAATTACGCGACGTGCAGCTCATTGGGCATAATTCCTACG  >  NC_000913/394635‑394923

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