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

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  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)
                                         aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCtt                                                                                                   >  2:313297/1‑141 (MQ=17)
                                           ttGGGGCATATTCTTGAACCCAACGGTATACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGt                                                                                                 <  1:447623/141‑1 (MQ=17)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                 <  1:85203/141‑1 (MQ=32)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCCACAGTGATAACATCCTTGt                                                                                                 <  1:55735/141‑1 (MQ=17)
                                           tgggggCATATTCTTGAACCCACCCGTAGACCGTTGAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                 <  1:162038/139‑1 (MQ=12)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                              >  1:320265/1‑141 (MQ=32)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATTTAGATCTTAGATAGTAGCCAACGGCTACAGTGCTACGATCCTTGTTaa                                                                                              >  1:22184/1‑141 (MQ=11)
                                               ggCCTATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATACCACGATAACTCAATGTATGTCTTAGATAGTAGCCACCGGCCACAGTGATAACCTCCCTGTTaca                                                                                             >  2:92202/1‑139 (MQ=11)
                                               ggCAAATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGGTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGCTaaa                                                                                             >  2:545530/1‑141 (MQ=12)
                                               agCATATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGATCCCTTCATATTTCAGAAATATCACGATAACTCACTGTAAATCTTAGATCATGGACAACGCCTCAATGAATAACATCCTTGTTcaa                                                                                             >  2:92214/2‑141 (MQ=11)
                                                 cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATATATGTCTTAGATATTAGCCAACGGCTACAGTGATAACATCCTAGTTAAAtt                                                                                           >  2:754851/1‑141 (MQ=11)
                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGGTtat                                                                                     >  2:2533/1‑141 (MQ=17)
                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATACCATCGTTGTTAAATTGTTtat                                                                                     <  1:405336/141‑1 (MQ=17)
                                                           aaCCCATCCGTAGATTGTTTAATGATGAACGTTTACACGACGTCCTCCTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                 <  2:537633/141‑1 (MQ=11)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                               >  2:394832/1‑141 (MQ=32)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAGATATTGCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGAGATAACATCCTTTTTAACTTGTTTATATCTGa                                                                               >  1:327676/1‑141 (MQ=11)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACAACACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATGTATACTTTAGATAGTAGGCAACGGCTACAGTGATAACATCCTTGTTAACTTGTTTAGATCTGa                                                                               >  2:329954/1‑141 (MQ=11)
                                                               cAACCGTAGACCCTTGAATAATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                             <  2:468732/141‑1 (MQ=12)
                                                                      aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                      <  2:96053/141‑1 (MQ=32)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                     >  1:303520/1‑141 (MQ=32)
                                                                            ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                                >  2:773191/1‑141 (MQ=32)
                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttt                                                     >  2:581495/1‑141 (MQ=17)
                                                                                           ttACACCACGTTCCTTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  2:610387/141‑1 (MQ=17)
                                                                                           ttACACCACCTTCCCTTAACATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:28769/141‑1 (MQ=17)
                                                                                           tgacacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:14448/139‑1 (MQ=17)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:124833/141‑1 (MQ=32)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:129676/141‑1 (MQ=32)
                                                                                               cccccGTTCTCTTAGTATTTCAAATATACCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                             <  1:266322/137‑1 (MQ=1)
                                                                                                ccacGTTCCCTTAATATTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                            <  1:23867/141‑1 (MQ=17)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                         >  2:18857/1‑141 (MQ=32)
                                                                                                      tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      >  2:497075/1‑141 (MQ=32)
                                                                                                      tCCCTTAACATATCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      <  2:507431/141‑1 (MQ=17)
                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  1:632673/141‑1 (MQ=35)
                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  2:157457/141‑1 (MQ=35)
                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                          >  1:440501/1‑141 (MQ=35)
                                                                                                                   cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                         >  2:773679/1‑141 (MQ=35)
                                                                                                                     aatatatCACGATAACTCAATGGATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                       <  2:144553/140‑1 (MQ=17)
                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATAATATAAATTTAACTTTGCa                     >  2:471689/1‑141 (MQ=255)
                                                                                                                       tatatCACGATAACTCAACGGATATCGTAGATAGAAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATGCAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACTTTGCa                     <  1:202110/141‑1 (MQ=11)
                                                                                                                       taaatCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  1:353346/138‑1 (MQ=21)
                                                                                                                         tatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                   <  1:416965/141‑1 (MQ=255)
                                                                                                                          atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                  <  2:577815/141‑1 (MQ=35)
                                                                                                                              cGATAACTCAATGTATATCTTAGATAGTAGCCATCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc              >  2:520093/1‑141 (MQ=25)
                                                                                                                               gatagcgccaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcc             <  2:456042/132‑1 (MQ=1)
                                                                                                                               gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAccc             >  2:193791/1‑141 (MQ=25)
                                                                                                                                 tagcTCAGTGCATATCTGATATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGGTTGGATCTTAAATAGTTCAGACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  1:550511/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:242193/138‑1 (MQ=12)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:76936/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACCCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAccacacctt           <  2:683328/138‑9 (MQ=11)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:756098/3‑141 (MQ=12)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:640460/3‑141 (MQ=9)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:231193/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:788733/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCAAACGGCTACAGTGATAACATCCTGGTTAAATTGTTTATATCTGAAATAGATCATACACAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACTGAACCGtat         >  2:190875/2‑141 (MQ=11)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:772242/1‑141 (MQ=17)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        <  2:702947/141‑1 (MQ=18)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:171321/1‑141 (MQ=16)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:710573/1‑141 (MQ=12)
                                                                                                                                     ccAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta       <  1:66445/140‑1 (MQ=40)
                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATAATTGCTGGACGTTGGATGTCTGAATCCTTGCATACAGAAGACTCCTTTTGGTTAGAAGGATACTATGAATTCAACTTTGCAACAGAACCGtattat      >  1:777847/1‑141 (MQ=9)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:60161/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  2:101883/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:518306/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:129214/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTACTTTTATTTGAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:346955/141‑3 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     >  2:583819/1‑138 (MQ=255)
                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    >  1:451434/1‑141 (MQ=255)
                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    <  1:713354/141‑1 (MQ=255)
                                                                                                                                        aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACCGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat    <  1:65944/141‑5 (MQ=255)
                                                                                                                                          gTACATCTTAGATAGTAGCCAACCGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATACTTCATACAGAAGACTCCTTCTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta  <  2:410062/141‑12 (MQ=11)
                                                                                                                                           |                                                                                                                                            
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  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: