Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 41.9% intergenic (+89/+19) USA300HOU_0025 → / ← USA300HOU_0026 hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8580TA41.9% 124.0 / 109.2 96intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (9/31);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 5.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATT  >  CP000730/35720‑35971
                                                                                                                                          |                                                                                                                  
aaGTATTTAATGTAGAAGCTATTGAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                                  >  2:548819/1‑141 (MQ=255)
aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCAtttt                                                                                                                  >  2:226169/1‑140 (MQ=255)
aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                                  >  1:534413/1‑141 (MQ=255)
aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACAACTAATTTa                                                                                                                  >  1:431520/1‑141 (MQ=255)
  gTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATc                                                                                                                >  1:187510/1‑141 (MQ=255)
    aTTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCAt                                                                                                              <  1:414221/141‑1 (MQ=255)
     tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAGAAATGATTACTTGAACAAGAGCTAACATAAAATTGTGATCTAATAAAAATAGAAGGTTGTGTTGCAACGTAAAAAAATATAGATAACCAGTAATTTATTATg                                                                                                             >  2:617607/1‑141 (MQ=255)
     tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                                             >  1:374363/1‑141 (MQ=255)
       tAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTAGCATGTc                                                                                                           <  2:424938/141‑1 (MQ=255)
        aaTGTAGAAGCTATTGAGTATAAAGCTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCa                                                                                                          >  1:386021/1‑141 (MQ=255)
        aaTGGAGAAGCTAATCAGTATAAAGGTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTAATAAAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCa                                                                                                          <  2:417546/141‑1 (MQ=37)
            tAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAgtgt                                                                                                      <  1:161337/141‑1 (MQ=255)
              gAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGCGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTCTTGCAATGTAAAAAAATATGGCTAACCACTAATGTATCATGACAGCGTTc                                                                                                    >  1:78298/1‑141 (MQ=255)
                aGCTACTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCGGTTGCAAAGTAAAAAAATATAGCTAACCTCTAATTTATCATGTCAGTGTTCGc                                                                                                  <  1:145417/141‑1 (MQ=255)
                 gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                                                 >  1:592981/1‑141 (MQ=255)
                    aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCAAAAATAAAATTGTGATCTAATAAAAATAGAAGGTACTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATATAAGTGTTCACATaa                                                                                              >  2:713768/1‑141 (MQ=37)
                       cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACtt                                                                                           >  2:403176/1‑141 (MQ=255)
                        aGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTg                                                                                          >  1:631679/1‑141 (MQ=255)
                         gTATAAAGTTTTTGGCGAAAAATGATTAGTTGAACAAGAGCTAATAGAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                         <  2:113515/141‑1 (MQ=255)
                                 tttttGGCGAAAAATGATTACTTGAACCAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTACTTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                                 >  2:89857/1‑141 (MQ=255)
                                        cGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatgatcac                                                                          >  2:377117/1‑136 (MQ=255)
                                           aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTc                                                                      <  1:404544/141‑1 (MQ=255)
                                             aaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGt                                                                    <  2:744828/141‑1 (MQ=255)
                                              aTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTg                                                                   <  1:125692/141‑1 (MQ=255)
                                               tGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTgg                                                                  >  1:238577/1‑141 (MQ=255)
                                                gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                                 >  1:117231/1‑141 (MQ=255)
                                                gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                                 >  1:367499/1‑141 (MQ=255)
                                                gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATACTAATTTCTTGGc                                                                 >  1:285700/1‑141 (MQ=255)
                                                 aTTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCa                                                                >  1:147266/1‑141 (MQ=255)
                                                     cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                            <  1:283214/141‑1 (MQ=255)
                                                     cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                            <  1:161676/141‑1 (MQ=255)
                                                     attaaaaaaGAGCTAAAATAAAATTGTTATCTAAAAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                            <  2:534380/134‑1 (MQ=255)
                                                             agagCTAAAATAAAATTGTGGTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                    >  2:14839/1‑141 (MQ=255)
                                                              gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                   <  2:210447/141‑1 (MQ=255)
                                                              gagCTAAAATAAAACTGGGATCTAATAAAAATCGAAGGTTCTGTTGCAACGACAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACGTGCTAGCATGATGCTACTTACGTGGccgtgcgagactccg                                                   >  2:172539/1‑127 (MQ=17)
                                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAGAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAAGTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                  >  2:336004/1‑141 (MQ=255)
                                                                gCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTGGCATGCTGCCCAATCTCTGGCATGGCGAAAATCCGTa                                                 <  1:184987/141‑1 (MQ=18)
                                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGCTc                                             >  1:170029/1‑141 (MQ=255)
                                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                             >  1:94181/1‑141 (MQ=255)
                                                                      acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCCGaa                                         <  1:245252/127‑1 (MQ=11)
                                                                       taaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATACCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGa                                          >  1:146200/1‑141 (MQ=255)
                                                                        aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                       <  2:325815/129‑1 (MQ=11)
                                                                         aaatTGTTATCTAATAAAAATAGAAGTCGCTGTGGCAAAGTAACAAAATATAGCTAACCACTAATTTATCATGTCAGTGTCCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                        <  2:85341/141‑1 (MQ=14)
                                                                         aaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                        <  1:669591/141‑1 (MQ=255)
                                                                         aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                      <  1:741271/130‑1 (MQ=11)
                                                                           atTGTGATCTAATAAAAATAGAAGGGTTTGTTGAAAAGAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                      <  2:273787/141‑1 (MQ=255)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  1:102241/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  1:640564/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  1:452464/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  1:104734/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  2:461944/132‑1 (MQ=11)
                                                                           agagtagttgaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    <  2:374421/131‑1 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                   <  2:240880/133‑1 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                   <  1:98840/133‑1 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGATAATCACTAATATATCCTGTCAGTGTTCGTTTAACTTACTAGCATACTGCTACTTTCGGGGCATGGCGAAAATCCGTAGATATGAAGAAAcc                                   >  2:747799/9‑141 (MQ=2)
                                                                             tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                    >  2:306295/1‑141 (MQ=255)
                                                                              gtgATCTAATACAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                   <  1:170459/141‑1 (MQ=255)
                                                                              gtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                   <  1:359196/141‑1 (MQ=255)
                                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                 >  2:529666/7‑141 (MQ=11)
                                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGATAATCACTAATATATAATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                 >  1:108390/7‑141 (MQ=11)
                                                                              gtattttaataaAAACGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCGAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                 <  1:23048/135‑1 (MQ=11)
                                                                              gtattttaataaAAACGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACATGCTAGCATGATGTTAATTTCGTGGCATGGCGAAAATTCGTAGATCTGAAGAGACCTg                                 <  1:466066/135‑1 (MQ=11)
                                                                               tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAACATCCGTAGATCTGAAGAGACCTGc                                >  2:329494/6‑141 (MQ=11)
                                                                               tattttaataaAAATGATGGTTCTGCAGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                <  1:587785/136‑1 (MQ=11)
                                                                               gattttaataaAAACGATGGGTCTTTTGCAAAGTAAAAAACATAGCTAATCACGAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                <  1:347449/136‑1 (MQ=11)
                                                                                gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAACTCAGTAGATCTAAAGAGACCTg                                 >  2:720003/1‑141 (MQ=255)
                                                                                gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGATGAGACCTg                                 >  2:525534/1‑141 (MQ=255)
                                                                                gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                 >  1:634764/1‑141 (MQ=255)
                                                                                attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAATTCCGTAGATCTGAAGAGACCTGCg                               <  2:589830/137‑1 (MQ=11)
                                                                                attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                               <  1:417923/137‑1 (MQ=11)
                                                                                 ttttaataaACATGATGGTTATGTTGCAGAGTAAAAAATATAGCTAATCACTAATATATAATGTCAGTGTGCGTTTAACTTGCTAGCATTATGCTATTTTCGTGGCATGGCGAAAATCAGTAGATCTGAAGATACCTGCgg                              >  2:238573/4‑141 (MQ=11)
                                                                                 ttttaataaAAATGATTGTTCTGTTGCAAATTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              <  2:705854/138‑1 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              <  2:54549/138‑1 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              >  2:326033/4‑141 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              <  1:71200/138‑1 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              <  1:298128/138‑1 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCGAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                              >  1:748206/4‑141 (MQ=11)
                                                                                 aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                <  1:130404/141‑1 (MQ=255)
                                                                                  ttaaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCCCTCATACATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                             <  1:260537/138‑1 (MQ=11)
                                                                                  tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                               <  1:19658/141‑1 (MQ=255)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCTGTAGATCTGAAGAGACCTGCGGtt                            <  1:717824/140‑1 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                            <  1:399858/140‑1 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                            <  2:41180/140‑1 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCAGtt                            <  1:551527/140‑1 (MQ=11)
                                                                                     aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAGCCCGTAGATCTGAAGAGACCTGCGGTTCt                          >  2:77445/1‑141 (MQ=11)
                                                                                     aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                          <  2:541797/141‑1 (MQ=11)
                                                                                     aataaAAATGATGGTTCTGTCGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                          <  2:578283/141‑1 (MQ=11)
                                                                                     aataaAAAGAGAAGGTTATGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                            <  2:224697/141‑1 (MQ=255)
                                                                                      ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                         <  1:581494/141‑1 (MQ=11)
                                                                                      ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                         <  2:624017/141‑1 (MQ=11)
                                                                                           aaTAGAAGGCTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                      <  1:239082/141‑1 (MQ=21)
                                                                                            aTAGACGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACGTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTta                     <  2:656717/141‑1 (MQ=17)
                                                                                             tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat                    <  1:165874/141‑1 (MQ=32)
                                                                                              aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTATCGTAACGTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                   <  1:557842/141‑1 (MQ=255)
                                                                                                aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCCGAAGAGACCTGCGGTTCTTTTtatata                 <  1:219886/141‑1 (MQ=11)
                                                                                                    ttcttgtGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCg            <  1:218007/135‑1 (MQ=1)
                                                                                                               aagtaacaaaTATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACAtt  >  2:90471/8‑141 (MQ=1)
                                                                                                                                          |                                                                                                                  
AAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATT  >  CP000730/35720‑35971

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: