Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,343,438 C→T 100% intergenic (+114/+250) USA300HOU_2200 → / ← alsD hypothetical protein/acetolactate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,343,4380CT87.0% 62.9 / 4.7 31intergenic (+114/+250)USA300HOU_2200/alsDhypothetical protein/acetolactate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (13/14);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTAA  >  CP000730/2343304‑2343572
                                                                                                                                      |                                                                                                                                      
aTACAAATTGTAGGTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATTGt                                                                                                                                  >  2:331616/1‑141 (MQ=255)
 tACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTa                                                                                                                                 >  2:5393/1‑141 (MQ=255)
  aCAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATGGTaa                                                                                                                                >  2:328013/1‑141 (MQ=255)
      aTTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAAttt                                                                                                                            <  1:250567/141‑1 (MQ=255)
                    ataatTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTa                                                                                                              <  2:242770/141‑1 (MQ=255)
                      aatTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTAcc                                                                                                            >  1:318308/1‑141 (MQ=255)
                              aaTAGTTTTATAGGGCAAGCGCTGAAAAATTTGGGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTta                                                                                                    <  2:287876/141‑1 (MQ=255)
                              aaTAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTta                                                                                                    >  2:103770/1‑141 (MQ=255)
                                tAGTTTTATAGGGCAAGAACCTCAAAGGGTGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTAGCACAATTtata                                                                                                  >  2:322822/1‑141 (MQ=255)
                                      tataGGGCAAGCGCTGAAAAATTGGCGTTTTTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTaa                                                                                            <  1:160120/141‑1 (MQ=255)
                                          gggCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACAt                                                                                        <  1:284284/141‑1 (MQ=255)
                                                      aaaaaTTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAAtct                                                                            <  2:218750/141‑1 (MQ=255)
                                                             ggCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTATCGTAGTAACAGATTGGCCTATGTAATAGCAAATTTTATAGGAGTGTGGCACCCAAATAAAAATCGAAACATACCAGCCAATCTCagaaaa                                                                     >  2:229408/1‑136 (MQ=255)
                                                                 tttGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTgg                                                                 >  1:98374/1‑141 (MQ=255)
                                                                       tttttACTTATTGACCTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTACATGAGTGTGTACCACAATTTATAATCTAAACATAACCGCAAATCTCAAATATCTGGAACCta                                                           <  2:154730/141‑1 (MQ=255)
                                                                        tcttAGGTGTTGAACTTGTAAACTATTGTTCTTCCCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAACTTTATATGAATGTCCACCACAATTTATAATCTAAACATAACAGCAAACCTCAAATATCTGAAACCtaa                                                          <  2:273761/139‑1 (MQ=255)
                                                                          ttAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCtaata                                                        >  1:15214/1‑141 (MQ=255)
                                                                             gTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAAt                                                     <  1:333869/141‑1 (MQ=255)
                                                                                aTTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATACATccc                                                  >  2:277692/1‑141 (MQ=255)
                                                                                  tGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCt                                                <  1:5393/141‑1 (MQ=255)
                                                                                   gAACTTGTAAATTATTGTTGTTTTCTTCTCGTAGGAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAAGCCCCTc                                               <  2:22855/141‑1 (MQ=255)
                                                                                   gAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTc                                               <  2:130907/141‑1 (MQ=255)
                                                                                                    ttgttTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGa                              >  1:134715/1‑141 (MQ=255)
                                                                                                         tGCTTCTCGTAGTAACCGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGa                         <  1:307396/141‑1 (MQ=255)
                                                                                                             tctcGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGattta                     <  2:98374/141‑1 (MQ=255)
                                                                                                                  aagtaACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGAtttatttta                <  1:104374/140‑1 (MQ=255)
                                                                                                                     taACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTttatta             >  1:56694/1‑141 (MQ=255)
                                                                                                                        cAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACtt          <  1:1510/141‑1 (MQ=255)
                                                                                                                           aTTGGCCTATATAATACTACAGTTTATATGAGTGTGTACCACACTTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTCTTTTATTACTTAGa       >  1:171310/1‑141 (MQ=255)
                                                                                                                               gCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTa   <  1:264257/141‑1 (MQ=255)
                                                                                                                                ccTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTaa  >  1:34007/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                      
ATACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTAA  >  CP000730/2343304‑2343572

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