Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 26.2% T25A (ACA→GCA)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC26.2% 107.9 / 46.3 86T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base C (9/13);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618
                                                                                                                                           |                                                                                                                                            
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtt                                                                                                                                             >  2:535678/1‑140 (MQ=255)
 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCATTGtaa                                                                                                                                            >  2:268896/1‑140 (MQ=255)
 aataaGCTTCTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtat                                                                                                                                            >  1:577815/1‑141 (MQ=255)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAGCGTTTACACCACGTTCCCTTAATATTTCAGATATATCATGATAACTCAATGTATATCt                                                                                                                                        <  1:473852/141‑1 (MQ=255)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                        <  2:533888/141‑1 (MQ=38)
            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACAGCACGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                                 >  1:582588/1‑141 (MQ=25)
            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                                 >  2:337639/1‑141 (MQ=35)
             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:467483/1‑141 (MQ=35)
             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:96053/1‑141 (MQ=35)
             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACAGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:512139/1‑141 (MQ=25)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                               >  1:94844/1‑141 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                               >  2:332236/1‑141 (MQ=35)
                gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                             >  1:564426/1‑141 (MQ=35)
                  tcttgcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCa                                                                                                                           <  1:241332/136‑1 (MQ=25)
                          ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   <  2:49111/141‑1 (MQ=21)
                          ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTTTATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   >  1:599114/1‑141 (MQ=21)
                          ccAAAGTTGATACAGAGTGTGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   <  2:123005/141‑1 (MQ=11)
                           ttttttttATACAAAATTGGGGCATAATCTTGTACCCAACGGTAGACCGTTGAATGGTGAACGTTTACACCACTTTTTTTTATTATTTCAGATATATGATGATAACTCAATTTATATTTTAGATAGTAGCCAACgggtgaa                                                                                                                  >  2:495916/5‑136 (MQ=11)
                                  gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAGGGCTACAGTgctca                                                                                                           >  1:521911/1‑137 (MQ=255)
                                         aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATGCtt                                                                                                   >  2:313297/1‑141 (MQ=17)
                                           ttGGGGCATATTCTTGAACCCAACGGTATACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:447623/141‑1 (MQ=17)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:85203/141‑1 (MQ=32)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCCACAGTG‑ATAACATCCTTGt                                                                                                 <  1:55735/141‑1 (MQ=17)
                                           tgggggCATATTCTTGAACCCACCCGTAGACCGTTGAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:162038/139‑1 (MQ=12)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTaa                                                                                              >  1:320265/1‑141 (MQ=32)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATTTAGATCTTAGATAGTAGCCAACGGCTACAGTG‑CTACGATCCTTGTTaa                                                                                              >  1:22184/1‑141 (MQ=11)
                                               ggCCTATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATACCACGATAACTCAATGTATGTCTTAGATAGTAGCCACCGGCCACAGTG‑ATAACCTCCCTGTTaca                                                                                             >  2:92202/1‑139 (MQ=11)
                                               ggCAAATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGGTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGCTaaa                                                                                             >  2:545530/1‑141 (MQ=12)
                                               agCATATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGATCCCTTCATATTTCAGAAATATCACGATAACTCACTGTAAATCTTAGATCATGGACAACGCCT‑CAATGAATAACATCCTTGTTcaa                                                                                             >  2:92214/2‑141 (MQ=11)
                                                 cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATATATGTCTTAGATATTAGCCAACGGCTACAGTG‑ATAACATCCTAGTTAAAtt                                                                                           >  2:754851/1‑141 (MQ=11)
                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTtat                                                                                     >  2:2533/1‑141 (MQ=17)
                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATACCATCGTTGTTAAATTGTTtat                                                                                     <  1:405336/141‑1 (MQ=17)
                                                           aaCCCATCCGTAGATTGTTTAATGATGAACGTTTACACGACGTCCTCCTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCt                                                                                 <  2:537633/141‑1 (MQ=11)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGa                                                                               >  2:394832/1‑141 (MQ=32)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAGATATTGCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGAG‑ATAACATCCTTTTTAACTTGTTTATATCTGa                                                                               >  1:327676/1‑141 (MQ=11)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACAACACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATGTATACTTTAGATAGTAGGCAACGGCTACAGTG‑ATAACATCCTTGTTAACTTGTTTAGATCTGa                                                                               >  2:329954/1‑141 (MQ=11)
                                                               cAACCGTAGACCCTTGAATAATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                             <  2:468732/141‑1 (MQ=12)
                                                                      aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                      <  2:96053/141‑1 (MQ=32)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                     >  1:303520/1‑141 (MQ=32)
                                                                            ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                                >  2:773191/1‑141 (MQ=32)
                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttt                                                     >  2:581495/1‑141 (MQ=17)
                                                                                           ttACACCACGTTCCTTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  2:610387/141‑1 (MQ=17)
                                                                                           ttACACCACCTTCCCTTAACATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:28769/141‑1 (MQ=17)
                                                                                           tgacacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:14448/139‑1 (MQ=17)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:124833/141‑1 (MQ=32)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:129676/141‑1 (MQ=32)
                                                                                               cccccGTTCTCTTAGTATTTCAAATATACCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                             <  1:266322/137‑1 (MQ=1)
                                                                                                ccacGTTCCCTTAATATTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                            <  1:23867/141‑1 (MQ=17)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                         >  2:18857/1‑141 (MQ=32)
                                                                                                      tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      >  2:497075/1‑141 (MQ=32)
                                                                                                      tCCCTTAACATATCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      <  2:507431/141‑1 (MQ=17)
                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  1:632673/141‑1 (MQ=35)
                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  2:157457/141‑1 (MQ=35)
                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                          >  1:440501/1‑141 (MQ=35)
                                                                                                                   cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                         >  2:773679/1‑141 (MQ=35)
                                                                                                                     aatatatCACGATAACTCAATGGATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                       <  2:144553/140‑1 (MQ=17)
                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATAATATAAATTTAACTTTGCa                     >  2:471689/1‑141 (MQ=255)
                                                                                                                       tatatCACGATAACTCAACGGATATCGTAGATAGAAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATGCAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACTTTGCa                     <  1:202110/141‑1 (MQ=11)
                                                                                                                       taaatCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  1:353346/138‑1 (MQ=21)
                                                                                                                         tatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                   <  1:416965/141‑1 (MQ=255)
                                                                                                                          atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                  <  2:577815/141‑1 (MQ=35)
                                                                                                                              cGATAACTCAATGTATATCTTAGATAGTAGCCATCGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc              >  2:520093/1‑141 (MQ=25)
                                                                                                                               gatagcgccaTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcc             <  2:456042/132‑1 (MQ=1)
                                                                                                                               gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAccc             >  2:193791/1‑141 (MQ=25)
                                                                                                                                 tagcTCAGTGCATATCTGATATAGGAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTTGGATCTTAAATAGTTCAGACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  1:550511/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:242193/138‑1 (MQ=12)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:76936/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCACATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACCCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAccacacctt           <  2:683328/138‑9 (MQ=11)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:756098/3‑141 (MQ=12)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:640460/3‑141 (MQ=9)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:231193/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:788733/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCAAACGGCTACAGTG‑ATAACATCCTGGTTAAATTGTTTATATCTGAAATAGATCATACACAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACTGAACCGtat         >  2:190875/2‑141 (MQ=11)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:772242/1‑141 (MQ=17)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        <  2:702947/141‑1 (MQ=18)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:171321/1‑141 (MQ=16)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:710573/1‑141 (MQ=12)
                                                                                                                                     ccAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta       <  1:66445/140‑1 (MQ=40)
                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATAATTGCTGGACGTTGGATGTCTGAATCCTTGCATACAGAAGACTCCTTTTGGTTAGAAGGATACTATGAATTCAACTTTGCAACAGAACCGtattat      >  1:777847/1‑141 (MQ=9)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:60161/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  2:101883/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:518306/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:129214/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTACTTTTATTTGAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:346955/141‑3 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     >  2:583819/1‑138 (MQ=255)
                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    >  1:451434/1‑141 (MQ=255)
                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    <  1:713354/141‑1 (MQ=255)
                                                                                                                                        aTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACCGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat    <  1:65944/141‑5 (MQ=255)
                                                                                                                                          gTACATCTTAGATAGTAGCCAACCGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATACTTCATACAGAAGACTCCTTCTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta  <  2:410062/141‑12 (MQ=11)
                                                                                                                                           |                                                                                                                                            
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618

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: ‑