Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,817 A→T 35.2% intergenic (+48/+60) 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,8170AT35.2% 154.3 / 86.6 99intergenic (+48/+60)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (8/26);  total (41/58)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATA  >  CP000730/35677‑35954
                                                                                                                                            |                                                                                                                                         
gCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGaa                                                                                                                                           >  1:54439/1‑141 (MQ=255)
 cAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAGAGTAATTAAGGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAg                                                                                                                                          >  1:407881/1‑141 (MQ=255)
 cAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGATAAATGATTACTTGAACAAGAGCTAAAATGAAATTGTGATCGAATAAAAATAGAAg                                                                                                                                          >  1:468502/1‑141 (MQ=255)
 aaaaTTATTTCGAATTTAAAGATGGAGATTATCTCATTGAAAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGCTGTTGGCGAAAAATGATTACTCGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAg                                                                                                                                          <  2:401102/140‑1 (MQ=255)
  aaaTTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAgg                                                                                                                                         >  1:284359/1‑141 (MQ=255)
    attattGTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTAAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGtt                                                                                                                                       <  2:74649/141‑1 (MQ=255)
    attattCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGtt                                                                                                                                       <  1:737065/141‑1 (MQ=255)
     ttattCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTc                                                                                                                                      <  1:369417/141‑1 (MQ=255)
     ttattCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTc                                                                                                                                      <  2:137797/141‑1 (MQ=255)
       attCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTg                                                                                                                                    <  1:497440/141‑1 (MQ=255)
         tCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGtt                                                                                                                                  >  2:52254/1‑141 (MQ=255)
                    aGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaa                                                                                                                       >  1:715768/1‑141 (MQ=255)
                        gTAGATGATTTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAat                                                                                                                   <  1:258680/141‑1 (MQ=255)
                            atgatCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTTTTTTTTTTTGATTACTTGAACAAGAGCTAAAATAAAATGGTGATGTAAGAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAg                                                                                                               >  2:692066/1‑141 (MQ=255)
                             tgatCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAACATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGc                                                                                                              >  1:106302/1‑141 (MQ=255)
                             agatTTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGc                                                                                                              <  2:34274/140‑1 (MQ=255)
                              aatCTCATTGACAAATTATTTAAAGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCt                                                                                                             <  2:284359/140‑1 (MQ=255)
                                tctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTCTTGCAAAGTAAAAAAATATAGCTaa                                                                                                           <  1:111965/141‑1 (MQ=255)
                                     ttGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACt                                                                                                      <  1:317586/141‑1 (MQ=255)
                                           aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCAGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                >  1:567609/1‑141 (MQ=255)
                                            aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAACATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAGTAATTTAt                                                                                               <  2:570793/141‑1 (MQ=255)
                                            aGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAGGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTCAt                                                                                               >  2:735321/1‑141 (MQ=255)
                                                  tAATGTAGAATTTATCAAGTATAGCGTTGTGTGCTAAAAAGGATTCCTTGAACAAGTCCTAAAATAAAATTGTAATCTAATACCCATATAAGGTTCTAAAAAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTc                                                                                         <  2:313693/141‑1 (MQ=12)
                                                  tAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAATAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTGACCACTAATTTATCaagac                                                                                         >  1:221268/1‑137 (MQ=255)
                                                      gTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCtgtg                                                                                     >  1:419738/1‑141 (MQ=255)
                                                          aaGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCg                                                                                 <  1:775174/141‑1 (MQ=255)
                                                           aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                                >  2:52105/1‑141 (MQ=255)
                                                            gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTACTAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                               >  1:463112/1‑141 (MQ=255)
                                                                 tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                          >  1:691791/1‑141 (MQ=255)
                                                                            tttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAACGAAACAAAATATAGCTAACCACTAATTGATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                               <  1:145606/141‑1 (MQ=255)
                                                                            tttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTGGCatga                                                               >  1:314522/1‑141 (MQ=255)
                                                                            gttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                               <  1:171141/140‑1 (MQ=255)
                                                                              tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAACTATAGCTAACCACTAATTTATCATGTGTGTGGTCGATTAACTTGCTAGCatgctg                                                             >  1:21021/1‑141 (MQ=25)
                                                                                     aaaaaTGATTACTTTAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAAttt                                                      >  1:426570/1‑141 (MQ=255)
                                                                                           gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                <  1:390212/141‑1 (MQ=255)
                                                                                           gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGc                                                >  1:287414/1‑141 (MQ=255)
                                                                                                cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                           <  1:646081/141‑1 (MQ=255)
                                                                                                cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGGAAAAAAATATAGCTTACCACTGATTTATCATGTCAGTGTTCGCGTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                           <  1:57565/141‑1 (MQ=255)
                                                                                                         gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTGATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                  >  1:522163/1‑141 (MQ=255)
                                                                                                         gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGCTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCACGGCGAGAATCCg                                  >  2:126760/1‑141 (MQ=255)
                                                                                                         gagCTAAAAGAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                  <  2:54439/141‑1 (MQ=255)
                                                                                                         gaacTAAAATAAAATTGTGATATAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                  <  2:266392/138‑1 (MQ=255)
                                                                                                          agCTAAAATAAAATTGTGATGTAATAAAAATAGAAGGTTCTGTTCCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGTATGATGCTAATTTCGTGGTATGGCGAAAATCCGt                                 >  1:469132/1‑141 (MQ=255)
                                                                                                          agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                 >  2:289826/1‑141 (MQ=255)
                                                                                                          agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                 >  2:624081/1‑141 (MQ=255)
                                                                                                          agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                 >  2:289982/1‑141 (MQ=255)
                                                                                                          agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                 >  1:418726/1‑141 (MQ=255)
                                                                                                            gtaaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTTCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAg                               <  2:36323/140‑1 (MQ=255)
                                                                                                            caaaaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTACGTGGCATGGCGAAAATCCGTAg                               <  2:553619/139‑1 (MQ=255)
                                                                                                            ataaaaaaaaaTTGTGATCTAATAAAAATAGAAGGCTCTGTTGCAAAGTAAAAAAAGATAGCTAACCACTAATTTATCATGGCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAg                               <  2:532436/140‑1 (MQ=255)
                                                                                                                 acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                        >  1:440208/15‑141 (MQ=11)
                                                                                                                 acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                        >  1:675577/15‑141 (MQ=11)
                                                                                                                    aaatTGTGATCTAATAAAAATAAAAAGTTCTTTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                       <  2:345091/141‑1 (MQ=255)
                                                                                                                    aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAAACACTAATATATCATGACAGAGTTCGTTTCACTTGCCAGCATAGTGCTACTTTCGTGGCATGGCGAAAAACGGTAGATCTCAaggg                     >  2:532023/12‑139 (MQ=2)
                                                                                                                     aatTGGTATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                      <  2:143856/141‑1 (MQ=255)
                                                                                                                      atagtagtttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAACATATCATGCCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCTTAGATCTGAAGAGAc                   <  2:488477/132‑1 (MQ=11)
                                                                                                                      agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:754523/132‑1 (MQ=11)
                                                                                                                      agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:263126/132‑1 (MQ=11)
                                                                                                                      agagtaattcaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATTTAGATTACCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:682058/131‑1 (MQ=11)
                                                                                                                       gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                  >  2:402488/9‑141 (MQ=11)
                                                                                                                       gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                  >  2:402460/9‑141 (MQ=11)
                                                                                                                       gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                  >  1:766973/9‑141 (MQ=11)
                                                                                                                        agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                 <  2:672182/134‑1 (MQ=11)
                                                                                                                         gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                <  2:185896/135‑1 (MQ=11)
                                                                                                                         gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                <  2:72156/135‑1 (MQ=11)
                                                                                                                         gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAACCCGTAGATCTGAAGAGACCTg                <  2:94221/135‑1 (MQ=11)
                                                                                                                          tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  2:10339/136‑1 (MQ=11)
                                                                                                                          tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCGAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  2:155131/136‑1 (MQ=11)
                                                                                                                          tattttaataaAAATGAGGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  1:679578/136‑1 (MQ=11)
                                                                                                                          ggATCGAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                 <  1:232169/140‑1 (MQ=255)
                                                                                                                           attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  1:642688/137‑1 (MQ=11)
                                                                                                                           attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  2:211134/137‑1 (MQ=11)
                                                                                                                           attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  1:659376/137‑1 (MQ=11)
                                                                                                                            ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg             <  2:36853/138‑1 (MQ=11)
                                                                                                                            ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg             >  2:524611/4‑141 (MQ=11)
                                                                                                                            aTCTAATAAAAGTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  2:114978/141‑1 (MQ=255)
                                                                                                                            aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  1:378033/141‑1 (MQ=255)
                                                                                                                             tttaataaAAATGATTGTTCTGTTGCAAAGTAAAAAATATACCCAATCCCAAATACTTCATGTAAGTCTTTGTTTAACTTGCTAGCATTATGCTAATTTAGTGGCATGGCGAAAATCCGTAGACCTGAAGAGACCTGCGGt            <  2:193526/139‑1 (MQ=2)
                                                                                                                             tttaataaAAATGATGTTTCTGTTGCAAAGTAACAAATATAGCTAATCCCCATTATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAACTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            <  1:349562/139‑1 (MQ=11)
                                                                                                                             tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            <  2:673339/139‑1 (MQ=11)
                                                                                                                             tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            <  2:646707/139‑1 (MQ=11)
                                                                                                                             ttgaataaAAATGATGCTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATAGCATGTCAGTGTTCGTTTAACTTGCTAGCATGTTGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            <  1:245161/138‑1 (MQ=11)
                                                                                                                              ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt           >  1:36603/2‑141 (MQ=11)
                                                                                                                              ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt           <  1:376605/140‑1 (MQ=11)
                                                                                                                              ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAAGCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt           <  2:682348/140‑1 (MQ=11)
                                                                                                                              cTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAAAGATATAGCTAACCACTAATTTATCAGTGCAGTGGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGAgg             >  2:261106/1‑141 (MQ=14)
                                                                                                                              cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTACTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg             >  1:287213/1‑141 (MQ=255)
                                                                                                                               taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            >  2:256714/1‑141 (MQ=40)
                                                                                                                               taataaAAATGATGGTTGTGTTGCAAAGTCAAAAATCTAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc          <  1:3529/141‑1 (MQ=11)
                                                                                                                               taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc          <  1:471129/141‑1 (MQ=11)
                                                                                                                               taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc          <  1:587593/141‑1 (MQ=11)
                                                                                                                               taataaAAATGATGGTTATGCTGCAAAGTAAAAAATATAGATAATCACTAAGATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCGGAAGAGACCTGCGGTTc          <  1:107998/141‑1 (MQ=11)
                                                                                                                                aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGACCTGAAGGGCCCTGCGGtt           >  2:142565/1‑141 (MQ=255)
                                                                                                                                aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCACTTAACTTGCTAGCATGATGCTAATTTCGTGGCAAGGAGAGAATCAGTAGTTGTTAAGAGACGTGCGAtt           >  2:36613/1‑141 (MQ=14)
                                                                                                                                aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt         <  1:592031/141‑1 (MQ=11)
                                                                                                                                  taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGCAAATCCGTAGATCTGAAGAGACCTGCGGTTCt         >  1:108364/1‑141 (MQ=38)
                                                                                                                                     aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt      <  1:541634/141‑1 (MQ=35)
                                                                                                                                     aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCGTAACTTGCTAGCATGATGCTAATTTCGTGGCATGTCGAAATTCCGTGTATCTGAAGAGACCTGCGGTTCtttt      <  1:735321/141‑1 (MQ=11)
                                                                                                                                      aaTAGAAGGTTATGATGCAAAGAACCAAAATATAGCTAACTACTAATTTATCATGTCAGGGTTCGCTTAATTTGCGAGCATGATTCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt     >  2:534753/1‑141 (MQ=11)
                                                                                                                                         aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata  >  1:687715/1‑141 (MQ=32)
                                                                                                                                            |                                                                                                                                         
GCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATA  >  CP000730/35677‑35954

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: