Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,113,305 G→T 7.3% D36E (GAC→GAA USA300HOU_RS10725 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,113,3050GT 7.3% 702.4 / 6.1 236D36E (GAC→GAAUSA300HOU_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.99e-01

Only 100 of 239 total aligned reads displayed.

GGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCTCTCT  >  USA300TCH1516_ALE/2113165‑2113443
                                                                                                                                             |                                                                                                                                          
ggATGTTGGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATgg                                                                                                                                            <  2:308161/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)
        gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTCTTTTTCGATGTa                                                                                                                                    >  2:881674/1‑141 (MQ=255)
        gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTa                                                                                                                                    >  2:635126/1‑141 (MQ=255)
        gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTa                                                                                                                                    >  1:723595/1‑141 (MQ=255)
          cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                                  <  2:492377/141‑1 (MQ=255)
          cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                                  >  2:1084129/1‑141 (MQ=255)
           aTTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATAATGGAATTCTGTAATTATTGTTTTTTCGTCCTGTTCGATGTAAAt                                                                                                                                 >  2:288322/1‑141 (MQ=255)
                ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTGCATAATTTTATCTTTTAAATCTTCTGACACGCATACCTTTAAATCATGGAATTCTGTATTTATTGTTTTTTGGTCTTGGTTGATGTAAATTTCTa                                                                                                                            >  2:155239/1‑141 (MQ=255)
                ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTa                                                                                                                            >  1:68963/1‑141 (MQ=255)
                 cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAAGGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATGCCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAGATTTCtgg                                                                                                                           >  2:598823/1‑139 (MQ=255)
                 cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTACTTATTGTTTTTTCGTCATGGTCGATGTGAATCTATAg                                                                                                                           >  1:135740/1‑141 (MQ=255)
                 cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                           <  1:623988/141‑1 (MQ=255)
                 cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                           >  2:749347/1‑141 (MQ=255)
                 cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                           >  1:407224/1‑141 (MQ=255)
                    ttcaccctcATCAATTCCAGCCAATTCATACATAACTCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGt                                                                                                                        <  2:740739/136‑1 (MQ=255)
                    ggcctcctCATCAAGTCCAGCCAATTCATACATACCACCGTCCATACTATGTTTTGTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGt                                                                                                                        <  2:387360/139‑1 (MQ=255)
                     tcctccTCATCTATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTa                                                                                                                       <  2:1024632/141‑1 (MQ=255)
                     tcctccTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAGATTTCTAGCGTa                                                                                                                       <  2:534079/141‑1 (MQ=255)
                             atcaATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGc                                                                                                               >  1:435802/1‑141 (MQ=255)
                               caATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTa                                                                                                             >  1:456933/1‑141 (MQ=255)
                               caATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTa                                                                                                             <  1:1074977/141‑1 (MQ=255)
                                aaTTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTaa                                                                                                            >  1:344442/1‑141 (MQ=255)
                                     cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                       <  2:1095287/141‑1 (MQ=255)
                                     cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                       >  1:327945/1‑141 (MQ=255)
                                      cgCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaa                                                                                                      <  1:118236/140‑1 (MQ=255)
                                       gCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTTGCGTACCATCAGCTAAATTaaaa                                                                                                     >  1:430957/1‑141 (MQ=255)
                                            ttCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTACTTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTATATTCAAAATAga                                                                                                >  2:865371/1‑140 (MQ=255)
                                              catacataACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCtt                                                                                              <  1:769320/141‑1 (MQ=255)
                                               atacataACTCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGCTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCttt                                                                                             <  2:279884/141‑1 (MQ=255)
                                               agacataACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTGCTAGCGTACGATCAGCTAAATTAGCAATAGCttt                                                                                             <  1:810118/139‑1 (MQ=255)
                                                     aaCGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTc                                                                                       <  1:64580/141‑1 (MQ=255)
                                                      aCGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCa                                                                                      <  1:578953/141‑1 (MQ=255)
                                                          tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTCTCTTCActtt                                                                                  >  2:84076/1‑141 (MQ=255)
                                                          tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAtttt                                                                                  >  2:972846/1‑141 (MQ=255)
                                                          tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAtttt                                                                                  >  1:872608/1‑141 (MQ=255)
                                                          tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAtttt                                                                                  >  1:876404/1‑141 (MQ=255)
                                                          tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCGTTATCTTCTtttt                                                                                  >  1:297138/1‑141 (MQ=255)
                                                           ctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTTTTTTTTTTTTCTGTAAATTTCTAGCGTACCATCAGTTAACTTAGAAATAGCTTTATCTCCGttttt                                                                                 >  2:781958/1‑140 (MQ=255)
                                                            tcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCt                                                                                >  1:19255/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                                                                             >  2:698115/1‑140 (MQ=255)
                                                                tACTTTTATCTTTTAAATATTCTGAAACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTAGTCATGTTCGCTGTAAATTTATAGCGAACCATCAACTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGt                                                                            >  2:41683/1‑141 (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:200614/1‑141 (MQ=255)
                                                                            ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtt                                                                <  1:384795/141‑1 (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)
                                                                              aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGAACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                              <  1:664744/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)
                                                                                        acacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTACTTTCTTCa                                                    >  1:83238/1‑141 (MQ=255)
                                                                                        acacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCa                                                    >  1:991390/1‑141 (MQ=255)
                                                                                         cacGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCaa                                                   >  2:943791/1‑141 (MQ=255)
                                                                                           cGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAg                                                 >  1:970098/1‑141 (MQ=255)
                                                                                           cGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTAGCCTTTTTCTTTTAATTGCTTCaaaa                                                 >  1:499792/1‑140 (MQ=255)
                                                                                                aCCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTTTTCAAAGATTCa                                            >  2:13921/1‑141 (MQ=255)
                                                                                                aCCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTAATGGTCAATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCa                                            <  2:501440/141‑1 (MQ=255)
                                                                                                 ccTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTTTTTTTTTTTCTTCAATGATTcct                                           >  2:985673/1‑139 (MQ=255)
                                                                                                 ccTTTAAATCAGGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAGATTTCTAGCGTACCATCAGCTGAATTAAAAGTAGCTTTATCATCATGTTCTTTGTGTCCTTTTTCTTTTAGTTTCTTCACAGATGCAt                                           >  2:24471/1‑141 (MQ=255)
                                                                                                   tttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTa                                         >  2:1040480/1‑141 (MQ=255)
                                                                                                    ttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTaa                                        <  1:303081/141‑1 (MQ=255)
                                                                                                    ttAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCAATACCGTTGTATCCTTTTTTTTTTAATTTTTTCAAAGATTCATTaa                                        >  2:433739/1‑141 (MQ=255)
                                                                                                       aaTCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattt                                     <  2:456933/141‑1 (MQ=255)
                                                                                                        aTCATGCAATTCTGTAATTTTTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAatttt                                    <  2:149945/141‑1 (MQ=255)
                                                                                                            tGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTCATTTCTTCAAAGACTCAttctttttcttt                                >  1:321690/1‑141 (MQ=255)
                                                                                                            tGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                                >  2:506039/1‑141 (MQ=255)
                                                                                                            tGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                                >  1:142988/1‑141 (MQ=255)
                                                                                                            tGGAATTCTGTAATTATTGATTTCTGGGGGGTGCCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCATTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                                >  1:534618/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)
                                                                                                                  tCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAGATAGCTTTATCTTCTTTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                          <  1:479517/141‑1 (MQ=255)
                                                                                                                   ttGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCCTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGt                         <  2:185155/140‑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: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)
                                                                                                                                   tttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCa          <  1:887134/141‑1 (MQ=255)
                                                                                                                                      cGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGc       >  1:487486/1‑141 (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)
                                                                                                                                         cATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTATGGTGAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCATGctct    >  1:804565/1‑141 (MQ=255)
                                                                                                                                           tGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGctctct  >  1:923808/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                          
GGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTATTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCTCTCT  >  USA300TCH1516_ALE/2113165‑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: