Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,024 G→A 12.2% intergenic (+40/+42) ychS → / ← rttR putative uncharacterized protein YchS/small RNA RttR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0240GA12.2% 101.3 / 6.0 41intergenic (+40/+42)ychS/rttRputative uncharacterized protein YchS/small RNA RttR
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base A (2/3);  total (18/23)
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

TGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGG  >  NC_000913/1286879‑1287162
                                                                                                                                                 |                                                                                                                                          
tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                                                                        >  1:1320243/1‑150 (MQ=39)
   aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                                                                     <  2:1320243/150‑1 (MQ=39)
               tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                                                                         <  2:1398220/150‑1 (MQ=255)
                   ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                                                                     <  1:1138121/150‑1 (MQ=255)
                      aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                                  >  2:180345/1‑150 (MQ=255)
                       ggatgagaagcttcgaccaaggttcgactcgagcgccagcgagagagcgttgccgcaggcaacgacccgaagggcgaagcgcgcagcgctgagtaATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                                                                 <  1:1254975‑M1/55‑1 (MQ=255)
                         aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAc                                                                                                               <  2:136829/150‑1 (MQ=255)
                               cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                                                                      <  2:1551718/121‑1 (MQ=40)
                               cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                                                                      >  1:1551718/1‑121 (MQ=40)
                                  gggtgggcgcTTCGCCCTTCGGGTCGTTGACTGCGGCAAGGCTCTCTCTCTGTCTCTCGAGCCGAACCTTTGGCGTAGCTGCTCATCGTTCCCCGCTTGGGCATAATATTTGATTGCGGATTCTCTTGAGAGTTCAGGGACTTTTGAAAg                                                                                                      <  2:1089742/144‑1 (MQ=14)
                                     tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                                                   <  1:988873/150‑1 (MQ=255)
                                                cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                                                                        <  1:789872/150‑1 (MQ=255)
                                                 actcgagcgccagcgagagagcgttgccgcaggcaacgacccgaagggcgaagcgcgcagcgctgagtaATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATTGTGGTggggg                                                                                       >  2:971001‑M1/70‑150 (MQ=255)
                                                     cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                                                                            >  2:340517/1‑93 (MQ=255)
                                                     cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                                                                            <  1:340517/93‑1 (MQ=255)
                                                          ccagcgagagagcgttgccgcaggcaacgacccgaagggcgaagcgcgcagcgctgagtaATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                                                                              >  2:1254975‑M1/61‑150 (MQ=255)
                                                           cagcgagagagcgttgccgcaggcaacgacccgaagggcgaagcgcgcagcgctgagtaATCCTTCCCCGCATGTGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                             >  2:507188‑M1/60‑150 (MQ=255)
                                                              gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                                                                          >  2:1134937/1‑150 (MQ=255)
                                                                 ggCGGCAACGCTCTCTCGCTGCCGCGCCAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                       <  2:886360/149‑1 (MQ=14)
                                                                             ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                                                                >  1:408081/1‑81 (MQ=255)
                                                                             ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                  <  1:604930/127‑1 (MQ=255)
                                                                             ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                  >  2:604930/1‑127 (MQ=255)
                                                                             ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCg                                                                                                                                <  2:408081/81‑1 (MQ=255)
                                                                                 cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                                                                                         <  1:890881/68‑1 (MQ=255)
                                                                                 cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                                                                                         >  2:890881/1‑68 (MQ=255)
                                                                                   cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                                     <  1:1134937/150‑1 (MQ=255)
                                                                                      ccGCTCCAGCCGAACCTTCGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTGATTGCGTATTCTCTTTAGAGTTCAGGTCCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGATTCGCCCTTCGGGt                                                  <  2:828476/149‑1 (MQ=14)
                                                                                       cGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAtt                                                                                                                                  <  2:1155395/69‑1 (MQ=25)
                                                                                       cGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAtt                                                                                                                                  >  1:1155395/1‑69 (MQ=25)
                                                                                                       tAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACg                                 <  1:180345/150‑1 (MQ=255)
                                                                                                          tCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                                                                          >  2:892823/1‑74 (MQ=255)
                                                                                                          tCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                                                                          <  1:892823/74‑1 (MQ=255)
                                                                                                          tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                                     >  1:65720/1‑127 (MQ=255)
                                                                                                          tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                                     <  2:65720/127‑1 (MQ=255)
                                                                                                           cGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                             >  2:359464/1‑150 (MQ=255)
                                                                                                                   ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                     >  1:987770/1‑150 (MQ=255)
                                                                                                                     cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                   <  1:359464/150‑1 (MQ=255)
                                                                                                                                 cATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTGAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                      <  2:808665/135‑1 (MQ=255)
                                                                                                                                 cATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                      >  1:808665/1‑135 (MQ=255)
                                                                                                                                     ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                 >  1:1068552/1‑72 (MQ=255)
                                                                                                                                     ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                 <  2:1068552/72‑1 (MQ=255)
                                                                                                                                     ggCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGATACTTTTAAAATAGATTATGGTAAGGGAATGATTACTAAGAGCTGAGCGATTAGCCCTTATGGTAGTTAAATGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg   >  2:375186/1‑150 (MQ=9)
                                                                                                                                      gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg  <  1:614884/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
TGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGG  >  NC_000913/1286879‑1287162

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