Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,850 C→T 42.3% 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,8500CT42.3% 100.2 / 80.0 89intergenic (+81/+27)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base T (13/25);  total (34/55)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATAC  >  CP000730/35710‑35977
                                                                                                                                            |                                                                                                                               
ctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAAc                                                                                                                                 <  1:472113/141‑1 (MQ=255)
ctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAAc                                                                                                                                 <  1:34909/141‑1 (MQ=255)
ctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAAc                                                                                                                                 <  1:748248/141‑1 (MQ=255)
ctcATTGACAAAGTATTTAATGTAGAAGCCATTCAGTATAAATTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAAc                                                                                                                                 <  1:748347/141‑1 (MQ=255)
       aCAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAt                                                                                                                          <  2:632199/141‑1 (MQ=255)
         aaaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAttt                                                                                                                        >  1:247519/1‑141 (MQ=255)
           aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATGAAAATGGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTAt                                                                                                                      >  2:801891/1‑141 (MQ=255)
            gTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACCAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAGAAAATATAGCTAACCACTAATTTATc                                                                                                                     >  1:514774/1‑141 (MQ=255)
                  aaTGTAGAAGCTATTCAGTATAAAGATTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCa                                                                                                               >  1:419807/1‑141 (MQ=255)
                   aTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                                                              <  1:519546/141‑1 (MQ=255)
                   aTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                                                              <  1:832205/141‑1 (MQ=255)
                         aaGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCg                                                                                                        >  1:718788/1‑141 (MQ=255)
                          aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAAAAAAATTGTGATCTCATAAAAATAGAAGGTTCTGAAACAAAGAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                                                       <  2:738458/141‑1 (MQ=38)
                            cTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACGAATTTATCATGTCAGTGTTCGCtt                                                                                                     <  1:610640/141‑1 (MQ=255)
                                 cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACtt                                                                                                <  1:358008/141‑1 (MQ=255)
                                   gTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                              >  2:147429/1‑141 (MQ=255)
                                     ataAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTa                                                                                            >  1:220302/1‑141 (MQ=255)
                                        aaGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCa                                                                                         >  1:22430/1‑141 (MQ=255)
                                             tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                    <  1:139089/141‑1 (MQ=255)
                                               tGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCt                                                                                  >  1:655806/1‑141 (MQ=255)
                                                 gCGAAGAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                                                <  2:110178/141‑1 (MQ=255)
                                                 gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                                                >  2:270981/1‑141 (MQ=255)
                                                 gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                                                >  1:170303/1‑141 (MQ=255)
                                                        aTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTg                                                                         <  2:247519/141‑1 (MQ=255)
                                                         tGATTACTTGAACAACAGCGGCGCAAAGATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTgg                                                                        >  2:181572/1‑141 (MQ=17)
                                                           aTTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCa                                                                      <  1:671244/141‑1 (MQ=255)
                                                            ttACTTGAACAAGAGCTAAAAGAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAGGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCag                                                                     >  1:436035/1‑140 (MQ=255)
                                                               cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAATTAAAAAAATATAGCTAACCACTAATTTATCAGGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                                  <  2:72996/141‑1 (MQ=255)
                                                               cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTGATCATGTCAGTGTTCGCTTAACTTGCTAGCGTGATGCTAATTTCGTGGCATGGc                                                                  <  2:562144/141‑1 (MQ=255)
                                                                      aaGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                                           <  1:258160/141‑1 (MQ=255)
                                                                        gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTTAGTGTTCGCTTAACTTGCTCGCATGATGATAATTTAGTGGCATGGCGAAAATCCg                                                         >  1:507400/1‑141 (MQ=255)
                                                                        gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                         >  1:551358/1‑141 (MQ=255)
                                                                        aagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                         <  2:710303/140‑1 (MQ=255)
                                                                             aaaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAt                                                    >  2:95281/1‑141 (MQ=255)
                                                                             aaaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAt                                                    >  1:821036/1‑141 (MQ=255)
                                                                               aataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCt                                                  >  1:632977/1‑141 (MQ=255)
                                                                                acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAGATCCGTAGATCTGaa                                               >  2:345561/15‑141 (MQ=11)
                                                                                 caaaagagtattCTACTAAAAATGATGGTTCTGTTGCGAAGTAAAAAATATAGCTAATCACTCATATAGCAGGTCAGTGTTCGTTTAACTTGCGAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                              <  2:517678/130‑1 (MQ=11)
                                                                                  aaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAAGTTGCTACAATCATGCTAATTTCGTGGCATGCCGAAAATCCGTAGATCTGaa                                               >  2:11244/1‑141 (MQ=17)
                                                                                  aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                             <  1:92811/129‑1 (MQ=11)
                                                                                   aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATACTAATTTCTTGGCAGCGTCAAAATCCGTAGATCTGAagag                                            <  1:616312/130‑1 (MQ=2)
                                                                                    aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                           <  1:833569/131‑1 (MQ=11)
                                                                                    aagagtattttaataaAAATGATGGTTCTGTCGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                           <  1:445313/131‑1 (MQ=11)
                                                                                     agggtattgtaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGGGTTCGTTTAACTTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  2:193729/132‑1 (MQ=11)
                                                                                     agcgtattttaataaAAATGATGGTTCTGTTGCAAAGAAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  2:706920/132‑1 (MQ=11)
                                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  1:778908/132‑1 (MQ=11)
                                                                                      gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                         >  1:665357/9‑141 (MQ=11)
                                                                                      gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                         <  2:537940/133‑1 (MQ=11)
                                                                                       tgtgCTCTGATAAAAATAGAAGGTTCTGTTACAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  2:839706/141‑1 (MQ=255)
                                                                                       tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  2:222094/141‑1 (MQ=255)
                                                                                       tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACGTTCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  2:760745/141‑1 (MQ=255)
                                                                                       tgtgATATAGTAAAAATAGCCGGTTATGAGGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                          <  1:393895/141‑1 (MQ=11)
                                                                                       agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                        <  1:413813/134‑1 (MQ=11)
                                                                                       agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                        <  1:375906/134‑1 (MQ=11)
                                                                                       agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                        <  1:413635/134‑1 (MQ=11)
                                                                                       agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATCCCGAAAATCCGTAGATCTGAAGAGACCt                                        <  2:464790/134‑1 (MQ=11)
                                                                                        gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                       <  1:733475/135‑1 (MQ=11)
                                                                                        gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                       >  1:496512/7‑141 (MQ=11)
                                                                                        gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                       <  1:698599/135‑1 (MQ=11)
                                                                                         tgATATAATAAAAATAGAAGGTTCTGTGGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                        <  2:463726/141‑1 (MQ=255)
                                                                                         tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATTGCGAAAATCCGTAGATCTGAAGAGACCTgg                                      >  2:668475/6‑140 (MQ=11)
                                                                                         tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                      <  1:797411/136‑1 (MQ=11)
                                                                                         tattgtaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCAAATCACTAATATATCATGTCAGTGTTCGTTTAACTTACTAGCATGATGCTAATTTCGTGGCATAGCGAAACTCCGTAAAACTGCATAGCCCTGc                                      >  2:437670/6‑141 (MQ=2)
                                                                                          taacTAATAAAAATAGAAGGTTCTATTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                       <  2:149111/138‑1 (MQ=255)
                                                                                          attttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                     <  2:483250/137‑1 (MQ=11)
                                                                                          attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                     <  2:206346/137‑1 (MQ=11)
                                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCGCTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                    >  2:278790/4‑141 (MQ=11)
                                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGCgg                                    >  1:341215/4‑141 (MQ=11)
                                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                    >  1:492142/4‑141 (MQ=11)
                                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTAAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                    <  2:134401/138‑1 (MQ=11)
                                                                                           ttttaataaAAATGATGGTTCTGCTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                    <  1:152758/138‑1 (MQ=11)
                                                                                           aTCTATTAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                      <  2:609559/141‑1 (MQ=255)
                                                                                           aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                      <  1:102870/141‑1 (MQ=255)
                                                                                           aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                      <  1:206777/141‑1 (MQ=255)
                                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAACCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCTACAATCCGTAGATCTGAAGAGGCGTGCggg                                   >  2:534774/3‑140 (MQ=11)
                                                                                             ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCGGTAGAACTGAAGAGACCGGCGGtt                                  >  2:496364/2‑141 (MQ=11)
                                                                                             ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                  <  2:49125/140‑1 (MQ=11)
                                                                                              taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                   <  2:655806/141‑1 (MQ=40)
                                                                                              taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                 <  1:427491/141‑1 (MQ=11)
                                                                                              taataaAAATGATGGTTCAGTTGCAAAGTAAACAATATCGCTACTCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTATCATGATGCTAATTTCGTGGCATGTCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                 <  2:827263/141‑1 (MQ=11)
                                                                                               aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                >  2:127452/1‑141 (MQ=11)
                                                                                               aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTATTCGTTTAACTTGCTAGCATGATGCTAATTTCGTCGCATGCCGAAAATCCGTAGATCTGAAGAGCCCTCGGGTTCt                                >  2:174776/1‑141 (MQ=11)
                                                                                                ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                               >  2:21499/1‑141 (MQ=11)
                                                                                                 taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                              <  1:84329/141‑1 (MQ=11)
                                                                                                 aaaaaaaTGAAGGTTCTGATGCAAAGTAAAAAATATAGCTAATCACTAATATAACATATAAGAGAACGTTTAACTTGCTAACATGTTGCTAATACCGAGGACTTGCGAAAATCCGTAGATCTGAAGAGAACTGCGGTTCtt                               >  2:250293/2‑141 (MQ=2)
                                                                                                    aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGTTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                             <  1:166336/141‑1 (MQ=25)
                                                                                                       tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCCGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat                          <  1:679145/141‑1 (MQ=17)
                                                                                                       tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat                          <  2:263973/141‑1 (MQ=32)
                                                                                                        tgAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACGAATTTGTCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                         <  2:643114/140‑1 (MQ=14)
                                                                                                         gAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCAAGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtctat                        >  1:478299/1‑141 (MQ=14)
                                                                                                           agcttCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCCTGCCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATag                      <  2:764303/138‑1 (MQ=0)
                                                                                                              ttatGTTGCAAAGTAAAAAATATAGGTAATCACTATTATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCg                  <  2:259477/138‑1 (MQ=1)
                                                                                                                               aaaagctaGCTAACCACTTATTTATCATGTCAGTAGTCGCTTAACTTGCTAGCATCATGCTAATTTCGTGTCATGGCGAAAATCCTTAGATCTGAAGAGACCTGCGGTTCTTTTTACATAGAGCGTAAATACACTCAATAc  <  2:671899/135‑1 (MQ=1)
                                                                                                                                            |                                                                                                                               
CTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATAC  >  CP000730/35710‑35977

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: