Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,843,140 G→A 100% Q566* (CAG→TAG)  gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,1400GA96.2% 78.0 / ‑4.2 26Q566* (CAG→TAG) gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/12);  total (13/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 = 8.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG  >  NC_002947/1842993‑1843283
                                                                                                                                                   |                                                                                                                                               
tCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTaca                                                                                                                                               <  1:820900/150‑1 (MQ=255)
     cacaGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAAc                                                                                                                                          <  2:842174/150‑1 (MQ=255)
     cacaGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAAc                                                                                                                                          <  1:103260/150‑1 (MQ=255)
      aaaGCACTACGCATTGGCAGTGGTAACGCTCGATTCCCCGATGATCGTGGCTCGGTTGCTCGGGCGATAGGTTACCCCGGTTGGCACCGCGAACCGCGCACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACTCCCTACAGCAACt                                                                                                                                         <  2:240769/148‑1 (MQ=255)
               cGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagg                                                                                                                                <  1:696711/150‑1 (MQ=255)
                  aTTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagga                                                                                                                               >  1:306064/1‑148 (MQ=255)
                  aTTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagga                                                                                                                               <  2:306064/148‑1 (MQ=255)
                                 gCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGACCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGc                                                                                                              >  2:389529/1‑150 (MQ=255)
                                                 atAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTc                                                                                              <  1:131137/150‑1 (MQ=255)
                                                                   gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc                                                                                                        >  1:31219/1‑122 (MQ=255)
                                                                   gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc                                                                                                        <  2:31219/122‑1 (MQ=255)
                                                                       gACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCAAACCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTgg                                                                        >  1:62081/1‑150 (MQ=255)
                                                                            gTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGgcgcgc                                                                   >  2:293130/1‑150 (MQ=255)
                                                                                 ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt                                                              >  2:772461/1‑150 (MQ=255)
                                                                                 ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt                                                              >  1:327979/1‑150 (MQ=255)
                                                                                     gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTg                                                          <  2:327979/150‑1 (MQ=255)
                                                                                      ttGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTgg                                                         <  1:293130/150‑1 (MQ=255)
                                                                                                    cgcgAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATggc                                           >  2:379037/1‑150 (MQ=255)
                                                                                                           tCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcgg                                    <  2:557116/150‑1 (MQ=255)
                                                                                                             gAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTTCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcggcg                                  >  1:401790/1‑150 (MQ=255)
                                                                                                                  gCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg                                                                                        >  1:521146/1‑91 (MQ=255)
                                                                                                                  gCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg                                                                                        <  2:521146/91‑1 (MQ=255)
                                                                                                                                        cTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGt               >  1:474035/1‑142 (MQ=255)
                                                                                                                                        cTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGt               <  2:474035/142‑1 (MQ=255)
                                                                                                                                             cccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTAACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCgg  >  1:629911/1‑150 (MQ=255)
                                                                                                                                              ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTgg             >  1:753951/1‑138 (MQ=255)
                                                                                                                                              ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTgg             <  2:753951/138‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG  >  NC_002947/1842993‑1843283

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