Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 868,373 A→T 5.5% intergenic (+56/+1266) smpB → / ← USA300HOU_RS04240 SsrA‑binding protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE868,3730AT 5.5% 673.0 / 20.7 199intergenic (+56/+1266)smpB/USA300HOU_RS04240SsrA‑binding protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (104/84);  new base T (8/3);  total (112/87)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 211 total aligned reads displayed.

GAGATGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAGGG  >  USA300TCH1516_ALE/868283‑868463
                                                                                          |                                                                                          
 agaTGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAAt                                                                                           >  2:811790/1‑91 (MQ=255)
  gaTGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATATGCTTTTGTAAAtt                                                                                          <  2:3586096/91‑1 (MQ=255)
   aTGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATGTGTGCTTTTGTAAATTa                                                                                         >  2:544950/1‑91 (MQ=255)
   aTGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTa                                                                                         >  1:3534124/1‑91 (MQ=255)
   aTGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTa                                                                                         >  2:739535/1‑91 (MQ=255)
       tGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACaaa                                                                                     >  2:2682621/1‑91 (MQ=255)
        gcgcgcGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACaaaa                                                                                    >  1:680144/1‑91 (MQ=255)
         cgcgcgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAg                                                                                   >  2:52548/1‑91 (MQ=255)
          gcgcgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGt                                                                                  >  2:3198038/1‑91 (MQ=255)
          gcgcgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGt                                                                                  >  2:2380754/1‑91 (MQ=255)
           cgcgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTa                                                                                 <  1:3436864/91‑1 (MQ=255)
            gcgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTAt                                                                                >  2:3124360/1‑91 (MQ=255)
             cgATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATg                                                                               >  1:1241220/1‑91 (MQ=255)
               atatGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAt                                                                             >  1:41563/1‑91 (MQ=255)
               atatGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAt                                                                             <  2:342718/91‑1 (MQ=255)
               atatGAAAGCCCGTTATTAAGCGATTTAGTTGCCTAATCGGGCTATCTTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAt                                                                             >  2:718995/1‑91 (MQ=255)
                tatGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtt                                                                            <  1:1342210/91‑1 (MQ=255)
                tatGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtt                                                                            <  1:1677012/91‑1 (MQ=255)
                 atGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGCttt                                                                           >  1:3011586/1‑91 (MQ=255)
                 atGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAttt                                                                           >  2:700376/1‑91 (MQ=255)
                 atGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAttt                                                                           >  2:1582414/1‑91 (MQ=255)
                 atGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAttt                                                                           >  1:818550/1‑91 (MQ=255)
                  tGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtttg                                                                          >  2:1317793/1‑91 (MQ=255)
                   gAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtttgt                                                                         >  1:2746348/1‑91 (MQ=255)
                   gAAAGCCCGTTATTAAGCGATTTAGTTGATTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtttgt                                                                         >  1:3288980/1‑91 (MQ=255)
                     aaGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtttgttt                                                                       <  1:3540038/91‑1 (MQ=255)
                     aaGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGAtttgttt                                                                       >  1:1886983/1‑91 (MQ=255)
                         ccGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGAtt                                                                   <  2:1821010/91‑1 (MQ=255)
                          cGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGAttt                                                                  <  2:1142243/91‑1 (MQ=255)
                          cGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGAttt                                                                  <  1:1780818/91‑1 (MQ=255)
                            ttattaAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGCGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATttat                                                                <  2:489689/91‑1 (MQ=255)
                             tattaAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATttatt                                                               >  1:3688032/1‑91 (MQ=255)
                              attaAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATttatta                                                              <  1:1492226/91‑1 (MQ=255)
                              attaAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATttatta                                                              >  2:12319/1‑91 (MQ=255)
                              attaAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATttatta                                                              >  2:3741954/1‑91 (MQ=255)
                                 aaGCGATTTAGTTGCTTAATCGGGCTATATTTGATATCGTTATATGTGCTTCTGTAAATTACAAATGTATGATTTGTTTGATTTACTAttt                                                           >  1:3605823/1‑91 (MQ=255)
                                 aaGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTAttt                                                           <  2:1398219/91‑1 (MQ=255)
                                    cGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCgg                                                        <  1:1715748/91‑1 (MQ=255)
                                    cGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCgg                                                        >  2:3820414/1‑91 (MQ=255)
                                     gATTTAGTTGCTTAATCGGGCTATATTTGAGATAGTTATATGTGCTTTTGTAAAGTACAAAAGTATGATTTGTTTGATTTATTATTTCggg                                                       <  2:3861594/91‑1 (MQ=255)
                                      aTTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCgggg                                                      <  1:3295478/91‑1 (MQ=255)
                                      aTTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCgggg                                                      >  1:2424950/1‑91 (MQ=255)
                                          aGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGt                                                  <  1:3819884/91‑1 (MQ=255)
                                           gTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGtt                                                 >  2:3444296/1‑91 (MQ=255)
                                           gTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGtt                                                 >  1:760054/1‑91 (MQ=255)
                                              gCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCAt                                              <  1:3372638/91‑1 (MQ=255)
                                              gCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCAt                                              >  2:849029/1‑91 (MQ=255)
                                               cTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATg                                             <  1:1010750/91‑1 (MQ=255)
                                               cTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATg                                             <  2:1700036/91‑1 (MQ=255)
                                               cTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATg                                             >  1:1286767/1‑91 (MQ=255)
                                               cTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATg                                             <  2:955709/91‑1 (MQ=255)
                                                 tAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGa                                           <  1:438769/91‑1 (MQ=255)
                                                  aaTCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGAt                                          <  2:2614887/91‑1 (MQ=255)
                                                    tCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTc                                        >  1:2082994/1‑91 (MQ=255)
                                                     cGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCg                                       >  2:1574431/1‑91 (MQ=255)
                                                       ggCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGAc                                     <  2:1915033/91‑1 (MQ=255)
                                                        gCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACa                                    >  2:314720/1‑91 (MQ=255)
                                                        gCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACa                                    >  1:2144468/1‑91 (MQ=255)
                                                         cTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAg                                   <  2:3651840/91‑1 (MQ=255)
                                                         cTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAg                                   <  2:1953331/91‑1 (MQ=255)
                                                         cTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAg                                   >  2:652332/1‑91 (MQ=255)
                                                           atatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAggg                                 <  1:1608525/91‑1 (MQ=255)
                                                            tatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAgggg                                >  2:293178/1‑91 (MQ=255)
                                                            tatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAgggg                                <  2:352174/91‑1 (MQ=255)
                                                            tatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAgggg                                <  1:3764282/91‑1 (MQ=255)
                                                            tatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAgggg                                <  1:3764457/91‑1 (MQ=255)
                                                            tatTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAgggg                                <  1:495700/91‑1 (MQ=255)
                                                              tttGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTc                              <  1:90121/91‑1 (MQ=255)
                                                                  atatAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTccccc                          <  2:1551778/91‑1 (MQ=255)
                                                                  atatAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTccccc                          <  1:3721343/91‑1 (MQ=255)
                                                                    ataGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGa                        <  1:267788/91‑1 (MQ=255)
                                                                    ataGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGa                        >  2:1594290/1‑91 (MQ=255)
                                                                    ataGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGa                        >  2:1594197/1‑91 (MQ=255)
                                                                     taGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAg                       <  2:1723316/91‑1 (MQ=255)
                                                                     taGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAg                       >  1:2920299/1‑91 (MQ=255)
                                                                      aGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGc                      <  1:1163753/91‑1 (MQ=255)
                                                                            atGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTa                >  1:899907/1‑91 (MQ=255)
                                                                              atacTTTTGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAg              >  1:555211/4‑91 (MQ=255)
                                                                               tgCTTTTGTAAATTACAACAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGc             >  2:3084717/1‑91 (MQ=255)
                                                                               tacTTTTGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGc             >  2:2552648/3‑91 (MQ=255)
                                                                               tacTTTTGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGc             >  1:1745321/3‑91 (MQ=255)
                                                                               tacTTTTGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGc             >  2:615106/3‑91 (MQ=255)
                                                                                gCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCg            >  2:3282097/1‑91 (MQ=255)
                                                                                gCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCg            <  2:432380/91‑1 (MQ=255)
                                                                                gCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCg            <  2:636231/91‑1 (MQ=255)
                                                                                 cTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgt           <  1:194167/91‑1 (MQ=255)
                                                                                 cTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgt           <  1:2487577/91‑1 (MQ=255)
                                                                                  ttttGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgtg          >  2:2174258/1‑91 (MQ=255)
                                                                                   tttGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgtgt         >  1:2605587/1‑91 (MQ=255)
                                                                                   tttGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgtgt         <  2:2255440/91‑1 (MQ=255)
                                                                                   tttGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCgtgt         <  2:1106723/91‑1 (MQ=255)
                                                                                    ttGTAATTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTc        >  2:2674832/1‑91 (MQ=255)
                                                                                      gTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCgg      >  1:3165559/1‑91 (MQ=255)
                                                                                       tAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGa     >  1:2540321/1‑91 (MQ=255)
                                                                                       tAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGa     >  1:2301961/1‑91 (MQ=255)
                                                                                       tAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGa     >  1:1854895/1‑91 (MQ=255)
                                                                                         aaTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAgg   >  1:2993934/1‑91 (MQ=255)
                                                                                         aaTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAgg   >  1:2994338/1‑91 (MQ=255)
                                                                                         aaTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAgg   <  2:1130048/91‑1 (MQ=255)
                                                                                         aaTTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAgg   >  2:1327761/1‑91 (MQ=255)
                                                                                          |                                                                                          
GAGATGTTGCGCGCGATATGAAAGCCCGTTATTAAGCGATTTAGTTGCTTAATCGGGCTATATTTGATATAGTTATATGTGCTTTTGTAAATTACAAAAGTATGATTTGTTTGATTTATTATTTCGGGGACGTTCATGGATTCGACAGGGGTCCCCCGAGCTCATTAAGCGTGTCGGAGGG  >  USA300TCH1516_ALE/868283‑868463

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