Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 52.7% 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,8580TA52.7% 76.5 / 153.1 116intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (16/45);  total (44/72)
Fisher's exact test for biased strand distribution p-value = 7.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

Only 100 of 116 total aligned reads displayed.

TTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35726‑35988
                                                                                                                                    |                                                                                                                                  
ttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGt                                                                                                                            >  1:576591/1‑141 (MQ=255)
  aaTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCa                                                                                                                          >  2:556568/1‑141 (MQ=255)
     gTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAgtg                                                                                                                       >  2:241817/1‑141 (MQ=255)
          aGCTAATAAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                                                                  <  2:562515/141‑1 (MQ=255)
           gCTATTAAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                                                                 <  2:383991/141‑1 (MQ=255)
                  aGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAATAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTg                                                                                                          >  1:162583/1‑141 (MQ=255)
                   gTATAAAGGTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAATTAGAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                                         <  1:129493/141‑1 (MQ=255)
                     ataAAGTTTTTGGCGAAAAATGATTACTTGGACAAGAGCTAAAATAAACTTGTGATCTAATAGAAATATAAGGTTCTATTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTTAGTGTTCGCTTAACTTGCTa                                                                                                       <  2:1971/141‑1 (MQ=25)
                            ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                                                                <  1:358663/141‑1 (MQ=255)
                            ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                                                                <  1:214941/141‑1 (MQ=255)
                             tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAACTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                               <  1:539352/141‑1 (MQ=255)
                             tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                               >  2:192729/1‑141 (MQ=255)
                              ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                                              <  1:139033/141‑1 (MQ=255)
                                ggCGAAAAATGATTACTTGAACAAGAGCTAAAATACAATTGTGAACTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGATCACTAACAATTGCTCATGTAAGTATACGCTTAGCTTGCTAGAGTGATGCTa                                                                                            >  2:398352/1‑141 (MQ=21)
                                  cGAAAAATGATTACTTGAACAAGAGCTAAAGTAAAATTGTAATCTAATAAACATAGAAGGTTATGTTGCAAAGAAAAAAAATATAGCAAACGACTGATTTATCATTAGAGTGTTCGCTTAACTTGCTATCATAGTACTAAt                                                                                          >  2:176185/1‑141 (MQ=21)
                                     aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGGTACTAATTTc                                                                                       >  1:107462/1‑141 (MQ=255)
                                     aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATACTAATTTc                                                                                       >  1:139794/1‑141 (MQ=255)
                                      aaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTTGGGCGTCCTGCTAGCATGATGCTAATTTCg                                                                                      <  1:159747/141‑1 (MQ=21)
                                      aaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTAGATCAACTTGCTAGCATGATGCTAATTTCg                                                                                      <  1:159721/141‑1 (MQ=255)
                                       aaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTATGTTGCAAAGTAAAAAAATATAGATAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGATAGGATCATGTTAATTTCGt                                                                                     >  2:385892/1‑141 (MQ=21)
                                              aCTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATATATCATGTCTATGACCACGGAGCTTGCTAGCATGATGCTAATTTCGTGGCATgg                                                                              <  1:657845/141‑1 (MQ=21)
                                               cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                                                             >  2:583744/1‑141 (MQ=255)
                                                   aaCAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTACCTTGGCATGGTTCGATCAACGGTGTAGCATGATGCTAATTTCGTGGCATGGCGaaa                                                                         <  1:157358/141‑1 (MQ=18)
                                                   aaCAAGAGCTAAAATAAAATCGTGATCTAATAAAAATAGAAGGTTAAGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAGTTTCGTGGCATGGAGaaa                                                                         >  2:508716/1‑141 (MQ=18)
                                                       agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                     >  1:667232/1‑141 (MQ=255)
                                                         agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                   <  1:192729/141‑1 (MQ=255)
                                                           cTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTAAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCCAAAAACCTTac                                                                 <  2:576591/141‑2 (MQ=255)
                                                            taaaataaaaTTGTGATCTAATAAAAATAGGAGGTTCTGTTGCAAAGTCAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGGGCAATTGAGTTGCTAGCATGATGCTAACTTCGTGGCAAGGGGAAAATTCGTAGa                                                                >  1:453624/1‑141 (MQ=17)
                                                            taaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGa                                                                >  1:537439/1‑141 (MQ=255)
                                                                acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                          >  2:147566/15‑141 (MQ=11)
                                                                 caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTACTCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  1:168558/128‑1 (MQ=11)
                                                                 caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAGATCCGTAGATCTGAag                                                         <  1:396503/128‑1 (MQ=11)
                                                                 caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  1:198259/128‑1 (MQ=11)
                                                                 caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  1:350364/128‑1 (MQ=11)
                                                                 caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  1:289535/128‑1 (MQ=11)
                                                                 aaaaatagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                         <  1:456635/128‑1 (MQ=11)
                                                                  aaaagagtattttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                                        <  1:559834/129‑1 (MQ=11)
                                                                  aaaagagtattttaataaAAATGATGGTTCTGTTGTTAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                                        <  1:513680/129‑1 (MQ=11)
                                                                   caatTGTGATCTACTAAAAGTAAAAGGTCCTGGTGCAAATTAGAAAAATATAGCTGTACACTAATTTATGATGTCAGTGTTCGGTAAACTTGCTAGTAAGATTCTAATTTCGTGGCATGGCGAAAATCCGTGGTTCTGaga                                                         >  2:322313/2‑139 (MQ=9)
                                                                   aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                       >  1:104310/12‑141 (MQ=11)
                                                                   aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                       >  1:641607/12‑141 (MQ=11)
                                                                   aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGAGGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                       >  2:536122/12‑141 (MQ=11)
                                                                    aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                      <  2:55909/131‑1 (MQ=11)
                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATTGCTAATCACTAATTTATCATGTCAGTTTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  2:296497/132‑1 (MQ=11)
                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:634279/132‑1 (MQ=11)
                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     >  1:603793/10‑141 (MQ=11)
                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:458466/132‑1 (MQ=11)
                                                                     agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:658730/132‑1 (MQ=11)
                                                                     agactattttaataaACATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTCAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  2:347629/132‑1 (MQ=11)
                                                                      aaagaatttaataaAAATGATGATTCTGTGTCAAAGAAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                    <  2:106640/133‑1 (MQ=11)
                                                                       tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCATTTATCATGTCAGTGTGCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                     <  1:254931/141‑1 (MQ=255)
                                                                        gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  2:107916/135‑1 (MQ=11)
                                                                        gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  <  1:552215/135‑1 (MQ=11)
                                                                         tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:663003/136‑1 (MQ=11)
                                                                         tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  2:280475/136‑1 (MQ=11)
                                                                         tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:541181/136‑1 (MQ=11)
                                                                         tattttaataaAAATGATGATTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:635151/136‑1 (MQ=11)
                                                                          gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                  >  1:554976/1‑141 (MQ=255)
                                                                          atttttataaAAATGATGGTTCTGTTGCGAAGTAAAAAATATAGCTAATCTCTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  2:27776/135‑1 (MQ=11)
                                                                          attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                >  2:331131/5‑141 (MQ=11)
                                                                          attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                >  2:442275/5‑141 (MQ=11)
                                                                          atttgaagaaaaaTGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  2:454307/133‑1 (MQ=11)
                                                                           ttttaataaAAATGATTGTTCTGTTCCAAAGTAAAAAATATAGCTAAACACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               <  2:135722/138‑1 (MQ=11)
                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGGGACATGCgg                                               >  2:275677/4‑141 (MQ=11)
                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               <  1:438781/138‑1 (MQ=11)
                                                                           ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAGCCGTAGATCTGAAGAGACCTGCgg                                               >  1:638948/4‑141 (MQ=11)
                                                                           agtcatattaaaTGGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCACGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCtacg                                                <  2:436980/132‑4 (MQ=11)
                                                                           aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                 <  1:531005/141‑1 (MQ=255)
                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCGTGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTGCGGGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:555063/139‑1 (MQ=11)
                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGAGCTGAAGAGACCTGCGGt                                              <  2:144755/139‑1 (MQ=11)
                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              >  2:394143/3‑141 (MQ=11)
                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:171271/139‑1 (MQ=11)
                                                                            tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  1:106057/139‑1 (MQ=11)
                                                                            tttaataaAAATGATGGTTCCTTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:362816/139‑1 (MQ=11)
                                                                            tttaagaaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                              <  2:374664/139‑1 (MQ=11)
                                                                            tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                <  1:512420/141‑1 (MQ=255)
                                                                             ttaataaAAATGATTGTTCTGGTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTACCTTGCCAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                             <  2:60192/140‑1 (MQ=11)
                                                                             cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                               >  2:515911/1‑141 (MQ=40)
                                                                              taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                            >  2:290786/1‑141 (MQ=11)
                                                                               aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAATAATTTATCATGTCAGTATTCGCTTAACTTGCTAGCATGAAGCTAATTTCGTGGCACGGCACACATCCGTAGACCTGAAGTGAGCcgcgcgc                                             >  2:348392/1‑138 (MQ=14)
                                                                               aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                           >  2:269111/1‑141 (MQ=11)
                                                                                ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                          >  1:567860/1‑141 (MQ=11)
                                                                                   aaaaaggatgggtCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                        <  2:162158/129‑1 (MQ=1)
                                                                                    aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                        >  2:598979/1‑141 (MQ=35)
                                                                                    aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                        >  2:55809/1‑141 (MQ=35)
                                                                                    aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                        >  1:537689/1‑141 (MQ=35)
                                                                                     aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCGCTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                       <  2:357565/141‑1 (MQ=21)
                                                                                       tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat                                     <  1:378710/141‑1 (MQ=32)
                                                                                        aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                    >  1:1103/1‑141 (MQ=32)
                                                                                        aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                    >  1:567056/1‑141 (MQ=32)
                                                                                          aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata                                  <  1:125442/141‑1 (MQ=11)
                                                                                          aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGCCGAAAATCCGTAGATCTGAAGAGACCTGCGTTTCtttttcttta                                  >  1:231785/1‑136 (MQ=2)
                                                                                               tcggTTGCAAAGTAAAAAATATGGCTCTTCAGTAATATATCATGGCAGTCTTCGTTTACCTTGCTAGCATGATGTTATTTTCGTGGCCTGGGGAAACTCCGTAGATCTGAAGAGACCTACGGTTCTCTTTATATAGCGCGt                            <  2:7276/138‑1 (MQ=0)
                                                                                                                   aatcgCTAACCACTAGTTTATCAGGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCttt         <  2:437230/137‑1 (MQ=1)
                                                                                                                      acccaatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTaa      <  1:217924/134‑1 (MQ=1)
                                                                                                                        ctcaCCTCTAACATATCATGTCAGTGTTCGTTTAACCTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAg    <  2:609990/138‑1 (MQ=1)
                                                                                                                         caatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGt   <  1:474603/137‑1 (MQ=1)
                                                                                                                         caatcACTAATATATCAGGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGt   <  2:407266/137‑1 (MQ=1)
                                                                                                                          aatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTATAGtt  >  2:16470/4‑140 (MQ=1)
                                                                                                                          aatcACTAATATATCATGTCAGTGTTAGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTa  <  1:408413/138‑1 (MQ=1)
                                                                                                                                    |                                                                                                                                  
TTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35726‑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: