Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,316,188 G→A M4I (ATG→ATA rcsB → DNA‑binding transcriptional activator RcsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,316,1880GA97.6% 146.8 / ‑5.2 42M4I (ATG→ATArcsBDNA‑binding transcriptional activator RcsB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/19);  minor base C (0/1);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGAAACTGGATG  >  NC_000913/2316071‑2316324
                                                                                                                     |                                                                                                                                        
aCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGAcc                                                                                                                    <  2:502629/140‑1 (MQ=255)
 cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCa                                                                                                                   >  2:461463/1‑140 (MQ=255)
    tGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGAt                                                                                                                        >  1:337786/1‑132 (MQ=255)
    tGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGAt                                                                                                                        <  2:337786/132‑1 (MQ=255)
                    tGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGt                                                                                                <  2:720935/140‑1 (MQ=255)
                      agaAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTAt                                                                                              >  1:205142/1‑140 (MQ=255)
                       gaAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGAAAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATCAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTAtt                                                                                             <  2:681230/140‑1 (MQ=255)
                         aGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCg                                                                                           <  2:283517/140‑1 (MQ=255)
                            aTGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                                                        <  2:557263/140‑1 (MQ=255)
                            aTGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                                                        <  2:275346/140‑1 (MQ=255)
                                    ggAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTg                                                                                <  1:431380/140‑1 (MQ=255)
                                               aTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACt                                                                                  <  2:574285/127‑1 (MQ=255)
                                               aTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACt                                                                                  >  1:574285/1‑127 (MQ=255)
                                                tACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCATATTGAGt                                                                    >  2:490445/1‑140 (MQ=255)
                                                                      aTGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAtt                                              >  1:394882/1‑140 (MQ=255)
                                                                       tGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAttt                                             >  2:243131/1‑140 (MQ=255)
                                                                            agagCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa                                        <  2:668648/140‑1 (MQ=255)
                                                                            agagCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa                                        <  1:461463/140‑1 (MQ=255)
                                                                            agagCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa                                        <  2:205142/140‑1 (MQ=255)
                                                                            agagCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa                                        >  2:451118/1‑140 (MQ=255)
                                                                            agagCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa                                        <  2:394882/140‑1 (MQ=255)
                                                                               gCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGActc                                     <  1:243131/140‑1 (MQ=255)
                                                                                  tGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAc                                  >  2:274157/1‑140 (MQ=255)
                                                                                  tGCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAc                                  >  1:575745/1‑140 (MQ=255)
                                                                                   gCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACa                                 >  2:162363/1‑140 (MQ=255)
                                                                                   gCTGTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACa                                 >  1:501552/1‑140 (MQ=255)
                                                                                      gTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCATATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCa                              <  1:490445/140‑1 (MQ=255)
                                                                                      gTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCa                              <  1:274157/140‑1 (MQ=255)
                                                                                      gTAGCAAGGTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCa                              <  1:154098/140‑1 (MQ=255)
                                                                                              gTAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaa                      >  1:160815/1‑140 (MQ=255)
                                                                                               tAGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaac                     >  2:199472/1‑140 (MQ=255)
                                                                                                aGCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaaca                    <  2:187822/140‑1 (MQ=255)
                                                                                                 gCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaacaa                   >  2:342671/1‑140 (MQ=255)
                                                                                                 gCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaacaa                   >  2:522431/1‑140 (MQ=255)
                                                                                                 gCCTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATcaacaa                   >  2:373921/1‑140 (MQ=255)
                                                                                                  ccTATTACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCaacaac                  >  2:513921/1‑140 (MQ=255)
                                                                                                       tACATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGcc             >  1:205809/1‑140 (MQ=255)
                                                                                                        aCATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCg            >  2:516886/1‑140 (MQ=255)
                                                                                                        aCATGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCg            >  2:652361/1‑140 (MQ=255)
                                                                                                           tGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGaaa         <  1:652361/140‑1 (MQ=255)
                                                                                                           tGAACAATATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGaaa         >  2:367704/1‑140 (MQ=255)
                                                                                                                aaTATAAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGAAACTGGa    <  1:516886/140‑1 (MQ=255)
                                                                                                                  tataAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGAAACTGGATg  <  1:367704/140‑1 (MQ=255)
                                                                                                                     |                                                                                                                                        
ACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGATCAACAACCTGCCGAAACTGGATG  >  NC_000913/2316071‑2316324

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