Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,018,612 G→A 100% R103* (CGA→TGA)  USA300HOU_1869 ← transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,018,6120GA97.6% 110.3 / ‑5.2 37R103* (CGA→TGA) USA300HOU_1869transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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.

ATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACA  >  CP000730/2018481‑2018737
                                                                                                                                   |                                                                                                                             
aTTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTc                                                                                                                      <  2:2762/141‑1 (MQ=255)
aTTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGGAGAATCATTCATCTTc                                                                                                                      <  1:234736/141‑1 (MQ=255)
 ttATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTAATCTTcc                                                                                                                     >  2:237207/1‑141 (MQ=255)
   aTCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATGATTCATCTTCCAt                                                                                                                   >  2:284267/1‑141 (MQ=255)
   aTCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCAt                                                                                                                   <  2:30963/141‑1 (MQ=255)
            ataatTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTTCATTGAAATaaa                                                                                                          >  1:246561/1‑141 (MQ=255)
            ataatTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATaaa                                                                                                          >  2:22183/1‑141 (MQ=255)
              aatTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAgg                                                                                                        <  1:123306/141‑1 (MQ=255)
               atTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAACTAAAGGc                                                                                                       >  1:194479/1‑141 (MQ=255)
               atTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGc                                                                                                       >  1:43106/1‑141 (MQ=255)
               atTTTTGACAACTTGCTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGc                                                                                                       >  1:292994/1‑141 (MQ=255)
                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGATGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAATGCa                                                                                                      >  2:42243/1‑141 (MQ=255)
                         aCTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACa                                                                                             >  1:207707/1‑141 (MQ=255)
                          cTTGTTGTGTAATGGTACTGCCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACAt                                                                                            >  1:295523/1‑141 (MQ=255)
                                       ggTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATGCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGtt                                                                               >  1:73611/1‑141 (MQ=255)
                                         tactacCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTAt                                                                             <  1:232272/141‑1 (MQ=255)
                                              ccaccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCt                                                                        <  1:173096/141‑1 (MQ=255)
                                              ccaccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCt                                                                        >  1:270453/1‑141 (MQ=255)
                                                 ccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGac                                                                     >  1:208416/1‑141 (MQ=255)
                                                        caTCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAac                                                              <  1:157487/141‑1 (MQ=255)
                                                          tctcTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAacta                                                            >  2:34090/1‑141 (MQ=255)
                                                          tctcTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAacta                                                            >  1:85693/1‑141 (MQ=255)
                                                                          ttGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTCCTTTTTTTTTCAAttt                                            >  2:329410/1‑141 (MQ=255)
                                                                            caaaaTAACGCTCTATTTGTCCCTTTCAACTCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACACTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAAttttt                                          <  2:255419/140‑1 (MQ=255)
                                                                                 tAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAg                                     <  1:343727/141‑1 (MQ=255)
                                                                                 tAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAg                                     <  1:389331/141‑1 (MQ=255)
                                                                                   aaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGtt                                   >  1:198771/1‑141 (MQ=255)
                                                                                   aaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTATATTTTCAATTTTTCTTAGtt                                   >  1:221677/1‑141 (MQ=255)
                                                                                            ttGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACAt                          <  2:367141/141‑1 (MQ=255)
                                                                                            ttGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACAt                          <  2:112872/141‑1 (MQ=255)
                                                                                            ctGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACAt                          <  2:235057/140‑1 (MQ=255)
                                                                                                    ttCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAGCTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCgt                  >  1:131815/1‑141 (MQ=255)
                                                                                                         aTCGTATACTTGATTATTGTCGAAGCGTGCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGa             <  2:295523/141‑1 (MQ=255)
                                                                                                          tCGAATCCATGATGATTGTAGAATCATTCATCTTCCAGTGAAATATCGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGAt            <  1:89723/141‑1 (MQ=255)
                                                                                                                  atgatgATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCCCCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATa    <  1:378946/141‑1 (MQ=255)
                                                                                                                   tgatgaTTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATAc   >  1:179646/1‑141 (MQ=255)
                                                                                                                    tatgatTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATTCTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACa  <  2:317323/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
ATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACA  >  CP000730/2018481‑2018737

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: