Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 469,786 A→G 100% G118G (GGA→GGG USA300HOU_0445 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730469,7860AG90.7% 47.1 / ‑2.3 21G118G (GGA→GGGUSA300HOU_0445hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (9/10);  total (10/11)
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.52e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAGATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAA  >  CP000730/469650‑469918
                                                                                                                                        |                                                                                                                                    
gataTGGCAAATTTAAAATGGGTTATAAAGGGATGGGGACCATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCat                                                                                                                                   <  1:262755/137‑1 (MQ=21)
  gaTGGCGAATTTAAACAGCGTGGTAAAGAGATGCGGACGATATATACAGATTGCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACatta                                                                                                                                 <  2:8995/140‑1 (MQ=255)
           tttAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCATAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTaa                                                                                                                        <  1:81014/140‑1 (MQ=255)
               aaaaGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTTGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACa                                                                                                                    <  1:82721/140‑1 (MQ=255)
                 aatggTGATAAAGGGATGTGGACGATATATACAAATTTCGCCGAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTATTATACTTAGATAGTAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACAtt                                                                                                                  <  2:136740/137‑1 (MQ=25)
                             gggATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTGTTTTGTTAAGAAATTCTATAATAAggg                                                                                                      >  1:39452/1‑140 (MQ=255)
                                   tgGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATACTATAATAAGGGCAAAtt                                                                                                >  1:102228/1‑140 (MQ=255)
                                     gACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTcc                                                                                              >  2:130694/1‑140 (MQ=255)
                                            tataCAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAg                                                                                       >  1:59305/1‑140 (MQ=255)
                                              taCAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGaa                                                                                     >  1:281385/1‑140 (MQ=255)
                                               aCAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAACTTCCCAGATAGaaa                                                                                    >  2:81885/1‑140 (MQ=255)
                                                             aaGTAATAAACAAGGTGGAGTGAGTAATGAAGGGATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtg                                                                      <  2:21719/140‑1 (MQ=255)
                                                                       cAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGaaaaa                                                            <  1:191025/140‑1 (MQ=255)
                                                                       cAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATACTATAATAAGATCAAATACCCAGATAGAAAAAACGATAAAGATAAAAAGaaaaa                                                            >  1:24887/1‑140 (MQ=255)
                                                                               aTTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATaaaa                                                    >  1:101855/1‑140 (MQ=255)
                                                                                ttGAGTTATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAAt                                                   <  1:226096/140‑1 (MQ=255)
                                                                                 tGAGTAATGAAGGTATTGTCTTATACTTAGATAGAAATACAAGGAATGCAAAGGGGCATTATTTTGTTAAGACATTCTATAAGAAATACAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAAtt                                                  <  1:139361/140‑1 (MQ=21)
                                                                                        tGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCttat                                           <  2:32558/140‑1 (MQ=255)
                                                                                            ggTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCCGATAGAAAAAATTATAAAGTTGAAATGAAAAATGATAAAATTATATTATTAGa                                       >  2:38679/1‑140 (MQ=255)
                                                                                                              gttagaAATACACGGACTGCAAAGGGGCATTATTTTGGTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATTATAAAATTATCTTATTAGATAAAGTAGAAGATACaaa                     <  2:102228/138‑1 (MQ=255)
                                                                                                                                 cAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGaa  >  2:136131/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GAGATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAA  >  CP000730/469650‑469918

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