Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 367,644 G→T 100% intergenic (‑96/‑20) nanA ← / → USA300HOU_0336 N‑acetylneuraminate lyase/possible transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730367,6440GT95.7% 82.8 / ‑5.5 29intergenic (‑96/‑20)nanA/USA300HOU_0336N‑acetylneuraminate lyase/possible transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (12/15);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGTT  >  CP000730/367512‑367779
                                                                                                                                    |                                                                                                                                       
gTAACGCTGCATATCAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAGACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAttt                                                                                                                                 >  2:191969/1‑141 (MQ=255)
   actcttcATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATg                                                                                                                              <  2:237313/135‑1 (MQ=255)
    cGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGt                                                                                                                             <  2:14279/141‑1 (MQ=255)
                  cTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGaaaa                                                                                                               <  1:231295/141‑1 (MQ=255)
                       aaaTCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATAttt                                                                                                          >  2:133437/1‑141 (MQ=255)
                         aTCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTAATATTTATGTATAATTAATGAAAATGttttt                                                                                                        >  1:105379/1‑141 (MQ=255)
                          tCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTc                                                                                                       <  2:245958/141‑1 (MQ=255)
                                   aTTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAAt                                                                                              <  2:10826/141‑1 (MQ=255)
                                   aTTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAAt                                                                                              <  1:45737/141‑1 (MQ=255)
                                      tttaaGTCACTCCAAATATCATTTGATAATATAGTTTATAGCTATAATTGTAAACGCTATTATAAAATGTAACAATATCTTTTTGACAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGa                                                                                           <  2:35866/138‑1 (MQ=255)
                                             cccTCCAAAAATCATTTGATAATATTGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgag                                                                                    <  1:252992/141‑1 (MQ=255)
                                                   aaataTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTTATATTTATGTATAATTGATGAAAATGTTTTTCAAAATCAATAGAAATTGAGTGAGTa                                                                              <  2:54829/141‑1 (MQ=255)
                                                   aaaaaTCATTTGATAATATAGTTTACAACTATAATTGTAAACGCTATCATAAAATGTAACAATATTTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTa                                                                              <  1:159173/141‑1 (MQ=255)
                                                     aaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAg                                                                            >  1:67899/1‑141 (MQ=255)
                                                       aTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTTTATATTTTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGGGGGAGTAAggg                                                                          >  2:240479/1‑140 (MQ=255)
                                                            ttGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTAtt                                                                     <  1:194870/141‑1 (MQ=255)
                                                                                  taATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAgg                                               >  1:251259/1‑141 (MQ=255)
                                                                                        tAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCa                                         >  1:29021/1‑141 (MQ=255)
                                                                                        tAAAAGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAAATCATGAAAATGTTTTTTAAAATCAATAGAAATTGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCa                                         <  2:42528/141‑1 (MQ=255)
                                                                                                aTCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTaaca                                 >  1:138349/1‑139 (MQ=255)
                                                                                                   aTAAAATGTAACAATATCTTTTTGAAAATTGTATTCATGTTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCgg                              >  2:170213/1‑141 (MQ=255)
                                                                                                    aaaaaTGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGc                             <  2:129248/140‑1 (MQ=255)
                                                                                                                ataaCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACCTCGCAATCGATCTTGCAGGCCATCACCTTAAATCGGCAGTTATTGATaa                 <  2:163977/141‑1 (MQ=255)
                                                                                                                 tatCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAg                >  2:59545/1‑141 (MQ=255)
                                                                                                                   tCTTTTTGAAAATTGTATTTATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCa              <  1:234160/141‑1 (MQ=255)
                                                                                                                     tttttGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAAGCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGAAGTAATTTATACGCact            <  2:129767/141‑3 (MQ=255)
                                                                                                                          gAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAat       <  1:129408/141‑1 (MQ=255)
                                                                                                                              aTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGt   >  2:237595/1‑141 (MQ=255)
                                                                                                                               ttGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGtt  >  1:208778/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGTT  >  CP000730/367512‑367779

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