Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,613,762 T→C 100% intergenic (‑31/‑37) scpA ← / → USA300HOU_1497 possible segregation and condensation protein A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,613,7620TC99.4% 86.7 / ‑6.8 27intergenic (‑31/‑37)scpA/USA300HOU_1497possible segregation and condensation protein A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/12);  minor base G (1/0);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAAATATCTATTTCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCATGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGAA  >  CP000730/1613625‑1613894
                                                                                                                                         |                                                                                                                                    
tAAATATCTATTTCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGaa                                                                                                                                   >  1:136899/1‑141 (MQ=255)
 aaaTATCTATTTCAAATTTCTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGaaa                                                                                                                                  >  1:237104/1‑141 (MQ=255)
                  ttttGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAAt                                                                                                                 >  2:233111/1‑141 (MQ=255)
                       gATAAGATGCAGCAATCAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATCATATCTACTTCGGATTGCGTCATATGTGCATAATTCTCTCAAAAATTATAGCAGGAAAATAACTTACTTTTAAATTgggg                                                                                                            >  2:134960/1‑141 (MQ=255)
                             aTGCAGCAATAAATCTAATGGTCCGTTGAAAGCATCTAATTTAACTTCATACAAAACATCTACCTCGTATTGCGGCATATGTGCATAATACTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAAttaaaaaaaaaa                                                                                                      <  2:161053/141‑11 (MQ=255)
                             aTGAAGCAATAAAACTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCCCCCAAATATTATATCACTAAAATAACTTACTTTTAAATTGGGGTGATaa                                                                                                      <  1:233111/141‑1 (MQ=255)
                              tGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAAc                                                                                                     <  2:388561/141‑1 (MQ=255)
                               gcagcaATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACa                                                                                                    <  2:209275/141‑1 (MQ=255)
                               gcagcaATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACa                                                                                                    <  2:201834/141‑1 (MQ=255)
                                  gaaaTAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACaaaa                                                                                                 <  1:176214/139‑1 (MQ=255)
                                     aTACATCTAATGGCCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGc                                                                                              <  2:258962/141‑1 (MQ=255)
                                     aTAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGc                                                                                              >  1:314706/1‑141 (MQ=255)
                                                   ccATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATa                                                                                >  2:257317/1‑141 (MQ=255)
                                                                          ttCATCCATAAGATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATAc                                                         <  1:16756/141‑1 (MQ=255)
                                                                                            cTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCt                                       <  1:40440/141‑1 (MQ=255)
                                                                                                   tGCGTCATATGTGCATGAGTCTCTAAGATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACCTTATTTCCTATGTCAt                                <  2:97829/141‑1 (MQ=255)
                                                                                                    gCGTAATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATg                               >  1:269112/1‑141 (MQ=255)
                                                                                                      gTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGAc                             <  1:39860/141‑1 (MQ=255)
                                                                                                       tCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACa                            >  1:327048/1‑141 (MQ=255)
                                                                                                           atGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACAtttt                        <  2:314706/141‑1 (MQ=255)
                                                                                                                cATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTagaag                   >  1:184379/1‑141 (MQ=255)
                                                                                                                        tctAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATgg           >  2:33714/1‑141 (MQ=255)
                                                                                                                          tAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGAGAACAAAATGCACAAAGCACAGATAAATTTCTAATATCAAATTCCTAATAAGAAAAATTATTTCCTATGTCAACACATTTTagaaattgattggaa         >  2:192301/1‑130 (MQ=255)
                                                                                                                             aaTATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGc      >  2:147027/1‑141 (MQ=255)
                                                                                                                                 ttatAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGTAGATGCATGGAAAGCTGaa  >  1:314914/1‑141 (MQ=255)
                                                                                                                                 ttatAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGaa  >  1:36927/1‑141 (MQ=255)
                                                                                                                                 ttatAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGaa  >  2:308012/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TAAATATCTATTTCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCATGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGAA  >  CP000730/1613625‑1613894

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: