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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,2690CT57.4% 80.0 / inf 346A95T (GCA→ACA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (77/71);  new base T (73/124);  total (151/195)
Fisher's exact test for biased strand distribution p-value = 6.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409
                                                                                                                                            |                                                                                                                                            
tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                                                                                                                              >  1:372603/1‑141 (MQ=32)
  ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATCTAATGTATGACCCTCTGCAt                                                                                                                                            >  2:139710/4‑141 (MQ=11)
  ccaCCTTCGTCGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCACATTGTTTAATGCGACGTTTGATAAACTCAGATGCTGAATTAGTATCTCGTTGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                                                                                                                            <  2:94132/139‑1 (MQ=11)
   tACCTTCGATGAAGGTGCCTGATCTTTAATTACCTTTTGATGTTTACCGAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                                                                                                                           <  1:110400/141‑1 (MQ=11)
     cctggtgtggaGGGGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGACACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                                                                                                                         <  2:174065/131‑1 (MQ=11)
      cTTCGTTGAAGTTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAATTTGGTTAATGATACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                                                                                                                        <  2:298450/141‑1 (MQ=11)
       ttCTTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                                                                                                                       <  2:301824/141‑1 (MQ=255)
        tCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATgg                                                                                                                                      <  2:226687/141‑1 (MQ=17)
         cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGTATATGCTGAATGATCATCTCGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                                                                                                                     <  2:390236/141‑1 (MQ=12)
         cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                                                                                                                     <  1:398934/141‑1 (MQ=32)
          gTTGAAGGTGCCTGATCTGTAACCACCTTTTGATGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCa                                                                                                                                    <  1:47687/141‑1 (MQ=11)
           ttGAAGGTGCCTCATCTAAGACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTACTATATGTACCTCTGTATCAATGGCAc                                                                                                                                   >  2:86283/1‑141 (MQ=11)
              aaGGTGCCTGATCTGTAACTACCTGTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                                                                                                                <  2:115127/141‑1 (MQ=255)
               aGGGGCCTGTTTGGTAATTACCTTTTGAGGGTTACCAAATGGTTTAAGGGGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGata                                                                                                                               <  2:41248/141‑1 (MQ=11)
                ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                                                                                                                              >  1:35830/1‑141 (MQ=32)
                 gTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAtata                                                                                                                             <  2:392516/141‑1 (MQ=32)
                    ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                                                                                                                          >  1:57507/1‑141 (MQ=12)
                        aTCTTTAACTACCTTGTGAGGCTTATCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCATCCAAATATCTAATGTATGTCCCTCTGTATCAATGTTACGATATAAATAGCt                                                                                                                      <  2:283660/141‑1 (MQ=11)
                        aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                                                                                                                      <  2:253442/141‑1 (MQ=17)
                          cttcaagtACCTTTTGAGGTTCACCAAATTGTTTAATGATCCGTTTGATAAACGCATATGCTGAATGATTAGCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                                                                                                                    <  1:13006/134‑1 (MQ=11)
                          cTGTACTTACCTTTGTATGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTACCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACCATATAAATAGCTcc                                                                                                                    <  1:213130/141‑1 (MQ=11)
                           tGTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa                                                                                                                   <  2:198305/141‑1 (MQ=17)
                             tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCACCAtt                                                                                                                 >  2:48715/1‑141 (MQ=17)
                                 tACCTTTTGATGTTTACCTAATTGTTTAATGAGACGTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTcc                                                                                                             <  1:86771/141‑1 (MQ=12)
                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt                                                                                                            >  1:152227/1‑141 (MQ=17)
                                  aCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt                                                                                                            <  2:202323/141‑1 (MQ=17)
                                   ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTAAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATGAATGGCACGATATAAATAGCTCCAATTTCCtt                                                                                                           >  1:267572/1‑141 (MQ=12)
                                   ccTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt                                                                                                           <  1:145100/141‑1 (MQ=17)
                                    cTTTTGTGGCTTACGAAATTGTTTAATGGGACGTTTGATAATCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                          <  1:134908/141‑1 (MQ=11)
                                    cTTTTGATGCTTACCAAATTGTTTAATGTGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                          <  1:100005/141‑1 (MQ=255)
                                     ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtttt                                                                                                         <  2:6420/141‑1 (MQ=17)
                                     tttcGAGCTTGACCACATTGTTGAATTGGACTTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtttt                                                                                                         <  2:375596/141‑1 (MQ=11)
                                         gAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTTATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt                                                                                                     >  1:330398/1‑141 (MQ=12)
                                          aGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATACGAATAGCACCCTCTTCCTTTTAtttt                                                                                                    >  2:292533/1‑141 (MQ=11)
                                           ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACCTCTGCATCAATTGCACGATATAAATAGCTCCATTTTCCTTTTATTTTg                                                                                                   >  1:91188/1‑141 (MQ=17)
                                             tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTTCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttttatttttat                                                                                                 >  1:346321/1‑141 (MQ=17)
                                             cttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                 <  2:259627/140‑1 (MQ=17)
                                             ctgaCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                 <  2:397932/138‑1 (MQ=17)
                                              ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAATCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATGGCTCCTTTTTCTTTTTATTttgttg                                                                                                >  1:221418/1‑141 (MQ=11)
                                                 ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATAc                                                                                             <  2:244195/141‑1 (MQ=17)
                                                 ccAAATAGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTACATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTAc                                                                                             >  2:345032/1‑141 (MQ=12)
                                                  cAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTCCg                                                                                            >  1:123660/1‑141 (MQ=17)
                                                       tgtttaattgtccGTTTGATATACGCATAGGCTGAATGATTATCTCGTCGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtca                                                                                       <  1:234116/129‑1 (MQ=11)
                                                       tGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtca                                                                                       <  1:315719/141‑1 (MQ=12)
                                                         tttAATGAGACGTTCGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcatc                                                                                     >  2:271157/1‑141 (MQ=12)
                                                          ttAATGAGACGTTTGATAAACGCACATGCTGAATGCTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcatca                                                                                    <  1:72735/141‑1 (MQ=255)
                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATGGCTCCCTTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                   >  2:332490/1‑141 (MQ=17)
                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGTCCCATTTTCCTTTTATTTTGATGTAGGTCACATCaa                                                                                   >  2:76155/1‑141 (MQ=255)
                                                            aaTGAGACGTTTGATAAACGCATCTGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCaaa                                                                                  >  1:36575/1‑140 (MQ=12)
                                                            aaTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAAt                                                                                  <  2:10472/141‑1 (MQ=17)
                                                            aaTGAGACCTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAAt                                                                                  <  1:367348/141‑1 (MQ=12)
                                                             aggagaCGTGTGAAAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                 <  1:3099/139‑1 (MQ=12)
                                                              tGAGATGTTTTATAAACCCATATGCTGCATGATTATCTCGTTGGTTACGCAACCAAATCTCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATAc                                                                                <  2:57631/141‑1 (MQ=11)
                                                              tGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCCTTTAACATGTTTTTTGATATACGTCTCATCAATAc                                                                                <  2:365776/141‑1 (MQ=11)
                                                                ggaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATTCAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                              <  2:381795/140‑1 (MQ=255)
                                                                atgcGTTTGATCAACGCAGATGCGGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATTTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                              <  1:339983/138‑1 (MQ=11)
                                                                 gaCGTTTGATAAACTCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGcc                                                                             <  1:158186/141‑1 (MQ=12)
                                                                  aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTGATCAATTCTCCa                                                                            >  1:388099/1‑141 (MQ=11)
                                                                     tttGATAAACGCATATGCTGAATGATTATGTCGTTGCTTACCAAACCAAAGATGTAATGTATGTCCCTCGGTATCAATGGCACGATATAAATAGCTACAGTTGCCTTTTATTTTGATATCAGTCTCATCAATACGCCAttt                                                                         >  2:163639/1‑141 (MQ=11)
                                                                      ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATCCGCCATTTg                                                                        >  1:321443/1‑141 (MQ=17)
                                                                        gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCAATATCTAATGTATGTCCCTCTGTATCAATGGCACGTTATAAATAGCTCCCTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAGTTGta                                                                      >  2:318755/1‑141 (MQ=11)
                                                                        gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGATATACGTCTCATCAATACGCCATTTTta                                                                      >  1:250626/1‑141 (MQ=255)
                                                                         aTAAACGCATATGCTGAATGATTATCTCGTTGCTTAGGCAACCAAATATCCAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaa                                                                     <  1:270640/141‑1 (MQ=255)
                                                                          tAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaat                                                                    <  2:50943/141‑1 (MQ=17)
                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTTCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaata                                                                   >  1:106740/1‑141 (MQ=17)
                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATGCGTCTCATCAATACGCCATTTTttata                                                                   >  1:162683/1‑141 (MQ=11)
                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCATTAAGCCATTTTtaata                                                                   >  1:186398/1‑141 (MQ=11)
                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaata                                                                   <  1:116435/141‑1 (MQ=18)
                                                                            aaCGCATATACTGAATGATTATCTAGATGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                  >  1:292971/1‑141 (MQ=11)
                                                                             aCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAg                                                                 <  1:334960/141‑1 (MQ=34)
                                                                                    tGCTGAATGATTATCTCGTTGCTTGCAACACCATATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                          >  1:193045/1‑141 (MQ=11)
                                                                                      cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTAt                                                        >  2:101154/1‑141 (MQ=37)
                                                                                          aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTAATCAATACGCCATTTGTAATAAGTTTTTTTATGAtt                                                    >  1:362779/1‑141 (MQ=14)
                                                                                             aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                 >  2:196778/1‑141 (MQ=255)
                                                                                                ttcccGTTGCTTGCGTAACCAAAGTTCTAATGTATTTCCCTGTGTATCAATGGTACGATATAAATAGCTCCATTTTCCTTGTATTTTGTGATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                              <  2:194900/137‑1 (MQ=11)
                                                                                                aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                              <  1:277827/141‑1 (MQ=255)
                                                                                                 tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTGTATGCTTTttcttc                                             <  1:209124/141‑1 (MQ=255)
                                                                                                    cGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCaa                                          <  2:191930/141‑1 (MQ=255)
                                                                                                      ttGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAt                                        <  2:152134/141‑1 (MQ=255)
                                                                                                      ttGCTCACGCAACCAAATATCTAATGCATGTCCCCCTGTATCAATGGCACGATATAAATAGCACCATTTTCCTTTTATTTTGATACCCGTCTCATCAATACGCCCTTTGTAATCACCTTCTTGATGCGTTTTCTTCCcaac                                        >  2:239559/1‑140 (MQ=11)
                                                                                                          ttACGCAACCAAATACCTAATGTATTTTCCTCTTCATCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                    <  2:1228/141‑1 (MQ=11)
                                                                                                           tACGCAACCAAATATATAATGTATGACCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCGTTTTTATGCTTTTTCTTGCAATTTTGAt                                   >  1:199677/1‑141 (MQ=11)
                                                                                                             cGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  1:377817/141‑2 (MQ=255)
                                                                                                              gCAAACAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                <  2:216560/141‑3 (MQ=255)
                                                                                                              ccAACCAAATATCTAATGTATGTCCCTCTGTATCAATGACACGAAATAAATAGCTCCATTTTCCTTTTATTTTGATCTACGTCTCATCAATACGCCGTTGTTAATAACTCTTTTGTTGCTTTTGCTTCATAAttcatgata                                >  2:267480/2‑134 (MQ=11)
                                                                                                                aaCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaa                              <  1:181301/141‑1 (MQ=255)
                                                                                                                  ccAGATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTATATTTTGATATACGTCTCATCAATAATCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACaggtg                            >  1:286375/1‑137 (MQ=14)
                                                                                                                  ccAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAt                            >  2:313882/1‑141 (MQ=38)
                                                                                                                   cAAATATCTAATGTATGTCCCTCTGTATCAATGACACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATCCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                           <  1:259232/141‑1 (MQ=11)
                                                                                                                    aaaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTg                          <  1:102167/141‑1 (MQ=255)
                                                                                                                     aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                         <  2:135787/141‑1 (MQ=255)
                                                                                                                            aaTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTATACGTCTCCTCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                  >  1:168316/1‑141 (MQ=11)
                                                                                                                             atgtttgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAtt                 <  2:337239/141‑1 (MQ=255)
                                                                                                                             atgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatact                 >  1:351913/1‑139 (MQ=255)
                                                                                                                             atgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTGGATGTAAAATTGGAGCATAtt                 >  1:280969/1‑141 (MQ=11)
                                                                                                                               gtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt               >  2:227065/1‑141 (MQ=255)
                                                                                                                                tatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCtt              <  2:170487/141‑1 (MQ=38)
                                                                                                                                           gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg   <  1:63103/141‑1 (MQ=18)
                                                                                                                                           gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACgg   <  2:106740/141‑1 (MQ=40)
                                                                                                                                            cATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGt  <  1:17110/141‑1 (MQ=18)
                                                                                                                                            |                                                                                                                                            
TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

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: