Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,814 Δ1 bp 42.7% intergenic (+45/+63) 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,8140A.42.7% 186.0 / 198.0 117intergenic (+45/+63)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base . (15/35);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 4.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 119 total aligned reads displayed.

TGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATA  >  CP000730/35676‑35954
                                                                                                                                          |                                                                                                                                            
tGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGa                                                                                                                                            >  1:60335/1‑141 (MQ=255)
 gCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGTTCTAATAAAAATAGaa                                                                                                                                           >  2:531005/1‑141 (MQ=255)
        attCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTg                                                                                                                                    >  2:512420/1‑141 (MQ=255)
           cTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTg                                                                                                                                 >  1:296030/1‑141 (MQ=255)
            tGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGc                                                                                                                                >  1:162033/1‑141 (MQ=255)
            tGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGc                                                                                                                                >  1:220818/1‑141 (MQ=255)
                   aaaGATCTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaa                                                                                                                         >  2:368924/1‑141 (MQ=255)
                     aGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaa                                                                                                                       >  2:129493/1‑141 (MQ=255)
                      gATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGGaaaaa                                                                                                                      >  2:321017/1‑141 (MQ=255)
                       aTGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaaaa                                                                                                                     >  2:177163/1‑141 (MQ=255)
                        tGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaaaaa                                                                                                                    <  1:665204/141‑1 (MQ=255)
                        tGTAGATGATCTCATTAACAAAGTATTTAAGGTAGAAGCTATCCAGTATAAATTTTTTGGCGAAAAAATATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaaaaa                                                                                                                    <  2:285771/141‑1 (MQ=255)
                                 tctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTaa                                                                                                           >  2:175675/1‑141 (MQ=255)
                                 tctcATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTaa                                                                                                           <  1:338017/141‑1 (MQ=255)
                                                  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)
                                                                  tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                          <  1:175675/141‑1 (MQ=255)
                                                                     gTATAAAGGTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAATTAGAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                       <  1:129493/141‑1 (MQ=255)
                                                                       aaaaacttataggCGAAAAATGATTACTTCAACAAGAGCTAAAATAAAAATTTTATCAAATAAATATATAATGTATTTTTAAAAAATAAAAAAATATAGCTAACCACTAATTAATCATGTCAGTGTTCGCTTAACTTGCTa                                                                     <  2:359211/130‑1 (MQ=12)
                                                                             tttttGGCGAAAAATGATTACTTGAACAAGGGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGcctca                                                               >  2:564176/1‑137 (MQ=255)
                                                                              ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                              <  1:214941/141‑1 (MQ=255)
                                                                              ttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgat                                                              <  1:358663/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)
                                                                                                 cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGc                                           >  2:583744/1‑141 (MQ=255)
                                                                                                     aaCAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaaa                                       <  2:162033/141‑1 (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)
                                                                                                              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:289535/128‑1 (MQ=11)
                                                                                                                   caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                       <  1:198259/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)
                                                                                                                     aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                     >  1:104310/12‑141 (MQ=11)
                                                                                                                     aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                     >  1:641607/12‑141 (MQ=11)
                                                                                                                     aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGAGGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                     >  2:536122/12‑141 (MQ=11)
                                                                                                                       agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATTGCTAATCACTAATTTATCATGTCAGTTTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  2:296497/132‑1 (MQ=11)
                                                                                                                       agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:634279/132‑1 (MQ=11)
                                                                                                                       agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:458466/132‑1 (MQ=11)
                                                                                                                       agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:658730/132‑1 (MQ=11)
                                                                                                                       agactattttaataaACATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTCAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  2:347629/132‑1 (MQ=11)
                                                                                                                         tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCATTTATCATGTCAGTGTGCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                   <  1:254931/141‑1 (MQ=255)
                                                                                                                         agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                 <  2:200808/134‑1 (MQ=11)
                                                                                                                         agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                 <  1:611719/134‑1 (MQ=11)
                                                                                                                          gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                <  1:437345/135‑1 (MQ=11)
                                                                                                                          gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                <  2:107916/135‑1 (MQ=11)
                                                                                                                          gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                <  1:552215/135‑1 (MQ=11)
                                                                                                                           tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  1:663003/136‑1 (MQ=11)
                                                                                                                           tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  1:541181/136‑1 (MQ=11)
                                                                                                                           tattttaataaAAATGATGATTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc               <  1:635151/136‑1 (MQ=11)
                                                                                                                            gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                >  1:659221/1‑141 (MQ=255)
                                                                                                                            gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                >  1:554976/1‑141 (MQ=255)
                                                                                                                            ccgtcatattaaaTAGCTGGTTCTGTTCCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCtac               <  1:527764/131‑3 (MQ=11)
                                                                                                                            atttttataaAAATGATGGTTCTGTTGCGAAGTAAAAAATATAGCTAATCTCTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  2:27776/135‑1 (MQ=11)
                                                                                                                            attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              >  2:442275/5‑141 (MQ=11)
                                                                                                                            attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              >  2:331131/5‑141 (MQ=11)
                                                                                                                            attttaataaAAAGGATGGTTCTGTTGCAAAGTCAAAAATAAAGATAATTAACTGTTTAGCATGTCAGTGTTCGTTTAACCTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  1:288862/137‑1 (MQ=2)
                                                                                                                            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)
                                                                                                                              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)
                                                                                                                              tttaagaaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt            <  2:374664/139‑1 (MQ=11)
                                                                                                                              tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg              <  1:512420/141‑1 (MQ=255)
                                                                                                                                taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc          >  2:290786/1‑141 (MQ=11)
                                                                                                                                 aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAATAATTTATCATGTCAGTATTCGCTTAACTTGCTAGCATGAAGCTAATTTCGTGGCACGGCACACATCCGTAGACCTGAAGTGAGCcgcgcgc           >  2:348392/1‑138 (MQ=14)
                                                                                                                                 aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt         >  2:269111/1‑141 (MQ=11)
                                                                                                                                 aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt         >  2:29345/1‑141 (MQ=11)
                                                                                                                                  ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt        >  1:567860/1‑141 (MQ=11)
                                                                                                                                      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:567056/1‑141 (MQ=32)
                                                                                                                                          aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata  >  1:1103/1‑141 (MQ=32)
                                                                                                                                          |                                                                                                                                            
TGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATA  >  CP000730/35676‑35954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: