Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,151,090 A→G 9.9% K66K (AAA→AAG USA300HOU_1076 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,151,0900AG 9.9% 104.1 / 3.9 41K66K (AAA→AAGUSA300HOU_1076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base G (2/2);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGTTTTTG  >  CP000730/1150966‑1151230
                                                                                                                             |                                                                                                                                            
cAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATCAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAATAATAAAACTagaagaaa                                                                                                                              >  2:445675/1‑140 (MQ=255)
       caTACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGacaacca                                                                                                                       >  2:380413/1‑138 (MQ=255)
       caTACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCg                                                                                                                       >  2:444432/1‑141 (MQ=255)
            aatatGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGtt                                                                                                                  <  1:552659/140‑1 (MQ=255)
                    gCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAaacaa                                                                                                          >  2:264782/1‑141 (MQ=255)
                      acaTTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACGACTAGAAGAAGAAAATCGACGTTGGAAaaaaca                                                                                                        >  2:162959/1‑141 (MQ=255)
                            aGTTGATGATAAAATAAAAGAATTAGGGCATAAAGCAGCAGGTTTAGATAATTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAAACGACGTTTAAAACAACAACTTca                                                                                                  >  2:149926/1‑141 (MQ=255)
                                     aaaaagaaaagAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTgca                                                                                         <  1:203261/140‑1 (MQ=255)
                                      aaaataaaaGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTgcag                                                                                        <  1:614097/141‑1 (MQ=255)
                                            agaggatgagagTATAAAGCAGCAGATTTAGATGCTTGAAGAGAGGCAATACTAACTACTGTGAATATTATGCATAAAAAAGTACTACTAAAAGAAGATAATCGAGGTTTGAAACAACAAATTCACAAATTGCAGCAGCGt                                                                                  >  2:341806/12‑141 (MQ=255)
                                               gAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTACGCATGAAAAAGTAATAGTAGAAGAAGAAATACGACGTTTGAAACAACAAATTCACCAATTGCAGCACGGTGAg                                                                               >  2:658034/1‑141 (MQ=255)
                                                      ggTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATaaa                                                                        >  2:332147/1‑141 (MQ=255)
                                                        tataAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATg                                                                      >  2:391695/1‑141 (MQ=255)
                                                        tataAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATAAAAACACGCTCGGGGATGATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATg                                                                      <  1:379919/141‑1 (MQ=255)
                                                              gcagcaGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCAtt                                                                <  1:367631/141‑1 (MQ=255)
                                                              gcagcaGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCAtt                                                                >  2:415062/1‑141 (MQ=255)
                                                                 gcagGTTTAGATACTTCAAGAACAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGAt                                                             <  1:424246/141‑1 (MQ=255)
                                                                   agGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGAttt                                                           >  1:436986/1‑141 (MQ=255)
                                                                   agGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATAGTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAGATTCACAAATTGGAGCAGCGTGACGAATAACTGGTCAttatatt                                                           >  2:514225/1‑136 (MQ=255)
                                                                           aTACTTCAAGAAAAGCACTACTAACAACGGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  2:472903/1‑141 (MQ=255)
                                                                           aTACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  1:412471/1‑141 (MQ=255)
                                                                           aTACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  1:107134/1‑141 (MQ=255)
                                                                           aTACTTCAAGAAAAGCAATACTAACTACTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   <  1:257553/141‑1 (MQ=255)
                                                                           aTACTGCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   <  1:546299/141‑1 (MQ=255)
                                                                            tACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTGTCATAAt                                                  >  1:598472/1‑141 (MQ=255)
                                                                            tACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAAt                                                  >  1:365474/1‑141 (MQ=255)
                                                                              cTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCa                                                >  2:243268/1‑141 (MQ=255)
                                                                                tCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCAtt                                              <  1:444432/141‑1 (MQ=255)
                                                                                tCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCAtt                                              <  1:380413/141‑1 (MQ=255)
                                                                                   aGACAAGCAATACTAACTGCTGTCAATATTATGCATGAAAACGTACTAATAGAAGCAGAACATAGACGTTTGAAACACCAAATTCTCCAATTGCGACAGGGTGAGCAATAAATGGCCATTGATTTTATCATAATCAtcttt                                           >  2:187680/1‑140 (MQ=255)
                                                                                     aaaaGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGACGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGGCATTGATTTTATCATAATCATTTTCtt                                         >  1:26528/1‑141 (MQ=255)
                                                                                      aaaGCAATACTAACTGCTGGGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGTTCATTGATTTTATCATAATCATTTTCttt                                        >  2:376588/1‑141 (MQ=255)
                                                                                       aaGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTtg                                       >  2:555746/1‑141 (MQ=255)
                                                                                        aGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTtgt                                      >  1:308918/1‑141 (MQ=255)
                                                                                          cAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTgtgt                                    >  1:343775/1‑141 (MQ=255)
                                                                                          cAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTgtgt                                    >  2:574434/1‑141 (MQ=255)
                                                                                             tACTAACTGCTGTGAATATTGTGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTAtt                                 >  1:468234/1‑141 (MQ=255)
                                                                                                             atatTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGAt                  <  1:243268/141‑1 (MQ=255)
                                                                                                              tattatGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGAtt                 >  2:254061/1‑141 (MQ=255)
                                                                                                                  atGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGa             >  2:526100/1‑141 (MQ=255)
                                                                                                                   tGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGAc            <  1:391695/141‑1 (MQ=255)
                                                                                                                           aaaGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGtttt    <  1:503388/141‑1 (MQ=255)
                                                                                                                            aaGTACTACTAGAAGAAGAAAATCGACGTTTGCAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTACCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGttttt   <  1:445675/141‑1 (MQ=255)
                                                                                                                             aGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCGTACGCGGTTTTTg  >  2:398802/1‑141 (MQ=255)
                                                                                                                             |                                                                                                                                            
CAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAACTGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGTTTTTG  >  CP000730/1150966‑1151230

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: