Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 44.1% 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,8580TA44.1% 96.0 / 95.1 86intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base A (8/30);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 2.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35720‑35988
                                                                                                                                          |                                                                                                                                  
aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCAGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                                                  >  1:567609/1‑141 (MQ=255)
 aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAACATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAGTAATTTAt                                                                                                                                 <  2:570793/141‑1 (MQ=255)
 aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAGGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTCAt                                                                                                                                 >  2:735321/1‑141 (MQ=255)
       tAATGTAGAATTTATCAAGTATAGCGTTGTGTGCTAAAAAGGATTCCTTGAACAAGTCCTAAAATAAAATTGTAATCTAATACCCATATAAGGTTCTAAAAAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTc                                                                                                                           <  2:313693/141‑1 (MQ=12)
       tAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAATAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTGACCACTAATTTATCaagac                                                                                                                           >  1:221268/1‑137 (MQ=255)
           gTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCtgtg                                                                                                                       >  1:419738/1‑141 (MQ=255)
               aaGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCg                                                                                                                   <  1:775174/141‑1 (MQ=255)
                aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                                                                  >  2:52105/1‑141 (MQ=255)
                 gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTACTAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                                                                 >  1:463112/1‑141 (MQ=255)
                      tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                                                            >  1:691791/1‑141 (MQ=255)
                                 tttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAACGAAACAAAATATAGCTAACCACTAATTGATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                                                 <  1:145606/141‑1 (MQ=255)
                                 tttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTGGCatga                                                                                                 >  1:314522/1‑141 (MQ=255)
                                 gttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                                                 <  1:171141/140‑1 (MQ=255)
                                   tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAACTATAGCTAACCACTAATTTATCATGTGTGTGGTCGATTAACTTGCTAGCatgctg                                                                                               >  1:21021/1‑141 (MQ=25)
                                          aaaaaTGATTACTTTAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAAttt                                                                                        >  1:426570/1‑141 (MQ=255)
                                                gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                                                  <  1:390212/141‑1 (MQ=255)
                                                gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                                                  >  1:287414/1‑141 (MQ=255)
                                                     cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                                             <  1:646081/141‑1 (MQ=255)
                                                     cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGGAAAAAAATATAGCTTACCACTGATTTATCATGTCAGTGTTCGCGTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                                             <  1:57565/141‑1 (MQ=255)
                                                              gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTGATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                                    >  1:522163/1‑141 (MQ=255)
                                                              gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGCTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCACGGCGAGAATCCg                                                                    >  2:126760/1‑141 (MQ=255)
                                                              gagCTAAAAGAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                                    <  2:54439/141‑1 (MQ=255)
                                                              gaacTAAAATAAAATTGTGATATAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                                    <  2:266392/138‑1 (MQ=255)
                                                               agCTAAAATAAAATTGTGATGTAATAAAAATAGAAGGTTCTGTTCCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGTATGATGCTAATTTCGTGGTATGGCGAAAATCCGt                                                                   >  1:469132/1‑141 (MQ=255)
                                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                   >  2:289826/1‑141 (MQ=255)
                                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                   >  1:418726/1‑141 (MQ=255)
                                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                   >  2:624081/1‑141 (MQ=255)
                                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                   >  2:289982/1‑141 (MQ=255)
                                                                 gtaaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTTCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAg                                                                 <  2:36323/140‑1 (MQ=255)
                                                                 caaaaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTACGTGGCATGGCGAAAATCCGTAg                                                                 <  2:553619/139‑1 (MQ=255)
                                                                 ataaaaaaaaaTTGTGATCTAATAAAAATAGAAGGCTCTGTTGCAAAGTAAAAAAAGATAGCTAACCACTAATTTATCATGGCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAg                                                                 <  2:532436/140‑1 (MQ=255)
                                                                      acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                          >  1:440208/15‑141 (MQ=11)
                                                                      acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                          >  1:675577/15‑141 (MQ=11)
                                                                         aaatTGTGATCTAATAAAAATAAAAAGTTCTTTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  2:345091/141‑1 (MQ=255)
                                                                         aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAAACACTAATATATCATGACAGAGTTCGTTTCACTTGCCAGCATAGTGCTACTTTCGTGGCATGGCGAAAAACGGTAGATCTCAaggg                                                       >  2:532023/12‑139 (MQ=2)
                                                                          aatTGGTATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                                        <  2:143856/141‑1 (MQ=255)
                                                                           atagtagtttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAACATATCATGCCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCTTAGATCTGAAGAGAc                                                     <  2:488477/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:754523/132‑1 (MQ=11)
                                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:263126/132‑1 (MQ=11)
                                                                           agagtaattcaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATTTAGATTACCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:682058/131‑1 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                    >  2:402460/9‑141 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                    >  1:766973/9‑141 (MQ=11)
                                                                            gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                    >  2:402488/9‑141 (MQ=11)
                                                                             agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                   <  2:672182/134‑1 (MQ=11)
                                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  2:72156/135‑1 (MQ=11)
                                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  2:185896/135‑1 (MQ=11)
                                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAACCCGTAGATCTGAAGAGACCTg                                                  <  2:94221/135‑1 (MQ=11)
                                                                               tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  2:10339/136‑1 (MQ=11)
                                                                               tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCGAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  2:155131/136‑1 (MQ=11)
                                                                               tattttaataaAAATGAGGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:679578/136‑1 (MQ=11)
                                                                               ggATCGAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                   <  1:232169/140‑1 (MQ=255)
                                                                                attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  1:659376/137‑1 (MQ=11)
                                                                                attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  1:642688/137‑1 (MQ=11)
                                                                                attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  2:211134/137‑1 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               >  2:524611/4‑141 (MQ=11)
                                                                                 ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               <  2:36853/138‑1 (MQ=11)
                                                                                 aTCTAATAAAAGTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  2:114978/141‑1 (MQ=255)
                                                                                 aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:378033/141‑1 (MQ=255)
                                                                                  tttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATACCCAATCCCAAATACTTCATGTAAGTCTTTGTTTAACTTGCTAGCATTATGCTAATTTAGTGGCATGGCGAAAATCCGTAGACCTGAAGAGACCTGCGGt                                              <  2:193526/139‑1 (MQ=2)
                                                                                  tttaataaAAATGATGTTTCTGTTGCAAAGTAACAAATATAGCTAATCCCCATTATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAACTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  1:349562/139‑1 (MQ=11)
                                                                                  tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:646707/139‑1 (MQ=11)
                                                                                  tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:673339/139‑1 (MQ=11)
                                                                                  ttgaataaAAATGATGCTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATAGCATGTCAGTGTTCGTTTAACTTGCTAGCATGTTGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  1:245161/138‑1 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                             >  1:36603/2‑141 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                             <  1:376605/140‑1 (MQ=11)
                                                                                   ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAAGCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                             <  2:682348/140‑1 (MQ=11)
                                                                                   cTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAAAGATATAGCTAACCACTAATTTATCAGTGCAGTGGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGAgg                                               >  2:261106/1‑141 (MQ=14)
                                                                                   cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTACTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               >  1:287213/1‑141 (MQ=255)
                                                                                    taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              >  2:256714/1‑141 (MQ=40)
                                                                                    taataaAAATGATGGTTGTGTTGCAAAGTCAAAAATCTAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                            <  1:3529/141‑1 (MQ=11)
                                                                                    taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                            <  1:471129/141‑1 (MQ=11)
                                                                                    taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                            <  1:587593/141‑1 (MQ=11)
                                                                                    taataaAAATGATGGTTATGCTGCAAAGTAAAAAATATAGATAATCACTAAGATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCGGAAGAGACCTGCGGTTc                                            <  1:107998/141‑1 (MQ=11)
                                                                                     aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGACCTGAAGGGCCCTGCGGtt                                             >  2:142565/1‑141 (MQ=255)
                                                                                     aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCACTTAACTTGCTAGCATGATGCTAATTTCGTGGCAAGGAGAGAATCAGTAGTTGTTAAGAGACGTGCGAtt                                             >  2:36613/1‑141 (MQ=14)
                                                                                     aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                           <  1:592031/141‑1 (MQ=11)
                                                                                       taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGCAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                           >  1:108364/1‑141 (MQ=38)
                                                                                          aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                        <  1:541634/141‑1 (MQ=35)
                                                                                          aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCGTAACTTGCTAGCATGATGCTAATTTCGTGGCATGTCGAAATTCCGTGTATCTGAAGAGACCTGCGGTTCtttt                                        <  1:735321/141‑1 (MQ=11)
                                                                                           aaTAGAAGGTTATGATGCAAAGAACCAAAATATAGCTAACTACTAATTTATCATGTCAGGGTTCGCTTAATTTGCGAGCATGATTCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                       >  2:534753/1‑141 (MQ=11)
                                                                                              aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                    >  1:687715/1‑141 (MQ=32)
                                                                                                   gtccTGTTGCAAAGTAAAAAATAGAGCTAATCACTAATAGTTCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGc                              <  2:725349/138‑1 (MQ=1)
                                                                                                        ggtgaaaaGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCGGCGGTTCTTTTTATATAGAGCGTaa                          <  2:57164/136‑1 (MQ=0)
                                                                                                                atataaaaaTATAGGCAACCACTATCTTATCATGTATGTGTCCGCTTAACTCTCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTc                  <  2:142673/137‑1 (MQ=1)
                                                                                                                agtaacaaaTAAAGCTAATCACTAAGATATCATGTCAGTGTTCGTTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACactc                  <  1:77199/135‑4 (MQ=1)
                                                                                                                                aatcACTAATATATCATGGCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa  <  1:5763/138‑1 (MQ=1)
                                                                                                                                          |                                                                                                                                  
AAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35720‑35988

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: