Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 707,910 A→G 100% intergenic (‑1/‑121) USA300HOU_0653 ← / → USA300HOU_0654 ABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730707,9100AG96.6% 91.0 / ‑4.5 29intergenic (‑1/‑121)USA300HOU_0653/USA300HOU_0654ABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/17);  minor base T (1/0);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAA  >  CP000730/707773‑708045
                                                                                                                                         |                                                                                                                                       
ttGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGt                                                                                                                                      <  1:45415/141‑1 (MQ=255)
    tAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtat                                                                                                                                  >  2:186255/1‑141 (MQ=255)
     aaGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtata                                                                                                                                 >  2:320557/1‑141 (MQ=255)
             aTTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTaa                                                                                                                         <  1:126188/141‑1 (MQ=255)
                tACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGCTAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTa                                                                                                                      >  1:122027/1‑141 (MQ=255)
                 aCCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTaa                                                                                                                     <  1:64560/141‑1 (MQ=255)
                        ccATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTa                                                                                                              <  1:345477/141‑1 (MQ=255)
                           ttCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGc                                                                                                           >  2:328675/1‑141 (MQ=255)
                                gCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAAcc                                                                                                      <  2:207942/141‑1 (MQ=255)
                                               tattTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTataa                                                                                       >  1:331932/1‑141 (MQ=255)
                                                        gcgTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTgag                                                                              <  2:16139/141‑1 (MQ=255)
                                                            aCTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAATCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAGGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAata                                                                          >  2:167732/1‑141 (MQ=255)
                                                                            agaTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTa                                                          <  2:41248/139‑1 (MQ=255)
                                                                              cggtgtttAAGTACATGCCGATTAGCTAAATACAGGTGTAAATCGTTTGTTGCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTTACATGTCAAGTaaa                                                        >  2:392609/6‑141 (MQ=255)
                                                                               tGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAAt                                                       >  2:46046/1‑141 (MQ=255)
                                                                               tGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATATATCGAAAGTGTTTATAGTTTAAATCATCTAAGAAAAAGTTTAGGCTAACTTAAATAATTGTATGATAAACTGAGAATATTTATCATGTCAAGTAAAt                                                       >  1:189949/1‑141 (MQ=255)
                                                                                      aaGTACATGCCTATTACCAAAAAACAGTTTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCCAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGat                                                <  1:190908/141‑1 (MQ=255)
                                                                                           cATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataat                                           <  1:320557/141‑1 (MQ=255)
                                                                                                     aCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatg                                 >  1:174538/1‑141 (MQ=255)
                                                                                                          aaaaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATCTATCATGTCAAGTAACTTCGTGATATAATATAGACAatgtatgt                            <  2:350322/141‑1 (MQ=255)
                                                                                                          aaaaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTGATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgtatgt                            <  1:30170/141‑1 (MQ=255)
                                                                                                          aaaaCAGATTTAAATCTTTTGGTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgtatgt                            <  1:14032/141‑1 (MQ=255)
                                                                                                             aCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGCTATACTATAGACAATGTATGTGAg                         >  2:126833/1‑141 (MQ=255)
                                                                                                                    tAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTGAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAg                  <  1:88785/141‑1 (MQ=255)
                                                                                                                    tAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAg                  <  1:125540/141‑1 (MQ=255)
                                                                                                                        tCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATg              <  1:46046/141‑1 (MQ=255)
                                                                                                                                 tcgaaCAAGCGTTTAGACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaa     <  1:346862/138‑1 (MQ=255)
                                                                                                                                   tAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaaga   >  1:107117/1‑141 (MQ=255)
                                                                                                                                    aaCAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaagaa  <  2:56677/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAA  >  CP000730/707773‑708045

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