Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 15.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,4790TC15.2% 156.0 / 24.8 86T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (9/4);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTAT  >  CP000730/36340‑36614
                                                                                                                                           |                                                                                                                                        
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGt                                                                                                                                          >  2:252146/1‑140 (MQ=255)
  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTAATCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCCTTGt‑at                                                                                                                                       >  1:303678/1‑140 (MQ=255)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATACATCACGATAACTCAATGT‑ATATc                                                                                                                                    <  1:628771/140‑1 (MQ=38)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGAGCATATTCTTAAACCAAACGGTAGACAGTTGAATGATGAACGTTTACACAACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATc                                                                                                                                    >  1:182491/1‑140 (MQ=11)
       cTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCtt                                                                                                                                  <  2:476748/140‑1 (MQ=37)
         tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAg                                                                                                                                >  2:51894/1‑140 (MQ=34)
            ttATGCTTTTTCTTCCAAATTTGATCCAAAATTGGGGCATATTCTTTAACCCAACGGTAAAACGTTGAATGATGAAAGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATAACTTAGAta                                                                                                                             >  1:217367/1‑140 (MQ=11)
            ttATGCTTTTTCTTCCAAATTTGATACATAATTGGGGCATATTCTTGTTCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAta                                                                                                                             <  1:184666/140‑1 (MQ=255)
             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtag                                                                                                                            >  1:389170/1‑140 (MQ=34)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                                                           <  2:105841/140‑1 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                                                           <  1:68217/140‑1 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                                                           <  1:163321/140‑1 (MQ=35)
              aTGCTTCTTCTTCCAAATTTGAGACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAGGT‑ATATCTTAGAtagt                                                                                                                           <  1:438547/140‑1 (MQ=11)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTCATATTTCAGATATATCACGATAACCCACTGT‑CTATCTTAGAGAGTAGc                                                                                                                        >  2:519474/1‑140 (MQ=11)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATATTAGc                                                                                                                        >  2:28334/1‑140 (MQ=21)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGc                                                                                                                        >  2:28385/1‑140 (MQ=34)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGc                                                                                                                        <  1:541309/140‑1 (MQ=35)
                        ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGc                                                                                                                 <  2:609109/140‑1 (MQ=37)
                            aaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAAATAGTAGCCAACGGCTACa                                                                                                             >  2:521121/1‑140 (MQ=25)
                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑AAATCTTAGATAGTAGCCAACGGCTACAGTGag                                                                                                        >  1:311526/1‑139 (MQ=255)
                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                    <  2:524853/140‑1 (MQ=34)
                                      cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                                                                   >  2:118336/1‑140 (MQ=33)
                                        aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                 >  1:14435/1‑140 (MQ=32)
                                        aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAAACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                 >  1:14716/1‑140 (MQ=17)
                                        aaaTTGGGGCAAATTCTTGAACTCAACGTGAAACCGATGAATGATGAACGTTTACACCACGTTCGCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGACTACAGTGATCACCTcc                                                                                                 >  2:659706/1‑140 (MQ=11)
                                           ttGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  2:630548/140‑1 (MQ=18)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  1:551904/140‑1 (MQ=33)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  1:415589/140‑1 (MQ=33)
                                            tGGGGCATATTCTTGAACCCAACGGTAGACGGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                             <  1:40507/140‑1 (MQ=16)
                                             ggggCATCTTGGTGAACCCAACGGTAGACCGTTGCACTGGGACCGTTTAAACGCCCACCTATAAGTATATCAAATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACGGTGATAACATCCTTGtt                                                                                            >  2:140449/1‑140 (MQ=2)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACGCAATGT‑AGATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                           <  2:283526/140‑1 (MQ=17)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:109408/1‑140 (MQ=32)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:552538/1‑140 (MQ=33)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATACTTGTTaa                                                                                          >  1:73872/1‑140 (MQ=18)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:207991/1‑140 (MQ=17)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTaa                                                                                          >  1:37007/1‑140 (MQ=17)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑TAATCTTGAATAGTAGCCAACGGCTACAGTGGTAAAATCCTTTTTaaa                                                                                         >  2:86028/1‑140 (MQ=11)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                         >  2:645068/1‑140 (MQ=32)
                                                 cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTAATAAAATACCTGTTAAAt                                                                                        >  1:565365/1‑140 (MQ=11)
                                                  atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAtt                                                                                       >  2:404552/1‑140 (MQ=33)
                                                    atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                     <  1:252146/140‑1 (MQ=32)
                                                            aCCCAACGGTAGACCGTTGAATGACGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                             <  2:204354/140‑1 (MQ=17)
                                                                aaCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                         >  2:156500/1‑140 (MQ=32)
                                                                aaCGCTAGTCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                         <  2:128669/140‑1 (MQ=18)
                                                                   ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAg                                                                      >  1:397917/1‑140 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                     >  2:275542/1‑140 (MQ=33)
                                                                      ctaCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                   >  1:413571/3‑140 (MQ=17)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                  >  1:530069/1‑140 (MQ=32)
                                                                           gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAAATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATATGAAATAGTACATAcc                                                              >  1:84679/1‑139 (MQ=12)
                                                                              gAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                           <  1:608582/140‑1 (MQ=33)
                                                                               aaTGATTAATGTTTACACCATGTTTTCTTAATATTTCAAATATATCACGATAGCTCAATGT‑GTATTTTAGTTAGTAGCCAACGGCTACAGTGATAACAGCCTTGCTCTCTTTATAACATCGGAAATAGTTCATACagaag                                                          <  1:159411/140‑1 (MQ=2)
                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                          >  1:535945/1‑140 (MQ=32)
                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                          >  1:536092/1‑140 (MQ=32)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCCACGGCTACAGTGATAACAACCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                        >  1:284007/1‑140 (MQ=17)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTCTATCTGAAATAGTTCATACAGAAGAc                                                        <  2:413571/140‑1 (MQ=17)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                        <  1:275542/140‑1 (MQ=33)
                                                                                    tgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTcc                                                     <  1:118336/140‑1 (MQ=32)
                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                             <  2:246545/140‑1 (MQ=33)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACTAGAGACGCCTTTTTGTTa                                            >  1:582957/1‑140 (MQ=255)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATAGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                            >  2:491271/1‑140 (MQ=17)
                                                                                                 cacGTTGTCTTAATATTTCAGATATATCACGATAGCTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                        <  1:261050/140‑1 (MQ=12)
                                                                                                   ctttCTCTTAATTATTCAGTGATAGCACGATAGCTC‑ATGTCATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                      <  2:582835/138‑1 (MQ=11)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                      <  2:534056/140‑1 (MQ=33)
                                                                                                            aaTATTTCATATATATCAAGATAATGCAATGT‑ATATCTTAGATTGTAGTTAATGTCTTCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                             <  1:18589/140‑1 (MQ=11)
                                                                                                            aaTATTTCAGAGATATCATGATAACTCAATGT‑ATATCTTAGATAGTAGTCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                             <  1:121266/140‑1 (MQ=255)
                                                                                                             atatTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                            <  2:562477/140‑1 (MQ=34)
                                                                                                                      atatatCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                   <  2:171479/140‑1 (MQ=35)
                                                                                                                        atatCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                 <  2:45790/140‑1 (MQ=35)
                                                                                                                          atCACGATAGCTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc               <  1:96730/140‑1 (MQ=21)
                                                                                                                             aCGATAGCTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCCTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa            <  2:94054/140‑1 (MQ=255)
                                                                                                                               gcaaaCACAATGC‑ATATCTTAGAAAGTAGCCAAACGCTACAGCGATACCATCCTTGTTAAATTCTTTATATCTGCAATAGTTCATACAGATGACTCCTTTATGTTAAAATTATACTATTAATTCAACTTTGCAACAgccc          <  2:297821/137‑4 (MQ=255)
                                                                                                                                taaaCTCTATGT‑ATATCTGAGTTATTAGCCGACGTCGACCGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCCACTTTGCAACAGAAcc         <  2:313885/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATTATTATATTATAAATTCAACTTTTCAACAGAActc        >  2:95744/4‑138 (MQ=11)
                                                                                                                                  agcTCAATGC‑ATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATATAGTTCATACAGAAGTTTCCTTTTTGTTAAGATTATATTATAAATTCAACTTTGCAACAGAACCGt       <  1:510464/138‑1 (MQ=11)
                                                                                                                                  agcTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATTTTTTTATTATAAAATCAGCTTTCCAACACAACCGt       >  2:250113/3‑140 (MQ=11)
                                                                                                                                  agcTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTAGTTAAATAGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt       >  1:201148/3‑140 (MQ=11)
                                                                                                                                   gcTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGAAACAGAACCGta      >  1:176325/2‑140 (MQ=12)
                                                                                                                                   aCTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGACCCGta      >  2:345241/1‑140 (MQ=38)
                                                                                                                                   aCTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGta      >  2:335032/1‑140 (MQ=39)
                                                                                                                                    cTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  2:133445/1‑140 (MQ=255)
                                                                                                                                    cTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtat     >  2:409100/1‑140 (MQ=255)
                                                                                                                                    cTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  2:407495/1‑140 (MQ=16)
                                                                                                                                    cTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  1:165123/1‑140 (MQ=12)
                                                                                                                                    cTCAATGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATATTATAAATTCAAATTTGCAACAGAACCGtat     >  1:652670/1‑140 (MQ=255)
                                                                                                                                      cAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta   >  1:217589/1‑140 (MQ=255)
                                                                                                                                       aaTGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATCTCTGAAATAGTTCATACAGAAGAATCCTTCCGATTAAAATAATATTATACATGCAACATTGCAACAGAACCGtattat  >  1:313865/1‑140 (MQ=11)
                                                                                                                                       aaTGG‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAAGATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAATAGAACCGtatttt  >  1:476929/1‑138 (MQ=11)
                                                                                                                                       aaTGC‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt  <  1:625514/140‑3 (MQ=255)
                                                                                                                                           |                                                                                                                                        
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTAT  >  CP000730/36340‑36614

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