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

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTAT  >  CP000730/36340‑36614
                                                                                                                                           |                                                                                                                                        
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGt                                                                                                                                          >  2:252146/1‑140 (MQ=255)
  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTAATCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCCTTGtat                                                                                                                                       >  1:303678/1‑140 (MQ=255)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATACATCACGATAACTCAATGTATATc                                                                                                                                    <  1:628771/140‑1 (MQ=38)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGAGCATATTCTTAAACCAAACGGTAGACAGTTGAATGATGAACGTTTACACAACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATc                                                                                                                                    >  1:182491/1‑140 (MQ=11)
       cTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCtt                                                                                                                                  <  2:476748/140‑1 (MQ=37)
         tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAg                                                                                                                                >  2:51894/1‑140 (MQ=34)
            ttATGCTTTTTCTTCCAAATTTGATCCAAAATTGGGGCATATTCTTTAACCCAACGGTAAAACGTTGAATGATGAAAGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAACTTAGAta                                                                                                                             >  1:217367/1‑140 (MQ=11)
            ttATGCTTTTTCTTCCAAATTTGATACATAATTGGGGCATATTCTTGTTCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                             <  1:184666/140‑1 (MQ=255)
             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                            >  1:389170/1‑140 (MQ=34)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                           <  2:105841/140‑1 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                           <  1:68217/140‑1 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                           <  1:163321/140‑1 (MQ=35)
              aTGCTTCTTCTTCCAAATTTGAGACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAGGTATATCTTAGAtagt                                                                                                                           <  1:438547/140‑1 (MQ=11)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTCATATTTCAGATATATCACGATAACCCACTGTCTATCTTAGAGAGTAGc                                                                                                                        >  2:519474/1‑140 (MQ=11)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATATTAGc                                                                                                                        >  2:28334/1‑140 (MQ=21)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                        >  2:28385/1‑140 (MQ=34)
                 cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                        <  1:541309/140‑1 (MQ=35)
                        ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                 <  2:609109/140‑1 (MQ=37)
                            aaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACGGCTACa                                                                                                             >  2:521121/1‑140 (MQ=25)
                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCTTAATATTTCAGATATATCACGATAACTCAATGTAAATCTTAGATAGTAGCCAACGGCTACAGTGag                                                                                                        >  1:311526/1‑139 (MQ=255)
                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                    <  2:524853/140‑1 (MQ=34)
                                      cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                                                                   >  2:118336/1‑140 (MQ=33)
                                        aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                 >  1:14435/1‑140 (MQ=32)
                                        aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAAACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                 >  1:14716/1‑140 (MQ=17)
                                        aaaTTGGGGCAAATTCTTGAACTCAACGTGAAACCGATGAATGATGAACGTTTACACCACGTTCGCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGACTACAGTGATCACCTcc                                                                                                 >  2:659706/1‑140 (MQ=11)
                                           ttGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  2:630548/140‑1 (MQ=18)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  1:551904/140‑1 (MQ=33)
                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                              <  1:415589/140‑1 (MQ=33)
                                            tGGGGCATATTCTTGAACCCAACGGTAGACGGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                             <  1:40507/140‑1 (MQ=16)
                                             ggggCATCTTGGTGAACCCAACGGTAGACCGTTGCACTGGGACCGTTTAAACGCCCACCTATAAGTATATCAAATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACGGTGATAACATCCTTGtt                                                                                            >  2:140449/1‑140 (MQ=2)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACGCAATGTAGATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                           <  2:283526/140‑1 (MQ=17)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:109408/1‑140 (MQ=32)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:552538/1‑140 (MQ=33)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATACTTGTTaa                                                                                          >  1:73872/1‑140 (MQ=18)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                          >  1:207991/1‑140 (MQ=17)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTaa                                                                                          >  1:37007/1‑140 (MQ=17)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTTAATCTTGAATAGTAGCCAACGGCTACAGTGGTAAAATCCTTTTTaaa                                                                                         >  2:86028/1‑140 (MQ=11)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                         >  2:645068/1‑140 (MQ=32)
                                                 cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTAATAAAATACCTGTTAAAt                                                                                        >  1:565365/1‑140 (MQ=11)
                                                  atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAtt                                                                                       >  2:404552/1‑140 (MQ=33)
                                                    atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                     <  1:252146/140‑1 (MQ=32)
                                                            aCCCAACGGTAGACCGTTGAATGACGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                             <  2:204354/140‑1 (MQ=17)
                                                                aaCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                         >  2:156500/1‑140 (MQ=32)
                                                                aaCGCTAGTCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                         <  2:128669/140‑1 (MQ=18)
                                                                   ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAg                                                                      >  1:397917/1‑140 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                     >  2:275542/1‑140 (MQ=33)
                                                                      ctaCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                   >  1:413571/3‑140 (MQ=17)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                  >  1:530069/1‑140 (MQ=32)
                                                                           gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATATGAAATAGTACATAcc                                                              >  1:84679/1‑139 (MQ=12)
                                                                              gAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                           <  1:608582/140‑1 (MQ=33)
                                                                               aaTGATTAATGTTTACACCATGTTTTCTTAATATTTCAAATATATCACGATAGCTCAATGTGTATTTTAGTTAGTAGCCAACGGCTACAGTGATAACAGCCTTGCTCTCTTTATAACATCGGAAATAGTTCATACagaag                                                          <  1:159411/140‑1 (MQ=2)
                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                          >  1:535945/1‑140 (MQ=32)
                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                          >  1:536092/1‑140 (MQ=32)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGATAACAACCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                        >  1:284007/1‑140 (MQ=17)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTCTATCTGAAATAGTTCATACAGAAGAc                                                        <  2:413571/140‑1 (MQ=17)
                                                                                 tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                        <  1:275542/140‑1 (MQ=33)
                                                                                    tgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTcc                                                     <  1:118336/140‑1 (MQ=32)
                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                             <  2:246545/140‑1 (MQ=33)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACTAGAGACGCCTTTTTGTTa                                            >  1:582957/1‑140 (MQ=255)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                            >  2:491271/1‑140 (MQ=17)
                                                                                                 cacGTTGTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                        <  1:261050/140‑1 (MQ=12)
                                                                                                   ctttCTCTTAATTATTCAGTGATAGCACGATAGCTCATGTCATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                      <  2:582835/138‑1 (MQ=11)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                      <  2:534056/140‑1 (MQ=33)
                                                                                                            aaTATTTCATATATATCAAGATAATGCAATGTATATCTTAGATTGTAGTTAATGTCTTCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                             <  1:18589/140‑1 (MQ=11)
                                                                                                            aaTATTTCAGAGATATCATGATAACTCAATGTATATCTTAGATAGTAGTCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                             <  1:121266/140‑1 (MQ=255)
                                                                                                             atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                            <  2:562477/140‑1 (MQ=34)
                                                                                                                      atatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                   <  2:171479/140‑1 (MQ=35)
                                                                                                                        atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                 <  2:45790/140‑1 (MQ=35)
                                                                                                                          atCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc               <  1:96730/140‑1 (MQ=21)
                                                                                                                             aCGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCCTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa            <  2:94054/140‑1 (MQ=255)
                                                                                                                               gcaaaCACAATGCATATCTTAGAAAGTAGCCAAACGCTACAGCGATACCATCCTTGTTAAATTCTTTATATCTGCAATAGTTCATACAGATGACTCCTTTATGTTAAAATTATACTATTAATTCAACTTTGCAACAgccc          <  2:297821/137‑4 (MQ=255)
                                                                                                                                taaaCTCTATGTATATCTGAGTTATTAGCCGACGTCGACCGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCCACTTTGCAACAGAAcc         <  2:313885/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATTATTATATTATAAATTCAACTTTTCAACAGAActc        >  2:95744/4‑138 (MQ=11)
                                                                                                                                  agcTCAATGCATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATATAGTTCATACAGAAGTTTCCTTTTTGTTAAGATTATATTATAAATTCAACTTTGCAACAGAACCGt       <  1:510464/138‑1 (MQ=11)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATTTTTTTATTATAAAATCAGCTTTCCAACACAACCGt       >  2:250113/3‑140 (MQ=11)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTAGTTAAATAGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt       >  1:201148/3‑140 (MQ=11)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGAAACAGAACCGta      >  1:176325/2‑140 (MQ=12)
                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGACCCGta      >  2:345241/1‑140 (MQ=38)
                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGta      >  2:335032/1‑140 (MQ=39)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  2:133445/1‑140 (MQ=255)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtat     >  2:409100/1‑140 (MQ=255)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  2:407495/1‑140 (MQ=16)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat     >  1:165123/1‑140 (MQ=12)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATATTATAAATTCAAATTTGCAACAGAACCGtat     >  1:652670/1‑140 (MQ=255)
                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta   >  1:217589/1‑140 (MQ=255)
                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATCTCTGAAATAGTTCATACAGAAGAATCCTTCCGATTAAAATAATATTATACATGCAACATTGCAACAGAACCGtattat  >  1:313865/1‑140 (MQ=11)
                                                                                                                                       aaTGGATATCTTAGATAGTAGCCAACGGCTACAGTGATAAGATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAATAGAACCGtatttt  >  1:476929/1‑138 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt  <  1:625514/140‑3 (MQ=255)
                                                                                                                                           |                                                                                                                                        
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTAT  >  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: