Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,876 C→T 61.2% intergenic (+107/+1) 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,8760CT61.2% 10.5 / 141.4 157intergenic (+107/+1)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base T (25/67);  total (59/98)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 160 total aligned reads displayed.

AGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACC  >  CP000730/35736‑36016
                                                                                                                                            |                                                                                                                                            
 gCTATTAAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                                                                                             <  2:383991/141‑1 (MQ=255)
      tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                                                                                        <  1:175675/141‑1 (MQ=255)
         gTATAAAGGTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAATTAGAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                                                                     <  1:129493/141‑1 (MQ=255)
                 tttttGGCGAAAAATGATTACTTGAACAAGGGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGcctca                                                                                                                             >  2:564176/1‑137 (MQ=255)
                  ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                                                                                            <  1:358663/141‑1 (MQ=255)
                  ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                                                                                            <  1:214941/141‑1 (MQ=255)
                   tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                                                           >  2:192729/1‑141 (MQ=255)
                      ggCGAAAAATGATTACTTGAACAAGAGCTAAAATACAATTGTGAACTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGATCACTAACAATTGCTCATGTAAGTATACGCTTAGCTTGCTAGAGTGATGCTa                                                                                                                        >  2:398352/1‑141 (MQ=21)
                           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)
                                         aaCAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaaa                                                                                                     <  2:162033/141‑1 (MQ=255)
                                         aaCAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTACCTTGGCATGGTTCGATCAACGGTGTAGCATGATGCTAATTTCGTGGCATGGCGaaa                                                                                                     <  1:157358/141‑1 (MQ=18)
                                         aaCAAGAGCTAAAATAAAATCGTGATCTAATAAAAATAGAAGGTTAAGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAGTTTCGTGGCATGGAGaaa                                                                                                     >  2:508716/1‑141 (MQ=18)
                                               agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                                               <  1:192729/141‑1 (MQ=255)
                                                  taaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGa                                                                                            >  1:537439/1‑141 (MQ=255)
                                                      acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                      >  2:147566/15‑141 (MQ=11)
                                                       caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAGATCCGTAGATCTGAag                                                                                     <  1:396503/128‑1 (MQ=11)
                                                       caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                     <  1:350364/128‑1 (MQ=11)
                                                       caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                     <  1:289535/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)
                                                          aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                                                  <  2:55909/131‑1 (MQ=11)
                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATTGCTAATCACTAATTTATCATGTCAGTTTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                 <  2:296497/132‑1 (MQ=11)
                                                           agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                 <  1:458466/132‑1 (MQ=11)
                                                           agactattttaataaACATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTCAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                 <  2:347629/132‑1 (MQ=11)
                                                             tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCATTTATCATGTCAGTGTGCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                 <  1:254931/141‑1 (MQ=255)
                                                             agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                                               <  1:611719/134‑1 (MQ=11)
                                                             agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                                               <  2:200808/134‑1 (MQ=11)
                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                              <  2:107916/135‑1 (MQ=11)
                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                              <  1:552215/135‑1 (MQ=11)
                                                              gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                              <  1:437345/135‑1 (MQ=11)
                                                               tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                             <  1:541181/136‑1 (MQ=11)
                                                               tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                             <  1:663003/136‑1 (MQ=11)
                                                               tattttaataaAAATGATGATTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                             <  1:635151/136‑1 (MQ=11)
                                                                gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                              >  1:659221/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)
                                                                 aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                             <  1:531005/141‑1 (MQ=255)
                                                                  tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGAGCTGAAGAGACCTGCGGt                                                                          <  2:144755/139‑1 (MQ=11)
                                                                  tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                          <  2:171271/139‑1 (MQ=11)
                                                                  tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                          <  1:106057/139‑1 (MQ=11)
                                                                  tttaataaAAATGATGGTTCCTTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                          <  2:362816/139‑1 (MQ=11)
                                                                   cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                           >  2:515911/1‑141 (MQ=40)
                                                                    taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                                                        >  2:290786/1‑141 (MQ=11)
                                                                     aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                       >  2:269111/1‑141 (MQ=11)
                                                                      ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                      >  1:567860/1‑141 (MQ=11)
                                                                          aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                                                    >  1:537689/1‑141 (MQ=35)
                                                                          aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                                                    >  2:55809/1‑141 (MQ=35)
                                                                          aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                                                    >  2:598979/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)
                                                                                aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata                                                              <  1:125442/141‑1 (MQ=11)
                                                                                     tcagTTCCAAACTAAAAAAACAGCGCTTGCCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCg                                                         <  2:520221/138‑1 (MQ=1)
                                                                                                         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)
                                                                                                                    cgttatatATCATGTCAGTGTTCGTTTAACTTGCTAGGATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCt                          <  2:368478/134‑1 (MQ=1)
                                                                                                                    cgtaatatatctgttCATTGTTCCTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGAGCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCt                          <  2:540384/127‑1 (MQ=1)
                                                                                                                        atataccATGTCAGTGGTCGGTTAACTGTCTAGCATGATGCTAATTTCGTAGCATGTCGAAAAACCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGc                      <  2:256031/135‑1 (MQ=1)
                                                                                                                        atatATCATGGCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTACATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGc                      <  1:364706/138‑1 (MQ=1)
                                                                                                                             tcttgccAGTTTTCGCTTAACTTGGGAGCATGATGATAATTTGGCGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTGTAGTATTCTTTGCCGTAt                 <  2:427353/135‑1 (MQ=1)
                                                                                                                                   cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAc           <  1:628985/141‑1 (MQ=2)
                                                                                                                                   cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc           <  2:356243/141‑1 (MQ=11)
                                                                                                                                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt          <  1:521340/141‑1 (MQ=255)
                                                                                                                                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt          <  1:336484/141‑1 (MQ=11)
                                                                                                                                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCAGGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt          >  2:332866/1‑141 (MQ=9)
                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTTCCGTATTGATACtt         <  1:435311/141‑1 (MQ=11)
                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt         >  1:16724/1‑141 (MQ=255)
                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTAGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACATTTTAAAGTATTATTTGCTGTATTAATACtt         >  2:579099/1‑141 (MQ=11)
                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGGTGCTAATTTCGTGGCATGGCGAAAATCCGTAGGTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt         >  2:589188/1‑141 (MQ=12)
                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTGTGTCGAGCGTAAATAGACACAATCCCTTTTAAAGTATTCTTTGCGGTATTTATACtt         >  1:285016/1‑141 (MQ=11)
                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt         >  1:430628/1‑141 (MQ=17)
                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGACACtt         <  1:454738/141‑1 (MQ=14)
                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTGTCTGTATTGATACtt         >  2:104631/1‑141 (MQ=12)
                                                                                                                                      tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt        <  2:500735/141‑1 (MQ=255)
                                                                                                                                      gtgtCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt        <  2:7907/138‑1 (MQ=11)
                                                                                                                                       gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg       >  2:381173/1‑141 (MQ=255)
                                                                                                                                       gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTCTTAAAGTATTCTTTGCCGTATTGATACTTTg       >  2:566958/1‑141 (MQ=11)
                                                                                                                                       gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg       >  2:18684/1‑141 (MQ=12)
                                                                                                                                        ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa      <  2:46492/141‑1 (MQ=17)
                                                                                                                                        ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATGCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa      <  2:664180/141‑1 (MQ=17)
                                                                                                                                        ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa      >  1:269120/1‑141 (MQ=17)
                                                                                                                                          cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa    >  2:44753/1‑141 (MQ=32)
                                                                                                                                            cTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc  <  1:570813/141‑1 (MQ=17)
                                                                                                                                            |                                                                                                                                            
AGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACC  >  CP000730/35736‑36016

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: