Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,909,214 C→A 54.7% Y7* (TAC→TAA USA300HOU_1771 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,909,2140CA54.7% ‑5.7 / 56.6 40Y7* (TAC→TAAUSA300HOU_1771hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base A (16/6);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAATTTTGTCTAT  >  CP000730/1909076‑1909354
                                                                                                                                          |                                                                                                                                            
tttttttttttttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAAtt                                                                                                                                            >  2:472663/10‑141 (MQ=255)
  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt                                                                                                                                          >  2:187205/1‑141 (MQ=255)
  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGAGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATTATATGATATGTAAATGATAATCATCCAGGAGGTCTAGTATGTCTTGTTCAAAAAAATAAtttt                                                                                                                                          >  1:204484/1‑141 (MQ=255)
     tataTTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACt                                                                                                                                       <  1:281724/141‑1 (MQ=255)
          tttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCt                                                                                                                                  <  1:322512/141‑1 (MQ=255)
              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCt                                                                                                                              >  2:316139/1‑141 (MQ=255)
              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATATAGCt                                                                                                                              >  1:343300/1‑141 (MQ=255)
                      aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGATGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGatt                                                                                                                      >  2:222925/1‑141 (MQ=255)
                         ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGattatt                                                                                                                   <  2:451795/141‑1 (MQ=255)
                                   ttGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCtttt                                                                                                         <  2:96968/141‑1 (MQ=255)
                                        tCTAAAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTAAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCttt                                                                                                    <  2:20242/141‑1 (MQ=255)
                                                ttAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACaaa                                                                                            <  1:89080/141‑1 (MQ=255)
                                                 taataaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATTATCCAGGAGGTCTAGTATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTCTTTTAAGCTTTTTCTTTAATACAAAt                                                                                           >  2:308142/1‑141 (MQ=255)
                                                   ataaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTATTACtaataa                                                                                         >  2:27591/1‑141 (MQ=255)
                                                   ataaAACTAATGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGAAACAAAAAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATaa                                                                                         <  2:478544/141‑1 (MQ=255)
                                                        aCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGCTTATTTTAAGCTTTTTCTTTAATACAAATAACGttt                                                                                    >  2:461068/1‑141 (MQ=255)
                                                            ttGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCt                                                                                >  2:240966/1‑141 (MQ=255)
                                                            ttGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCt                                                                                >  2:496474/1‑141 (MQ=255)
                                                                tttATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGt                                                                            <  2:204660/141‑1 (MQ=255)
                                                                tttATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGt                                                                            >  1:171714/1‑141 (MQ=255)
                                                                 ttATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTTCTAAGtt                                                                           >  1:224954/1‑140 (MQ=255)
                                                                        tGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACtttt                                                                    <  2:406660/141‑1 (MQ=255)
                                                                        tGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACtttt                                                                    >  1:212247/1‑141 (MQ=255)
                                                                                aTACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCtt                                                            <  2:412749/141‑1 (MQ=255)
                                                                                   cGATATGTAAATGATAATCATCGAGGCGGTCTAGCATGTAACGTTCAAAAAAATACTATTACTTATCTAGCTTAAAGATCATTTTAAGCTTTTTCTTTAATACAAATAATGTTTTCCTAAGTGGACATTTTAATTCttttc                                                         >  2:120130/1‑140 (MQ=255)
                                                                                      tatGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTttatta                                                      >  1:11902/1‑141 (MQ=255)
                                                                                        tGTAAATGATAATCATCTAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTattaaa                                                    <  2:450912/141‑1 (MQ=255)
                                                                                        tGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTattaaa                                                    >  1:118456/1‑141 (MQ=255)
                                                                                          tAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTattaaatt                                                  <  1:472663/141‑1 (MQ=255)
                                                                                           aaaTGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTattaaatta                                                 >  2:229243/1‑141 (MQ=255)
                                                                                            aaTGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTattaaattaa                                                >  1:390358/1‑141 (MQ=255)
                                                                                              tGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATa                                              <  1:342444/141‑1 (MQ=255)
                                                                                                     catcCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCt                                       >  1:237026/1‑141 (MQ=255)
                                                                                                       tcCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGc                                     <  2:204484/141‑1 (MQ=255)
                                                                                                                        gTCTCGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAAt                    >  1:180369/1‑141 (MQ=255)
                                                                                                                         tctcGTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAAtt                   >  1:77911/1‑141 (MQ=255)
                                                                                                                             gTTCAAAAAAATAATTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTAc               >  1:88344/1‑141 (MQ=255)
                                                                                                                                  aaaaaaTACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAAtt          >  1:308599/1‑141 (MQ=255)
                                                                                                                                        tACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAATTTTGTCt    <  2:91903/141‑1 (MQ=255)
                                                                                                                                         aCTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAATTTTGTCta   >  1:218008/1‑141 (MQ=255)
                                                                                                                                          cTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAATTTTGTCtat  >  2:180580/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCTAGCTTAATGATTATTTTAAGCTTTTTCTTTAATACAAATAACGTTTTCCTAAGTGGACTTTTTAATTCTTTTATTAAATTAATACTTTTCTGCAGTGTTATTAACTCAATTGTACTAATTTTGTCTAT  >  CP000730/1909076‑1909354

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