Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 46.6% 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,8580TA46.6% 98.0 / 110.3 90intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (19/29);  new base A (14/28);  total (33/57)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35719‑35988
                                                                                                                                           |                                                                                                                                  
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:832205/141‑1 (MQ=255)
          aTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                                                                         <  1:519546/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                                                                                           >  1:170303/1‑141 (MQ=255)
                                        gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                                                           >  2:270981/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                                                               >  1:821036/1‑141 (MQ=255)
                                                                    aaaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAt                                                               >  2:95281/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:413635/134‑1 (MQ=11)
                                                                              agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                   <  1:375906/134‑1 (MQ=11)
                                                                              agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                   <  1:413813/134‑1 (MQ=11)
                                                                              agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATCCCGAAAATCCGTAGATCTGAAGAGACCt                                                   <  2:464790/134‑1 (MQ=11)
                                                                               gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  1:733475/135‑1 (MQ=11)
                                                                               gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  1:698599/135‑1 (MQ=11)
                                                                               gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  >  1:496512/7‑141 (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)
                                                                                                                                caatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGGGCGTAAATACACTCAATACCTTTTAAAGt   <  2:323360/137‑1 (MQ=1)
                                                                                                                                 aatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCGAGCATGAAGCTCATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa  <  1:110147/138‑1 (MQ=1)
                                                                                                                                 aatcACTAATATATAATGTCAGTGTTCGTTTAACTTGTTAGTATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa  <  1:57181/138‑1 (MQ=1)
                                                                                                                                           |                                                                                                                                  
AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35719‑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: