Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,113,304 G→T 7.4% H37N (CAT→AAT)  USA300HOU_RS10725 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,113,3040GT 7.4% 693.7 / 7.5 236H37N (CAT→AAT) USA300HOU_RS10725hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (118/95);  new base T (21/2);  total (139/97)
Fisher's exact test for biased strand distribution p-value = 5.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 240 total aligned reads displayed.

AGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCTCTCT  >  USA300TCH1516_ALE/2113164‑2113443
                                                                                                                                             |                                                                                                                                           
aGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCATCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATg                                                                                                                                             <  1:783211/141‑1 (MQ=255)
   aTGTTAGTCACTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTc                                                                                                                                          <  2:869460/141‑1 (MQ=255)
      ttatctaTTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGAt                                                                                                                                       <  2:379366/135‑1 (MQ=255)
      ttAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGAt                                                                                                                                       <  2:795008/141‑1 (MQ=255)
        aGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGt                                                                                                                                     >  1:966268/1‑141 (MQ=255)
         gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTCTTTTTCGATGTa                                                                                                                                    >  2:881674/1‑141 (MQ=255)
         gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTa                                                                                                                                    >  2:635126/1‑141 (MQ=255)
         gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTa                                                                                                                                    >  1:723595/1‑141 (MQ=255)
           cagtttCCTGTTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGACTAACTTTAAAACATGATACGGGGGCATTATTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                                  >  2:499544/4‑141 (MQ=255)
           cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                                  <  2:492377/141‑1 (MQ=255)
           cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAAAGCTACAAAGACTTTTATATATTCCTACCATTGACACGAATACCTTAAAATCATGTAATTCTGTTATTATTGTTTTTTCGTCATGGTCGGGGTaaa                                                                                                                                  >  1:808752/1‑141 (MQ=255)
            aTTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATAATGGAATTCTGTAATTATTGTTTTTTCGTCCTGTTCGATGTAAAt                                                                                                                                 >  2:288322/1‑141 (MQ=255)
                 ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTGCATAATTTTATCTTTTAAATCTTCTGACACGCATACCTTTAAATCATGGAATTCTGTATTTATTGTTTTTTGGTCTTGGTTGATGTAAATTTCTa                                                                                                                            >  2:155239/1‑141 (MQ=255)
                 ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTa                                                                                                                            >  1:630376/1‑141 (MQ=255)
                  cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTACTTATTGTTTTTTCGTCATGGTCGATGTGAATCTATAg                                                                                                                           >  1:135740/1‑141 (MQ=255)
                  cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                           >  2:749347/1‑141 (MQ=255)
                     ttcaccctcATCAATTCCAGCCAATTCATACATAACTCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGt                                                                                                                        <  2:740739/136‑1 (MQ=255)
                     ggcctcctCATCAAGTCCAGCCAATTCATACATACCACCGTCCATACTATGTTTTGTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGt                                                                                                                        <  2:387360/139‑1 (MQ=255)
                      tcctccTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAGATTTCTAGCGTa                                                                                                                       <  2:534079/141‑1 (MQ=255)
                      gcctcctCATCCATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTa                                                                                                                       <  1:192591/140‑1 (MQ=255)
                              atcaATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGc                                                                                                               >  1:435802/1‑141 (MQ=255)
                                caATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTa                                                                                                             <  1:1074977/141‑1 (MQ=255)
                                caATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTa                                                                                                             >  1:456933/1‑141 (MQ=255)
                                 aaTTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTaa                                                                                                            >  1:344442/1‑141 (MQ=255)
                                      cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                       >  1:327945/1‑141 (MQ=255)
                                      cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                       <  2:1095287/141‑1 (MQ=255)
                                       cgCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaa                                                                                                      <  1:118236/140‑1 (MQ=255)
                                        gCCAATTCATACATCACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaaa                                                                                                     <  1:906568/141‑1 (MQ=255)
                                        gCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTTGCGTACCATCAGCTAAATTaaaa                                                                                                     >  1:430957/1‑141 (MQ=255)
                                        gCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCGACTAAATTaaaa                                                                                                     >  1:653207/1‑141 (MQ=255)
                                             ttCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTACTTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTATATTCAAAATAga                                                                                                >  2:865371/1‑140 (MQ=255)
                                                atacataACTCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGCTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCttt                                                                                             <  2:279884/141‑1 (MQ=255)
                                                  acataACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTAt                                                                                           >  1:70807/1‑141 (MQ=255)
                                                           tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTCTCTTCActtt                                                                                  >  2:84076/1‑141 (MQ=255)
                                                           tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAtttt                                                                                  >  1:876404/1‑141 (MQ=255)
                                                           tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAtttt                                                                                  >  2:972846/1‑141 (MQ=255)
                                                           tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCGTTATCTTCTtttt                                                                                  >  1:297138/1‑141 (MQ=255)
                                                            ctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAACTCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTc                                                                                 <  2:130983/141‑1 (MQ=255)
                                                            ctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTTTCTGTAAATTTCTAGCGTACCATCAGTTAACTTAGAAATAGCTTTATCTCCGttttt                                                                                 >  2:781958/1‑140 (MQ=255)
                                                            ctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTCTAAATCATGGAATTCTGTACTTATTGTTTTTTCGTCATGGTCGATGTCAATTTCTAGCGTACCATCAGCTAACTTAAACATAGCATTCTCTTCCTTTTc                                                                                 >  2:203019/1‑141 (MQ=255)
                                                              ccATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTTGATGTAAATTTCTGGCGTACCGTCAGCTAAATTAAAATTAGCTTTATTTTAATTTTCtt                                                                               >  2:509553/1‑141 (MQ=255)
                                                              ccATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCtt                                                                               <  2:819638/141‑1 (MQ=255)
                                                               cATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCttt                                                                              >  1:931966/1‑141 (MQ=255)
                                                                aTACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCtttt                                                                             >  1:865033/1‑140 (MQ=255)
                                                                aTACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCtttt                                                                             >  2:698115/1‑140 (MQ=255)
                                                                      ttATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTCTGTATCCt                                                                       >  1:560420/1‑141 (MQ=255)
                                                                         tCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATTTTCATTTTCTTTGTATCCtttt                                                                    >  2:706731/1‑141 (MQ=255)
                                                                            tttACATCTTCTGACACGAATACCTTTAAATCATGCAACTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCAGCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCt                                                                 <  1:94016/141‑1 (MQ=255)
                                                                             ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtt                                                                <  1:384795/141‑1 (MQ=255)
                                                                             ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtt                                                                >  1:200614/1‑141 (MQ=255)
                                                                             ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGGCCTGGTCGATGTAAATTTCTCGCGTACCATCAGCTACATTAAAAATAGCTTTATCTTCATTTTCTTTGTTTCCGTTTTCtt                                                                >  1:58306/1‑141 (MQ=255)
                                                                              tAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCttt                                                               <  1:880320/141‑1 (MQ=255)
                                                                              tAAATCTTCTGACACGAATACCTTTAAATCAGGGAATTCTAATTCACCTACTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCttt                                                               <  1:520905/141‑1 (MQ=255)
                                                                               aaaTCTTCTGACACGGATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTGTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                              <  1:900801/141‑1 (MQ=255)
                                                                               aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTTTTTGTAAATTTCTAGCTTACCATCAGCTAAATTAGACTTAGCTTTTCCTTCGTTTTCTTTGTATCTtttttttttt                                                              >  2:613354/1‑141 (MQ=255)
                                                                               aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                              <  1:48998/141‑1 (MQ=255)
                                                                               aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGAACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                              <  1:664744/141‑1 (MQ=255)
                                                                                  tcttctGACACGGATACCTTTTCATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAAt                                                           <  2:302517/141‑1 (MQ=255)
                                                                                    ttctGACCCGAATAGCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCTTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAAttt                                                         <  2:191648/141‑1 (MQ=255)
                                                                                     tctGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAttttt                                                        >  2:554383/1‑140 (MQ=255)
                                                                                     tctGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCTTTTTTTTTTATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCCTTTTCTTTGTATCCTTTTTTTTTTTAtttt                                                        >  2:939743/1‑140 (MQ=255)
                                                                                      ctGACACGAATACCTCTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATttct                                                       >  1:553765/1‑141 (MQ=255)
                                                                                         acacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTACTTTCTTCa                                                    >  1:83238/1‑141 (MQ=255)
                                                                                         acacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCa                                                    >  1:991390/1‑141 (MQ=255)
                                                                                          cacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAGTTAAAAATAGCTTTATCTTCATTTTCTTTGTTTCCTTTTTCTTTTATTTTGTTCaa                                                   >  2:107887/1‑141 (MQ=255)
                                                                                          cacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCaa                                                   >  2:943791/1‑141 (MQ=255)
                                                                                            cGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTAGCCTTTTTCTTTTAATTGCTTCaaaa                                                 >  1:499792/1‑140 (MQ=255)
                                                                                                  ccTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTTTTTTTTTTTCTTCAATGATTcct                                           >  2:985673/1‑139 (MQ=255)
                                                                                                    tttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTa                                         >  2:1040480/1‑141 (MQ=255)
                                                                                                     ttAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTTGATGTAAATTTCTAGCGTACCATCAGCTAAGTTAAAAATAGCTTTATGTTCGTGTTCTTTGTATCCTTTTTCTTTTATTTTATTCAGAGTTttattat                                        >  2:1013671/1‑140 (MQ=255)
                                                                                                     ttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTaa                                        <  1:303081/141‑1 (MQ=255)
                                                                                                     ttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAATACCGTTGTATCCTTTTTTTTTTAATTTTTTCAAAGATTCATTaa                                        >  2:433739/1‑141 (MQ=255)
                                                                                                      tAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAat                                       >  2:399721/1‑141 (MQ=255)
                                                                                                        aaTCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattt                                     <  2:456933/141‑1 (MQ=255)
                                                                                                         aTCATGCAATTCTGTAATTTTTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAatttt                                    <  2:149945/141‑1 (MQ=255)
                                                                                                            aTGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAGGATTCATTAtttttcat                                 >  1:42127/1‑141 (MQ=255)
                                                                                                             tGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                                >  2:506039/1‑141 (MQ=255)
                                                                                                             tGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                                >  1:142988/1‑141 (MQ=255)
                                                                                                                 aTTCTGTAATTATTGTTTTTTTTTTTTTTTTGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAATGATTTATTTatttttattttct                            >  2:677443/1‑140 (MQ=255)
                                                                                                                 aTTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTcc                            >  1:252198/1‑141 (MQ=255)
                                                                                                                 aTTCTGGAATTATTGTTTTTTCTTCTTGTTCGAAGTAAATTTCTAGTGTACCATCAGCTCAATTAAAAATAGCTTTCTTTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCCTTATTTTTCATTTTcc                            <  2:898571/141‑1 (MQ=255)
                                                                                                                  ttCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCt                           <  1:609144/141‑1 (MQ=255)
                                                                                                                   tCTTTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                          <  1:323777/141‑1 (MQ=255)
                                                                                                                   tCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAGATAGCTTTATCTTCTTTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                          <  1:479517/141‑1 (MQ=255)
                                                                                                                    cTGTAATTATTGTTTTTTCGTCATGGTCGATGTACATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGt                         <  2:597133/141‑1 (MQ=255)
                                                                                                                    cTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTAATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGt                         <  1:793615/141‑1 (MQ=255)
                                                                                                                       tAATTATTGTTTTTTTTTTATTTTTTATGTAAATTACTAGCGTACCATCAGCTAAATTAGAATTAGCTTTATCTTCATTTTCTTTGTATCCTTTTTTTTTTAATTTACTCAATGATTCATTACTTTTCATTTTCCTGTtcc                      >  2:618784/1‑141 (MQ=255)
                                                                                                                       tAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTtcc                      >  2:652386/1‑141 (MQ=255)
                                                                                                                         attattGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCCTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTtcctc                    <  2:327945/141‑1 (MQ=255)
                                                                                                                         attactGTTGTCTCGTCATGGTCGATGTAAATTTCTAGCGTACCATAAAATAAATAAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTtcctc                    <  2:365996/141‑1 (MQ=255)
                                                                                                                           tattGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTcctcct                  <  1:881674/141‑1 (MQ=255)
                                                                                                                           tattGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTcctcct                  <  1:635126/141‑1 (MQ=255)
                                                                                                                                 ttttttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCatat             <  1:1093875/141‑1 (MQ=255)
                                                                                                                                  tttttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCGTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATAtt            <  2:945380/141‑1 (MQ=255)
                                                                                                                                    tttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCa          <  1:396640/141‑1 (MQ=255)
                                                                                                                                        gTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGct      <  2:482974/141‑1 (MQ=255)
                                                                                                                                         tCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGctc     >  1:592790/1‑141 (MQ=255)
                                                                                                                                         tCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTCTCCCGTCCGTCCTCACATTCAAGctc     >  2:1056183/1‑141 (MQ=255)
                                                                                                                                          cATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCATCATATTCAAActct    >  1:670581/1‑141 (MQ=255)
                                                                                                                                           aTGGTCGATGTCAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTACCTTCATTTTCTTTGTATCCTTTTACTTGTGGTTTCTTCTAAGATTCATTATTTTTCATTTTCCTGTTCCTCCTCATATTCAAGctctc   <  2:176858/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCTCTCT  >  USA300TCH1516_ALE/2113164‑2113443

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