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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,2690CT54.6% ‑4.4 / inf 621A95T (GCA→ACA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (128/152);  new base T (162/174);  total (293/328)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 625 total aligned reads displayed.

GCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36130‑36409
                                                                                                                                           |                                                                                                                                            
tccaCCTTCGTTGAAGGTGCCTGATATGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGTGACGTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                                                                                                                             <  2:522798/138‑1 (MQ=255)
gCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                                                                                                                             <  1:601273/141‑1 (MQ=32)
  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTGGTTGCTTACGCAACCCAATATCTAATGTATGTACCTCTGCATc                                                                                                                                           >  2:660714/1‑141 (MQ=12)
  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCGGGATc                                                                                                                                           >  2:544211/1‑141 (MQ=17)
  tACCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCCCTGTATc                                                                                                                                           >  2:509353/1‑141 (MQ=255)
   gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                                                                                                                          >  1:478424/2‑141 (MQ=255)
   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                                                                                                                          <  2:80064/141‑1 (MQ=32)
   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                                                                                                                          >  1:192524/1‑141 (MQ=32)
   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                                                                                                                          <  2:487120/141‑1 (MQ=32)
     cTTCGTTGAATGTGACTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                                                                                                                        <  2:360736/141‑1 (MQ=12)
          ttGAAGGTGCCTGATCTGTAATTACCTCTTGAAGTTTACCAAATTGTTTAATGAGACGTTTGATATACGCATATGCTGAATGACTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                                                                                                                   <  2:669008/141‑1 (MQ=255)
               ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAGCAATGGCACGatat                                                                                                                              <  2:454066/141‑1 (MQ=255)
                gTGCCTGATCTGTAACTACCTCTTGAGGCTTACCAAATTGTTTAATGAGACGTCTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtata                                                                                                                             <  2:218699/141‑1 (MQ=11)
                      taTCTGTAACTACCTTTTGAGGCTTACCCAATTGTTTAATGAGATGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACTAAATTTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                                                                                                                       <  2:141570/140‑1 (MQ=11)
                       aTCTGTAATTACCTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                                                                                                                      <  1:286034/141‑1 (MQ=17)
                            tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGGTCCAtt                                                                                                                 <  2:369609/141‑1 (MQ=17)
                            tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAtt                                                                                                                 >  2:598454/1‑141 (MQ=32)
                            tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAACAAATATCTAATCTATGTCCCTCTATATCAATGGCACGATATATATAGCTCCAtt                                                                                                                 >  1:514479/1‑141 (MQ=11)
                             aactACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATATAATGTATGTCCCTCTGTATCAATGGCACGCTATACATAGATCCAttt                                                                                                                >  2:339672/4‑141 (MQ=11)
                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCCCGAGATAAATAGCCCCAtttt                                                                                                               >  1:114589/1‑141 (MQ=12)
                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATAGAAATAGCCCCATTTTCCt                                                                                                            >  1:592778/1‑141 (MQ=17)
                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAGATAAATAGCTCCCTTTTCCt                                                                                                            >  1:438095/1‑141 (MQ=17)
                                  ccTTTTGCGGCTTACCAAATTTTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt                                                                                                           <  1:641120/141‑1 (MQ=255)
                                   cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                          <  1:462654/141‑1 (MQ=17)
                                       tGAGGTTTGCCAAATTGTTTAATGAGACGTTTCATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGAAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt                                                                                                      <  2:164594/141‑1 (MQ=12)
                                         aGGCTTACCAAATTGTTTAATGATACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt                                                                                                    <  1:155730/141‑1 (MQ=12)
                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGAATCATTGGCACGATATAAGTAGCTTCATTTTCCTttttttttt                                                                                                   >  2:246325/1‑140 (MQ=11)
                                            tttACCAAATTGTTTAGTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                 >  1:342761/1‑141 (MQ=17)
                                            tttACCAAATTGTGTAACGAGCCGGTCAATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATAAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                 <  2:134616/141‑1 (MQ=11)
                                             ttACAAAATTGTTTAATGAGGCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCATGGCACGATATAAATAGCTCCATTTTCGTTTTATTTGGATg                                                                                                <  2:323821/141‑1 (MQ=11)
                                               aCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACTATATAAATAGCTCCATTTTCCTTTTATTTTGAtata                                                                                              >  1:370544/1‑141 (MQ=12)
                                                ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGTTCCATTTTCCTTTTATTTTGATATAc                                                                                             >  2:634015/1‑141 (MQ=12)
                                                 cAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTATTTTTGATATAct                                                                                            >  1:176435/1‑140 (MQ=255)
                                                    aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtct                                                                                         >  1:430949/1‑141 (MQ=17)
                                                        tttAATTAGACGTTTTATTAACGCTTATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatc                                                                                     <  1:508621/141‑1 (MQ=255)
                                                        tttAATGAGACGTTGGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcatc                                                                                     <  1:382426/141‑1 (MQ=12)
                                                          tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGCATCAATGGCAAGATACAAAAAGATCCATTTTCCTTTTATTTTGATGTACGtctctttca                                                                                   >  2:209792/1‑136 (MQ=11)
                                                           attGAGACGATTGATAGACGAATATGCTGAATGATTATCTCGCTGCTTAGCCAACCAAATATCTAATGTATGGCCCTCTGAATCCATGGCACGATATAAATAGCTCTATTTTCCTTTTATTTTGGTATACGtctattctgg                                                                                  <  2:278653/139‑8 (MQ=11)
                                                           aaTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                  <  1:537363/141‑1 (MQ=32)
                                                            aTGAGACGTTTGATAAACGCATATTCTGAATGATTATCTCTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                 <  2:535963/141‑1 (MQ=17)
                                                            aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                 <  2:230824/141‑1 (MQ=12)
                                                            aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAGGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                 <  1:56558/141‑1 (MQ=255)
                                                              gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACg                                                                               <  2:578755/141‑1 (MQ=17)
                                                               ggtctgttGATAACCGCATATGCTGAATGATTATCTCGTCGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                              <  2:55593/135‑1 (MQ=11)
                                                               ggaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGc                                                                              <  1:617247/140‑1 (MQ=17)
                                                               agaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                              <  2:472624/141‑1 (MQ=12)
                                                               agaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCCCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGc                                                                              <  1:24168/141‑1 (MQ=17)
                                                                  cGTTTGATGGACGCATATGCCGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCGCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAt                                                                           <  1:12557/141‑1 (MQ=11)
                                                                  cGTTTGAAAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGCTATACGTCTCTTCAATACGCCAt                                                                           >  1:416648/1‑141 (MQ=11)
                                                                     ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                        >  2:29547/1‑141 (MQ=17)
                                                                     ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATAGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                        >  2:433414/1‑141 (MQ=255)
                                                                      tGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                       <  1:377793/141‑1 (MQ=17)
                                                                       gctAAACGCATATTCTGAATGATTATCTCGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGta                                                                      <  1:277912/139‑1 (MQ=11)
                                                                       gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATCCGCCCTTTGta                                                                      >  2:118441/1‑141 (MQ=17)
                                                                          aaccGCATATGCTTAATGATTATCTCGTTGATTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaata                                                                   <  2:443461/138‑1 (MQ=12)
                                                                          aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTATCATCAAAAGACCATTTGtaata                                                                   >  2:602946/1‑141 (MQ=11)
                                                                          aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaata                                                                   <  1:573308/141‑1 (MQ=32)
                                                                             cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAGATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGc                                                                >  2:609051/1‑141 (MQ=255)
                                                                             cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGATGTACCGATCCTCAATACGCATTTTGtactaacc                                                                >  1:98346/1‑139 (MQ=11)
                                                                                  atGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCTTCGATACGCCATTTGTAATAAGGttttt                                                           >  1:187791/1‑141 (MQ=11)
                                                                                    gCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTACTTTTATTTTGATATACGTCTAATCTATACGCCATTTGTAATGAACttttttt                                                         >  1:137670/1‑140 (MQ=11)
                                                                                      tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                       <  1:660218/141‑1 (MQ=255)
                                                                                      tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                       <  1:409752/141‑1 (MQ=38)
                                                                                       gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGCTATAAATAGCTCCATTTTCCTTTTATTTTGATATCCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                      >  2:378813/1‑141 (MQ=255)
                                                                                       caaacaTTACCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                      <  1:310103/136‑1 (MQ=14)
                                                                                          tGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTTTCATCAATACGCCATTTGTGTTAAGCTTTTTTATGtttt                                                   >  1:92445/1‑141 (MQ=14)
                                                                                          tGAATATCTCGTTGCTTACGCAACCAAATATCTAATGTACGCCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTCCGTCTCATCAATACGCCATTTGTAATACGCTCTTTTCTGCttt                                                   <  2:468851/141‑1 (MQ=11)
                                                                                           gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                  >  2:255970/1‑141 (MQ=39)
                                                                                            aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATCCGTCTCCTCAATATGCCATTTGTAATAAGCTTTTTTATGCttttt                                                 >  2:154832/1‑141 (MQ=14)
                                                                                             ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATGTAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTGATGCattttt                                                >  2:193219/1‑140 (MQ=11)
                                                                                             ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATATGTCTCATCAATACGCCATTTGTAATATGTTTTTTTATGCtttttt                                                >  1:400293/1‑140 (MQ=14)
                                                                                             ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCATTACGCCATTTGTAATAAGCTTTTTTATGCtttttt                                                >  2:501789/1‑140 (MQ=255)
                                                                                               aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGACCTTTGGACTAAGATGTTTTCTGATTTttcct                                              >  1:229346/1‑139 (MQ=11)
                                                                                                ttccgtagattACCCAATCAAATGTCTAATGTATGGACCTCTTCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                             <  1:16377/132‑1 (MQ=14)
                                                                                                tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                             <  1:331138/141‑1 (MQ=255)
                                                                                                tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                             <  1:256932/141‑1 (MQ=40)
                                                                                                 ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                            <  1:66356/141‑1 (MQ=255)
                                                                                                 cgcGTTGCTTACGCAACCAAATTTCTATTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                            <  1:20077/139‑1 (MQ=255)
                                                                                                   cGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGTTTTTTTATGCTTTTTCTTCCaa                                          >  2:269732/1‑141 (MQ=255)
                                                                                                     ttGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAATTTGCTCCATTTTCCTTTTTTTTTGATATAAGTCTCATCACTACGCCATTTGTTATACGCTTTTTTATGCTTTTACTTCCCAAt                                        >  1:55376/1‑141 (MQ=14)
                                                                                                      tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                       <  1:27030/141‑1 (MQ=39)
                                                                                                      tGCTTACGCAACCAAATATCTAATGTATGTCCCGCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                       <  2:407847/141‑1 (MQ=255)
                                                                                                           aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATCTAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACTCTATTTGTAATCAGCTTTGTTCTGCTTTTTCTTCCACATTTGaca                                  >  1:218062/1‑139 (MQ=11)
                                                                                                           aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGAtt                                  >  1:300690/1‑140 (MQ=39)
                                                                                                            cGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  1:100285/141‑2 (MQ=255)
                                                                                                            cGCAACCAAATATCTAAGGTATTTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTATACGTCTCATCAATACGCCATTTGTAATACGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  2:648119/141‑2 (MQ=11)
                                                                                                              cAACCAAATTGCTAATGTATGTCCCTCGGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaa                               <  2:311242/141‑1 (MQ=11)
                                                                                                              cAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATataa                               <  1:457663/141‑4 (MQ=255)
                                                                                                                  cAAATATCTAATGTATGTCCCTCTGCATCAATGGCCCGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTACAAATTTGATATAAAAtt                           >  1:295109/1‑141 (MQ=255)
                                                                                                                     atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTGTTTCTACCAAATTTGATATAAAATTggg                        >  1:493128/1‑141 (MQ=11)
                                                                                                                       atCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACACCATTTGTAATAGGCTTTTTTATGCATTTTCTTCACAATTTGATAAACAATTgggcc                      >  1:74285/1‑139 (MQ=14)
                                                                                                                        tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAATCTTTTTTATCATTTTTTTTCCAAATTTTCTATAAAATTtggcaa                     >  1:316804/1‑138 (MQ=11)
                                                                                                                        tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTGCGAATTTGATATAAAATTggggga                     >  2:399491/1‑139 (MQ=11)
                                                                                                                          tAATGTATTTCCTTTTGCATCAATTTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCata                   <  1:63181/141‑1 (MQ=14)
                                                                                                                            atgtatgtCACTCTGGATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAtt                 <  1:583615/141‑1 (MQ=255)
                                                                                                                              gtatgtCCCTCTGTATCAATGGCACGTTATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt               <  2:416134/141‑1 (MQ=255)
                                                                                                                                       tctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTAACATAAAATTGGGGCATattattctacctaa      >  1:525168/1‑133 (MQ=11)
                                                                                                                                          gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACgg   <  1:577629/141‑1 (MQ=39)
                                                                                                                                           cATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCCATTTTGCTATACATTTAGGGCTTAATCTTGAACCCAACgat  >  1:356831/1‑139 (MQ=11)
                                                                                                                                           cATCAATGGCACGATATAAATAGCTCCATTATCCTTTTATNTTGATGTACATCTCATCAATACGCCATATGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATNTGACATAAAATTCAGCCATATCCTTGACCCCCACgtc  >  2:351250/1‑139 (MQ=11)
                                                                                                                                           |                                                                                                                                            
GCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36130‑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: