Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,615,390 A→G 7.7% S815S (AGT→AGC)  USA300HOU_2480 ← helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,615,3900AG 7.7% 141.1 / 3.3 52S815S (AGT→AGC) USA300HOU_2480helicase
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base G (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

AAAAATTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTAGAGCATTAGTA  >  CP000730/2615253‑2615528
                                                                                                                                         |                                                                                                                                          
aaaaaTTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTa                                                                                                                                         >  2:218164/1‑141 (MQ=255)
 aaaaTTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTaa                                                                                                                                        >  2:3079/1‑141 (MQ=255)
 aaaaTTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTaa                                                                                                                                        >  2:350156/1‑141 (MQ=255)
   aaTTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTaaaa                                                                                                                                      <  1:128686/141‑1 (MQ=255)
               tatTTTGTGCTATCACTNATATCCTCATGCTTATCATNAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATTTTCTCTTGaa                                                                                                                          >  1:295017/1‑141 (MQ=255)
                     gtgCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAAtattta                                                                                                                    >  2:57319/1‑141 (MQ=255)
                        cTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATttatta                                                                                                                 >  1:339209/1‑141 (MQ=255)
                        cTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATttatta                                                                                                                 >  1:108266/1‑141 (MQ=255)
                         tATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTtattat                                                                                                                <  1:313130/141‑1 (MQ=255)
                                 atatCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATATTATAAAATc                                                                                                        >  2:159622/1‑141 (MQ=255)
                                  tatCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCa                                                                                                       >  2:247237/1‑141 (MQ=255)
                                     taagcatgaggataTAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTCAAATCTTCTCTTGAATATTTATTATATAAAATCCGTc                                                                                                    >  2:23350/14‑141 (MQ=255)
                                          tGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTCTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTGCTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCAttt                                                                                               >  1:352674/1‑141 (MQ=255)
                                          tGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTCTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTGCTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCAttt                                                                                               >  1:19274/1‑141 (MQ=255)
                                            cTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCAttttt                                                                                             >  2:109326/1‑141 (MQ=255)
                                            cTTATCATAAGTCATAAATCTTGACATTTCTTAGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCAttttt                                                                                             <  2:433506/141‑1 (MQ=255)
                                                 cATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCtt                                                                                        >  2:115658/1‑141 (MQ=255)
                                                  aTAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCttt                                                                                       <  2:254392/141‑1 (MQ=255)
                                                      gTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTgg                                                                                   <  1:218164/141‑1 (MQ=255)
                                                        gATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTa                                                                                 <  1:298736/141‑1 (MQ=255)
                                                          tAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTAtt                                                                               <  1:311333/141‑1 (MQ=255)
                                                            aaTATTGGCATTTCTAGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGTtttttt                                                                             >  2:81494/1‑141 (MQ=255)
                                                               atTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCa                                                                          <  1:295776/141‑1 (MQ=255)
                                                                ttGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTTGTATTTTTCAt                                                                         >  1:409229/1‑141 (MQ=255)
                                                                ttGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCAt                                                                         >  1:90181/1‑141 (MQ=255)
                                                                  ggCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTg                                                                       >  1:340299/1‑141 (MQ=255)
                                                                     aTTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGttt                                                                    >  1:273932/1‑140 (MQ=255)
                                                                        tCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTGTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGTTTTTTTTCTTTGTTTTAc                                                                 >  2:437970/1‑141 (MQ=255)
                                                                            ggCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTg                                                             >  1:393851/1‑141 (MQ=255)
                                                                             gCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGa                                                            <  1:205369/141‑1 (MQ=255)
                                                                               tttttATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATa                                                          <  2:130134/141‑1 (MQ=255)
                                                                                    aTCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCTTTTTTATACTTTGATAATTCa                                                     >  2:197898/1‑141 (MQ=255)
                                                                                          tATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTGCTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGa                                               <  2:352674/141‑1 (MQ=255)
                                                                                                  gATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCCATAAGATCTTCTac                                       <  2:172789/141‑2 (MQ=255)
                                                                                                  gATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTaa                                       <  2:176532/141‑1 (MQ=255)
                                                                                                  gATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTaa                                       >  1:282984/1‑141 (MQ=255)
                                                                                                     tACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTaaaaa                                    <  1:174125/141‑1 (MQ=255)
                                                                                                        aGACGAACTATTTTTATTCCAATTAAATATTTTACTAATATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTATTACCTTTTTGGTATTTTTCACTGTTATACTTTGATATTTCACTAAGATCTTCTACAAAcct                                 <  2:207097/141‑3 (MQ=255)
                                                                                                            gAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTa                             <  1:64686/141‑1 (MQ=255)
                                                                                                             aaCCATTTTTATTCCAATTAAATATTTTGCTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAATATTTTaa                            <  2:19274/141‑1 (MQ=255)
                                                                                                                  tttttATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGtt                       >  2:219476/1‑141 (MQ=255)
                                                                                                                   ttttATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGttt                      <  2:288935/141‑1 (MQ=255)
                                                                                                                      tATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGAt                   >  1:390502/1‑141 (MQ=255)
                                                                                                                       aTTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGAtt                  >  2:213331/1‑141 (MQ=255)
                                                                                                                          ccAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACt               <  1:147054/141‑1 (MQ=255)
                                                                                                                          ccAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACt               <  1:457264/141‑1 (MQ=255)
                                                                                                                           cAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACtt              <  2:83505/141‑1 (MQ=255)
                                                                                                                              ttAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTTTTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTata           >  2:442514/1‑139 (MQ=255)
                                                                                                                                  atatTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTAGAGCAt       >  2:301574/1‑141 (MQ=255)
                                                                                                                                     ttttACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTAATTTTCATTGTTATACTTTGATAATTTAATAAGATTTTCTAAAAACATATTAAATTTTTGATTAATTATAGCATTat    >  2:143994/1‑140 (MQ=255)
                                                                                                                                      tttACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTAGAGCATTAGt   <  1:437100/141‑1 (MQ=255)
                                                                                                                                       ttACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTAGAGCATTAGTa  <  1:304925/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
AAAAATTCATCTTCATATTTTGTGCTATCACTTATATCCTCATGCTTATCATAAGTGATAAATATTGGCATTTCTTGGCTTTTTATCATGTATCCCATGATTACAGACGAACCATTTTTATTCCAATTAAATATTTTACTAAAATCTTCTCTTGAATATTTATTATATAAAATCAGTCCATTTTTACCTTTTTGGTATTTTTCATTGTTATACTTTGATAATTCAATAAGATCTTCTAAAAACATATTAAATGTTTGATTACTTAGAGCATTAGTA  >  CP000730/2615253‑2615528

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