Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,850 C→T 40.2% intergenic (+81/+27) 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,8500CT40.2% 101.2 / 69.9 86intergenic (+81/+27)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (7/27);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 5.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTC  >  CP000730/35714‑35972
                                                                                                                                        |                                                                                                                          
ttGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACt                                                                                                                        <  1:317586/141‑1 (MQ=255)
      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:624081/1‑141 (MQ=255)
                                                                     agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                   >  2:289982/1‑141 (MQ=255)
                                                                     agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                   >  1:418726/1‑141 (MQ=255)
                                                                     agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                   >  2:289826/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:675577/15‑141 (MQ=11)
                                                                            acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                          >  1:440208/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                                    >  1:766973/9‑141 (MQ=11)
                                                                                  gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                    >  2:402460/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:642688/137‑1 (MQ=11)
                                                                                      attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                <  1:659376/137‑1 (MQ=11)
                                                                                      attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                <  2:211134/137‑1 (MQ=11)
                                                                                       ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                               <  2:36853/138‑1 (MQ=11)
                                                                                       ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                               >  2:524611/4‑141 (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:376605/140‑1 (MQ=11)
                                                                                         ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                             >  1:36603/2‑141 (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)
                                                                                                                                        |                                                                                                                          
TTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTC  >  CP000730/35714‑35972

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: