Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,079,114 T→C 69.9% H201R (CAC→CGC)  USA300HOU_1936 ← ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,079,1140TC69.9% 40.6 / 19.2 33H201R (CAC→CGC) USA300HOU_1936ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base C (17/6);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

GTACGCTTATGTTCTTCATAAATAAGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGATCAAGACCTGAAG  >  CP000730/2078979‑2079249
                                                                                                                                       |                                                                                                                                       
gTACGCTTATGTTGTTCATAAATATGGAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATg                                                                                                                                    >  2:191627/1‑141 (MQ=255)
                   aaataaGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACAACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCaaccctcgtgt                                                                                                                 >  2:395677/1‑131 (MQ=255)
                    aataaGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTc                                                                                                                >  2:279865/1‑141 (MQ=255)
                      taaGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTa                                                                                                              >  2:4815/1‑141 (MQ=255)
                        aGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTaaa                                                                                                            >  1:219229/1‑141 (MQ=255)
                          tAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACATTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCGAAATc                                                                                                          >  1:82831/1‑141 (MQ=255)
                          tAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATc                                                                                                          >  1:324135/1‑141 (MQ=255)
                           agagATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCt                                                                                                         >  1:257563/1‑141 (MQ=255)
                            gagaTGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTg                                                                                                        <  2:51208/141‑1 (MQ=255)
                             agaTGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTga                                                                                                       <  2:282187/141‑1 (MQ=255)
                                  cTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATaa                                                                                                  >  2:30007/1‑141 (MQ=255)
                                     aTTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCCCCGTCACTTAAGCGCATGATATCATCAGCTATTTTTTCTAAATCTGAGATAAtgt                                                                                               >  1:52200/1‑141 (MQ=255)
                                     aTTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGACATacatt                                                                                               >  2:275672/1‑137 (MQ=255)
                                       tcatcatcaAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAAtgtgt                                                                                             >  1:245294/1‑141 (MQ=255)
                                               aaGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACAt                                                                                     <  1:56494/141‑1 (MQ=255)
                                                    ccTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAata                                                                                >  1:258375/1‑141 (MQ=255)
                                                        aaTAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGt                                                                            <  1:385662/141‑1 (MQ=255)
                                                              accacTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATATGACATAATGTGTGTCGACATGAATATTGTTTTAtt                                                                      >  2:191489/1‑141 (MQ=255)
                                                                      caACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATGATGTGTGTCGACATGCATATTGTTTTATTTTCATaaa                                                              >  2:16154/1‑141 (MQ=255)
                                                                      caACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGACTATTGTTTTATTTTCATCaa                                                              <  1:239238/141‑1 (MQ=255)
                                                                                tCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTg                                                    >  1:136181/1‑141 (MQ=255)
                                                                                tCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGGCGACATGAATATTGTTTTATTTTCATAAATTAATTATTg                                                    >  2:179825/1‑141 (MQ=255)
                                                                                  tttGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCt                                                  <  2:308343/141‑1 (MQ=255)
                                                                                    tGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGa                                                >  1:178723/1‑141 (MQ=255)
                                                                                        aCAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTTTTTGTTTTTGATAATGTGTGTCGACATGATTATTGTTTTTTTTTATTGGATTTATTCTTGTTGTAttt                                            >  2:286035/1‑140 (MQ=255)
                                                                                                        aCCATTTAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCAtt                            >  1:266314/1‑141 (MQ=255)
                                                                                                            tttAAGATAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTTTTTCTCTGATAATGTTTGTCGACATGTATATTGTTTTATTTTCGTCACTTAATTATTGTTGAGTTTTTGCCATTAGTTCATTTCtt                        >  2:206316/1‑140 (MQ=255)
                                                                                                                 taTAATTTCACCGTTACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTaa                   <  2:193949/140‑1 (MQ=255)
                                                                                                                   tAACTTCACCGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAgg                 <  2:393101/141‑1 (MQ=255)
                                                                                                                     aCTTCACCGTAACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGa               <  1:223029/141‑1 (MQ=255)
                                                                                                                            cGTCACTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGATCAAGAc        >  2:66543/1‑141 (MQ=255)
                                                                                                                                aCTTAAGCGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGATCAAGACCTGa    >  1:27100/1‑141 (MQ=255)
                                                                                                                                  ttAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGATCAAGACCTGAAg  <  1:72392/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GTACGCTTATGTTCTTCATAAATAAGTAGAGATGCTAATTCATCATCAAGGTCCTCAATAGCACCACTCACAACTTGATATCTTTGTAACAGTTGTTCTTTTGAACCATTTAAGATAACTTCACCGTCACTTAAGTGAATGATATAATCAGCTATTTTTTCTAAATCTGAGATAATGTGTGTCGACATGAATATTGTTTTATTTTCATCAATTAATTCTTGCTGAATTATCTCCAATAGCTCATTTCTAGCTAAGGGATCAAGACCTGAAG  >  CP000730/2078979‑2079249

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: