Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 757,864 T→C 9.1% T380T (ACT→ACC USA300HOU_0707 → ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730757,8640TC 9.1% 141.3 / 6.8 57T380T (ACT→ACCUSA300HOU_0707ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base T (33/18);  new base C (3/2);  total (37/20)
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.84e-01

GATAGATTGGTATTAAATGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGTAGG  >  CP000730/757730‑758004
                                                                                                                                      |                                                                                                                                            
gATAGATTGGTATTAAATGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATTCCCGCTTGTACATCTTATTGCAGGAGTTGATCAGCCAACAATAGGTACTATAAGt                                                                                                                                        >  1:549937/1‑141 (MQ=255)
       tGGTATTAAATGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAAcc                                                                                                                                 <  2:294412/141‑1 (MQ=255)
        ggTATTAAATGATTTAACTTTGGAAATAGTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATGAGCCAGCAATAGGTACTATAAATAGAAAACa                                                                                                                                >  1:147625/1‑141 (MQ=255)
           aTTAAATGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCg                                                                                                                             <  2:138054/141‑1 (MQ=255)
               aaTGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGAt                                                                                                                         <  2:430114/141‑1 (MQ=255)
               aaTGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGAt                                                                                                                         <  1:92400/141‑1 (MQ=255)
                aTGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGAtt                                                                                                                        >  1:205409/1‑141 (MQ=255)
                   aTTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTGAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTaa                                                                                                                     >  2:16074/1‑141 (MQ=255)
                   aTTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTaa                                                                                                                     >  1:439666/1‑141 (MQ=255)
                   aTTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGGGATTTaa                                                                                                                     >  2:234672/1‑141 (MQ=255)
                       aaaTTTGGAAATATTTAAAGGTGATCAACTTGCACTTGTTGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCCAGCGTTTATCAGACAACAATAGGTACTATAAGTACAAACCCGCGTGATTTAAatat                                                                                                                 <  2:2328/141‑1 (MQ=255)
                        aaTTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAAtata                                                                                                                >  1:55306/1‑141 (MQ=255)
                          tttGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAgg                                                                                                              <  2:429716/141‑1 (MQ=255)
                            tGGAAATATGTGAAGGTGATCAAAGTGCTCTTGTAGGTCCAAGCGGGGCCGGTAACTCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCGACAATAGGGGCTATAAGTACAAACCAACGTGATTTAAGTATAGGaa                                                                                                            >  1:358278/1‑141 (MQ=255)
                             ggAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAGCCCAGCGTGATTTAAATATAAGAAt                                                                                                           >  2:419660/1‑141 (MQ=255)
                             ggAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACCATAAGTACAAACCAGCGTGATTTAAATATAGGAAt                                                                                                           >  2:318685/1‑141 (MQ=255)
                             ggAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACCATAAGTACAAACCAGCGTGATTTAAATATAGGAAt                                                                                                           >  2:318581/1‑141 (MQ=255)
                              gAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGGGGTTATCAGCCAACAATAGGGACTATAAGTACAAACCAGCGAGATTTAAAGATAGGAATa                                                                                                          >  1:197528/1‑141 (MQ=255)
                                        aaGGTGATCAAATTGCATTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGCAACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCATCGTGATTTAAATATAGGAATACTTAGTCAAc                                                                                                <  1:88350/141‑1 (MQ=255)
                                            tGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGcc                                                                                            <  1:329254/141‑1 (MQ=255)
                                             gATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTCAATCCACTTTGACACATCTTATTGCAGGTGTTTCTCAGCCAACAATAGGTACTATAAGTACAAACCGGCGCGGTTTAAATATAGGCATACTTAGTCAAAAGccc                                                                                           >  2:149625/1‑140 (MQ=255)
                                             gATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCa                                                                                           >  1:136844/1‑141 (MQ=255)
                                             gATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATAAGCCAACAATAGGTACTGTAAGTACAACCCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCa                                                                                           >  2:249830/1‑141 (MQ=255)
                                              aTCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTACATATAGGAATAATTAGTCAACAGCCaa                                                                                          >  1:154547/1‑140 (MQ=255)
                                              aTCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATAGTTAGTCAACAGCCat                                                                                          >  1:359338/1‑141 (MQ=255)
                                                        aCTTGAAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAg                                                                                <  2:404781/141‑1 (MQ=255)
                                                                gTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAAGAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGAGCTTCta                                                                        >  2:194465/1‑141 (MQ=255)
                                                                  ccAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCGACGGTAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCtata                                                                      <  2:539255/141‑1 (MQ=255)
                                                                  ccAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCtata                                                                      >  1:431040/1‑141 (MQ=255)
                                                                    aaGCGGGGCAGGTAAATCCACTTTGACACACCTTATTGCCGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATaaa                                                                    <  1:186653/141‑1 (MQ=255)
                                                                     aGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATaaaa                                                                   >  2:58940/1‑141 (MQ=255)
                                                                        ggggCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAagag                                                                >  1:46928/1‑141 (MQ=255)
                                                                         gggCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAgcga                                                               >  1:383886/1‑141 (MQ=255)
                                                                            cAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTAAAAAAGAGAata                                                            >  1:318911/1‑141 (MQ=255)
                                                                              ggTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACGCAAATGCGTGATTTAAATATATGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATAtt                                                          <  1:268498/141‑1 (MQ=255)
                                                                                  aaTCGACTTTGAGACATCTTATTGCAGGTGTTTATCAGCCAACAATCGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCATCAGCCGTATATTTTCAGTGCTTCTATAGAAGAGAATATTACGa                                                      >  1:395715/1‑141 (MQ=255)
                                                                                              cacaTCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGata                                          >  1:306014/1‑141 (MQ=255)
                                                                                               acaTCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTAAAAACGAGAATATTACGATGTTTAAAGatat                                         >  2:45912/1‑141 (MQ=255)
                                                                                                 aTCTTATTGCAAGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAg                                       <  2:519965/141‑1 (MQ=255)
                                                                                                   cTTATTGCAGGTGTTTATCAGCCAACAATAGGTACCATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGaa                                     <  1:318581/141‑1 (MQ=255)
                                                                                                   cTTATTGCAGGTGTTTATCAGCCAACAATAGGTACCATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGaa                                     <  1:318685/141‑1 (MQ=255)
                                                                                                    ttATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGaaa                                    >  2:209469/1‑141 (MQ=255)
                                                                                                      aTTGCAGGTGTTTATCAGCCAACAATAGGTACCATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGTTATAGaaaaa                                  >  2:52841/1‑140 (MQ=255)
                                                                                                        tGCAGGTGTTTATCAGCCAACCATAGGTACTAAAAGTACAAACCAGCGAGATAAAAAGATAGGAATACTTAGTAAACAGCAAAATATTTTCAGTGCTTCCATAAAAGAGAAGATTACGACGTTTAAAGATCTACAAaataa                                >  2:194462/1‑141 (MQ=255)
                                                                                                            ggTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACt                            <  1:399180/141‑1 (MQ=255)
                                                                                                               gtTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGAGATAGAAAATAATACGAtt                         >  2:243463/1‑141 (MQ=255)
                                                                                                                     cAGCCAACAATAGCTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTgaagaa                   <  2:549937/141‑1 (MQ=255)
                                                                                                                      aGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGaagaag                  <  2:348216/141‑1 (MQ=255)
                                                                                                                         caacaaTAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGc               >  1:10582/1‑141 (MQ=255)
                                                                                                                          aacaaTAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATGGGAAATAATACTATTGAAGAAGTGca              >  1:154117/1‑140 (MQ=255)
                                                                                                                          aacaaTAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAATAATGCTATTGAAGAAGTGCTa             >  1:154171/1‑141 (MQ=255)
                                                                                                                             aaTAGGTACTATAAGTACACACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGa           >  1:335675/1‑141 (MQ=255)
                                                                                                                               tAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACg         <  1:398910/141‑1 (MQ=255)
                                                                                                                               tAGGTACAATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACg         >  1:350435/1‑141 (MQ=255)
                                                                                                                                   tACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGt     >  1:457108/1‑141 (MQ=255)
                                                                                                                                   tACTATAAGGACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGt     <  2:54769/141‑1 (MQ=255)
                                                                                                                                     cTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATCGAAAATAATACTATTGAAGAAGTGCTAGACGAAGTAg   >  1:290597/1‑141 (MQ=255)
                                                                                                                                      tataAGTACAAACCAGCTTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGTAgg  <  1:525167/141‑1 (MQ=255)
                                                                                                                                      tataAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGTAgg  <  2:295763/141‑1 (MQ=255)
                                                                                                                                      tataAGTACAAACCAGCGTGATTTAAAAATAGGAATACTTAGGCAACAGCCATATATTTTCAGTGCTTCTCTAAACGAGAATCATCCAAAGTTTAAAGCAGTACAGAATAATACGGTTGAAGAACAGATAGCGGAAGTAgg  >  2:246758/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                            
GATAGATTGGTATTAAATGATTTAAATTTGGAAATATTTAAAGGTGATCAAATTGCACTTGTAGGTCCAAGCGGGGCAGGTAAATCCACTTTGACACATCTTATTGCAGGTGTTTATCAGCCAACAATAGGTACTATAAGTACAAACCAGCGTGATTTAAATATAGGAATACTTAGTCAACAGCCATATATTTTCAGTGCTTCTATAAAAGAGAATATTACGATGTTTAAAGATATAGAAAATAATACTATTGAAGAAGTGCTAGACGAAGTAGG  >  CP000730/757730‑758004

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: