Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,779,810 A→T 8.0% intergenic (‑116/+26) USA300TCH1516_01664 ← / ← rplT hypothetical protein/50S ribosomal protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,779,8100AT 8.0% 859.1 / 44.7 266intergenic (‑116/+26)USA300TCH1516_01664/rplThypothetical protein/50S ribosomal protein L20
Reads supporting (aligned to +/- strand):  ref base A (106/138);  new base T (11/11);  total (117/149)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

TTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCGA  >  USA300TCH1516_ALE/1779720‑1779900
                                                                                          |                                                                                          
  gttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATg                                                                                          >  2:778578/1‑91 (MQ=255)
    gctgttAAATGCATTTTAGCAAATATGTCATATTGTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                        <  2:360422/87‑1 (MQ=255)
     ctctTAAATGCATTTTAGCAAGTATGTCATATTTTAGCGATTAAAGGTGGCATGTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                       <  2:1740608/91‑1 (MQ=255)
     ctctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                       <  1:1191278/91‑1 (MQ=255)
     ctctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                       >  1:2640398/1‑91 (MQ=255)
     ctctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                       <  1:3018159/91‑1 (MQ=255)
      tctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAt                                                                                      >  1:670071/1‑91 (MQ=255)
       ctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTAAGAAAAATAAAAAATAAAAAAGAGATGACAATGCCAtc                                                                                     <  2:3955895/91‑1 (MQ=255)
        ttAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtct                                                                                    >  1:1942701/1‑91 (MQ=255)
        ttAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtct                                                                                    >  1:2544741/1‑91 (MQ=255)
         tAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCAtctc                                                                                   <  2:2685128/91‑1 (MQ=39)
            aggCATTTTAGCAAATAGGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCttt                                                                                <  2:2254466/89‑1 (MQ=37)
            aTGCATTTTAGCAATTATGTCTTATTTTAGCGATTGATTGTGGCATTTTGCTATTTAAAAAATAAAATAGAGATCGCATTGTCATCTCttt                                                                                <  1:2862671/91‑1 (MQ=25)
             tGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtttt                                                                               <  1:2167543/91‑1 (MQ=38)
               cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                             >  2:1250174/1‑91 (MQ=255)
               cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                             <  2:2530462/91‑1 (MQ=255)
                aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTta                                                                            >  2:3471727/1‑91 (MQ=255)
                   ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                         >  1:2040246/1‑91 (MQ=255)
                     aGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgt                                                                       <  1:723500/91‑1 (MQ=255)
                      gCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGAAAATGCCATCTCTTTTTTATAAgtg                                                                      >  2:3676820/1‑91 (MQ=255)
                       cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                     <  2:313595/91‑1 (MQ=255)
                       cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                     >  1:435678/1‑91 (MQ=255)
                       cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                     <  1:2958777/91‑1 (MQ=255)
                        aaaTATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtt                                                                    <  1:1695283/91‑1 (MQ=255)
                          atatGTCATATTTTAGCGATTAAATGTGTCATTATACAATTAAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttat                                                                  <  2:3900738/91‑1 (MQ=38)
                           tatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatt                                                                 <  1:1597147/91‑1 (MQ=255)
                            atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                                >  2:26678/1‑91 (MQ=255)
                              gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTattatt                                                              >  2:1051192/1‑91 (MQ=255)
                                cATATTTTAGCGATTACATGTGGCATTTTGCTAGTTAAAAAATAAAAAGGAGATGACAATGCCATCCCTTGCTTATAAGTGTTGTTAtttt                                                            >  2:6552/1‑91 (MQ=255)
                                     tttAGCGAGTAAATGTGGGATGTTGCTATTTAAAAAATGAAAAAGAGGTGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGc                                                       <  2:3391419/91‑1 (MQ=255)
                                      ttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCa                                                      >  2:2479003/1‑91 (MQ=255)
                                       tAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCAt                                                     <  1:363415/91‑1 (MQ=255)
                                        aGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                    <  1:1862083/91‑1 (MQ=255)
                                         gCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATct                                                   >  1:1069746/1‑91 (MQ=255)
                                          cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                                  <  2:3815480/91‑1 (MQ=255)
                                          cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                                  >  1:3878530/1‑91 (MQ=255)
                                          cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                                  <  2:2117787/91‑1 (MQ=255)
                                            aTTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACCATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATcttta                                                >  2:2947589/1‑91 (MQ=255)
                                            aTTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATcttta                                                <  1:2079589/91‑1 (MQ=255)
                                                   gtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctttag                                         >  1:3208202/1‑91 (MQ=255)
                                                    tgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGt                                        <  2:682377/91‑1 (MQ=255)
                                                     ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                       >  1:2501637/1‑91 (MQ=255)
                                                     ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                       >  1:2330484/1‑91 (MQ=255)
                                                     ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                       >  2:3218884/1‑91 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                      <  1:2366928/91‑1 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                      <  1:70484/91‑1 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                      <  2:2910546/91‑1 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                      <  1:1738918/91‑1 (MQ=255)
                                                       cATTTTGCTGTTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                     <  2:681999/91‑1 (MQ=255)
                                                       cATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                     <  2:2181643/91‑1 (MQ=255)
                                                       cATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                     <  1:74744/91‑1 (MQ=255)
                                                       cATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                     <  2:3771156/91‑1 (MQ=255)
                                                        aTTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACt                                    >  2:382132/1‑91 (MQ=255)
                                                         ttttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTa                                   >  2:1956893/1‑91 (MQ=255)
                                                          tttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTaa                                  >  1:851794/1‑91 (MQ=255)
                                                           ttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAAt                                 <  2:3464288/91‑1 (MQ=255)
                                                             gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                               >  1:1536418/1‑91 (MQ=255)
                                                             gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                               >  2:71957/1‑91 (MQ=255)
                                                             gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                               >  1:571486/1‑91 (MQ=255)
                                                              cTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGa                              <  1:244571/91‑1 (MQ=255)
                                                              cTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGa                              <  1:900459/91‑1 (MQ=255)
                                                              cTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGa                              >  1:504684/1‑91 (MQ=255)
                                                               tATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAg                             <  1:919597/91‑1 (MQ=255)
                                                               tATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAg                             >  1:874685/1‑91 (MQ=255)
                                                                 tttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAATTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCa                           >  2:2502639/1‑91 (MQ=255)
                                                                 tttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCa                           <  2:1118927/91‑1 (MQ=255)
                                                                  ttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaa                          <  1:2731211/91‑1 (MQ=255)
                                                                  ttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaa                          >  2:2792112/1‑91 (MQ=255)
                                                                   taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                         <  1:597896/91‑1 (MQ=255)
                                                                   taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                         <  2:752578/91‑1 (MQ=255)
                                                                   taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                         <  2:3263493/91‑1 (MQ=255)
                                                                     aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                       <  2:1814429/91‑1 (MQ=255)
                                                                      aaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGc                      >  2:891329/1‑91 (MQ=255)
                                                                       aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTATTTGAGCAAATGCt                     >  2:2457331/1‑91 (MQ=255)
                                                                       aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCt                     <  2:1638855/91‑1 (MQ=255)
                                                                            aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTATTTGAGCAAATGCTTTTTc                >  2:3533593/1‑91 (MQ=255)
                                                                             aaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg               <  2:1582376/91‑1 (MQ=255)
                                                                             aaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg               <  2:1581938/91‑1 (MQ=255)
                                                                             aaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg               <  1:3579374/91‑1 (MQ=255)
                                                                             aaaaGAGATGACAATGCCATCTCCTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg               >  1:3133472/1‑91 (MQ=255)
                                                                              aaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt              >  2:2566618/1‑91 (MQ=38)
                                                                              aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt              <  2:200768/91‑1 (MQ=255)
                                                                              aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt              <  1:2285037/91‑1 (MQ=255)
                                                                               aaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc             >  1:2804244/1‑91 (MQ=38)
                                                                               aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc             >  2:3058818/1‑91 (MQ=255)
                                                                               aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTCTAGTTACTAATTGAGCACATGCTTTTTCGTc             >  2:582503/1‑91 (MQ=255)
                                                                                agagATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa            <  2:1380964/91‑1 (MQ=255)
                                                                                 gagaTGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAg           >  2:2037786/1‑91 (MQ=38)
                                                                                    aTGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGaaa        <  1:963569/91‑1 (MQ=255)
                                                                                     tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt       >  1:1949789/1‑91 (MQ=255)
                                                                                     tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt       <  2:1978058/91‑1 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTGCTAATTGAGCAAATGCTTTTTCGTCAGAAAtt      >  1:3377418/1‑91 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt      <  1:1756520/91‑1 (MQ=255)
                                                                                       aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg     >  1:3764937/1‑91 (MQ=255)
                                                                                        cAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATCGc    <  1:1385747/91‑1 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgca   <  1:3834260/91‑5 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgca   <  1:413915/91‑5 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgca   <  1:806268/91‑5 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCg   >  1:1999618/1‑91 (MQ=255)
                                                                                          aTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCGa  >  1:330966/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCGA  >  USA300TCH1516_ALE/1779720‑1779900

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