Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,877,011 G→A intergenic (‑59/+56) yidA ← / ← yidB sugar phosphate phosphatase/DUF937 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,0110GA78.6% 52.6 / 9.5 28intergenic (‑59/+56)yidA/yidBsugar phosphate phosphatase/DUF937 family protein
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (12/10);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGC  >  NC_000913/3876870‑3877157
                                                                                                                                             |                                                                                                                                                  
gCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaaca                                                                                                                                             >  5:139220/1‑149 (MQ=255)
      ttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGa                                                                                                                                        <  5:93180/148‑1 (MQ=255)
           ggcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGaa                                                                                                                                   >  6:25748/1‑148 (MQ=255)
           ggcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGaa                                                                                                                                   >  5:150768/1‑148 (MQ=255)
                 tGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggg                                                                                                                             >  6:115797/1‑148 (MQ=255)
                         tgtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAgg                                                                                                                               >  8:130566/1‑138 (MQ=255)
                         tgtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAgg                                                                                                                               <  7:130566/138‑1 (MQ=255)
                                   cagcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGt                                                                                                          <  5:25748/149‑1 (MQ=255)
                                                     catATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGc                                                                                        <  5:115797/149‑1 (MQ=255)
                                                     catATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGc                                                                                        <  6:150768/149‑1 (MQ=255)
                                                      atatCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCg                                                                                       >  7:64217/1‑149 (MQ=255)
                                                            aTAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACACCACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAagag                                                                                 >  8:36519/1‑149 (MQ=255)
                                                                  aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCt                                                                                            >  2:18101/1‑131 (MQ=255)
                                                                  aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCt                                                                                            <  1:18101/131‑1 (MQ=255)
                                                                                    aGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGc                                                          <  8:64217/148‑1 (MQ=255)
                                                                                          ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                        <  7:37517/128‑1 (MQ=255)
                                                                                          ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                        >  8:37517/1‑128 (MQ=255)
                                                                                            aGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGa                                                   >  4:113462/1‑147 (MQ=255)
                                                                                              gTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGtt                                                                                                              >  4:138864/1‑86 (MQ=255)
                                                                                              gTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGtt                                                                                                              <  3:138864/86‑1 (MQ=255)
                                                                                                 gATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAg                                                >  3:127445/1‑145 (MQ=255)
                                                                                                  aTTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGAt                                           <  7:36519/149‑1 (MQ=255)
                                                                                                                           gATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACtt                   >  4:70681/1‑148 (MQ=255)
                                                                                                                                cGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCtt             >  5:169613/1‑149 (MQ=255)
                                                                                                                                    cAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGc                                               >  4:135099/1‑111 (MQ=255)
                                                                                                                                    cAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGc                                               <  3:135099/111‑1 (MQ=255)
                                                                                                                                         aTCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTgaga    <  3:70681/149‑1 (MQ=255)
                                                                                                                                            aGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGGGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAgc  <  6:139220/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGC  >  NC_000913/3876870‑3877157

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