Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,409 G→A 5.6% intergenic (‑260/‑29) USA300HOU_RS13435 ← / → USA300TCH1516_02485 putative oxidoreductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,608,4090GA 5.6% 631.5 / 16.3 196intergenic (‑260/‑29)USA300HOU_RS13435/USA300TCH1516_02485putative oxidoreductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (81/104);  new base A (5/6);  total (86/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 201 total aligned reads displayed.

TTTTATTTTTAAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGATTT  >  USA300TCH1516_ALE/2608320‑2608499
                                                                                         |                                                                                          
ttttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGa                                                                                           <  2:2955236/91‑1 (MQ=255)
ttttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGa                                                                                           <  2:2953839/91‑1 (MQ=255)
ttttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGa                                                                                           <  2:1801575/91‑1 (MQ=255)
 tttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAt                                                                                          <  1:1005993/91‑1 (MQ=255)
 tttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAt                                                                                          >  1:494135/1‑91 (MQ=255)
  ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                         >  2:1546439/1‑91 (MQ=255)
  ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                         >  2:2793956/1‑91 (MQ=255)
  ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                         >  2:1547064/1‑91 (MQ=255)
   tatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAGAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTg                                                                                        <  2:3415556/91‑1 (MQ=255)
   tatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTg                                                                                        <  1:675444/91‑1 (MQ=255)
   tatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTg                                                                                        <  2:222198/91‑1 (MQ=255)
    atttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGa                                                                                       <  1:41487/91‑1 (MQ=255)
    atttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGa                                                                                       <  1:41551/91‑1 (MQ=255)
     tttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGaa                                                                                      <  1:1990927/91‑1 (MQ=255)
       tttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATg                                                                                    >  2:1717014/1‑91 (MQ=255)
       tttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATg                                                                                    >  2:3831697/1‑91 (MQ=255)
        ttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCAATCGAATGa                                                                                   <  1:177731/91‑1 (MQ=255)
         taAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGAt                                                                                  >  2:1203874/1‑91 (MQ=255)
             aaTGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTaga                                                                              <  2:1661816/91‑1 (MQ=255)
               tGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgagg                                                                            <  1:208626/91‑1 (MQ=255)
               tGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgagg                                                                            >  2:3684322/1‑91 (MQ=255)
                   gAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgagggagg                                                                        <  1:1832172/91‑1 (MQ=255)
                   gAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgagggagg                                                                        >  1:15648/1‑91 (MQ=255)
                    aaaaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGggagga                                                                       <  2:1286536/91‑1 (MQ=255)
                    aaaaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGggagga                                                                       <  2:1404683/91‑1 (MQ=255)
                     aaaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGaa                                                                      >  1:3514239/1‑91 (MQ=255)
                      aaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAAt                                                                     <  1:301050/91‑1 (MQ=255)
                      aaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAAt                                                                     <  2:47020/91‑1 (MQ=255)
                       aaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATa                                                                    <  2:1867539/91‑1 (MQ=255)
                        aTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaa                                                                   >  1:993264/1‑91 (MQ=255)
                        aTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaa                                                                   <  1:713400/91‑1 (MQ=255)
                         tCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaaa                                                                  >  1:1345834/1‑91 (MQ=255)
                            tCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACg                                                               <  2:3160203/91‑1 (MQ=255)
                            tCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACg                                                               <  1:1399113/91‑1 (MQ=255)
                              aCTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTg                                                             <  1:1080826/91‑1 (MQ=255)
                                  cttTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACat                                                         <  1:1120686/91‑1 (MQ=255)
                                  cttTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACat                                                         <  1:904195/91‑1 (MQ=255)
                                        tGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATg                                                   <  2:2727898/91‑1 (MQ=255)
                                         gTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGa                                                  >  2:2034252/1‑91 (MQ=255)
                                           aaaaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGc                                                <  1:3475049/91‑1 (MQ=255)
                                           aaaaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGc                                                >  1:3644419/1‑91 (MQ=255)
                                            aaaaCATTATATTAGTGATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                               <  2:2977818/91‑1 (MQ=255)
                                            aaaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                               >  1:292784/1‑91 (MQ=255)
                                            aaaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                               >  1:292819/1‑91 (MQ=255)
                                             aaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgt                                              <  1:2515675/91‑1 (MQ=255)
                                               aCATTATATTGGGAATAAAGTTAATGCTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                            <  2:2468135/91‑1 (MQ=255)
                                               aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                            <  2:2319160/91‑1 (MQ=255)
                                               aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                            <  2:849456/91‑1 (MQ=255)
                                                 aTTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGttt                                          <  2:3109982/91‑1 (MQ=255)
                                                  ttATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGtttt                                         <  1:2643902/91‑1 (MQ=255)
                                                  ttATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGtttt                                         >  2:397535/1‑91 (MQ=255)
                                                  ttATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGtttt                                         >  2:814487/1‑91 (MQ=255)
                                                   tataTTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTa                                        <  1:1758758/91‑1 (MQ=255)
                                                    atatTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGGGACATATCATGAGCGTGTTTTAg                                       >  2:933093/1‑91 (MQ=255)
                                                     tatTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGc                                      <  1:359388/91‑1 (MQ=255)
                                                      atTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCa                                     <  2:3285036/91‑1 (MQ=255)
                                                      atTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTCACATATCATGAGCGTGTTTTAGCa                                     >  2:3617796/1‑91 (MQ=255)
                                                         aGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTa                                  <  2:1164544/91‑1 (MQ=255)
                                                         aGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTa                                  <  2:612342/91‑1 (MQ=255)
                                                            aataaAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAaga                               <  1:550413/91‑1 (MQ=255)
                                                             ataaAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAagag                              >  2:1487620/1‑91 (MQ=255)
                                                                aaGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAg                           <  2:3865325/91‑1 (MQ=255)
                                                                 aGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGa                          <  2:777676/91‑1 (MQ=255)
                                                                  gTTAATACTGTGATTTATCATTCGATTGAATGATTAGCGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaa                         >  1:3277867/1‑91 (MQ=255)
                                                                  gTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaa                         >  1:2106491/1‑91 (MQ=255)
                                                                  gTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaa                         <  2:2891891/91‑1 (MQ=255)
                                                                   ttAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaaa                        <  2:3025136/91‑1 (MQ=255)
                                                                   ttAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaaa                        >  2:1924491/1‑91 (MQ=255)
                                                                      aTACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTa                     >  2:51313/1‑91 (MQ=255)
                                                                       tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                    <  2:3884222/91‑1 (MQ=255)
                                                                       tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                    >  1:3623471/1‑91 (MQ=255)
                                                                       tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                    <  2:1953920/91‑1 (MQ=255)
                                                                       tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                    <  1:1422165/91‑1 (MQ=255)
                                                                         cTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaaaa                  <  2:3645511/91‑1 (MQ=255)
                                                                           gtgATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGa                >  2:881478/1‑91 (MQ=255)
                                                                            tgATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGaa               >  1:3579211/1‑91 (MQ=255)
                                                                             gATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAAc              <  2:2746596/91‑1 (MQ=255)
                                                                              aTTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAccc             >  2:3546808/1‑91 (MQ=255)
                                                                               tttATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCg            <  2:1310412/91‑1 (MQ=255)
                                                                               tttATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCg            >  1:2309106/1‑91 (MQ=255)
                                                                                ttATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGc           >  1:3843756/1‑91 (MQ=255)
                                                                                 tATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCa          <  2:2844978/91‑1 (MQ=255)
                                                                                 aaTCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCa          <  2:3849074/90‑1 (MQ=255)
                                                                                   tCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAtt        >  1:454002/1‑91 (MQ=255)
                                                                                   tCATTCAATCGAATGATAAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAtt        >  2:1008039/1‑91 (MQ=255)
                                                                                    cATTCGATTGAATGATTAGAGGGAGGAATTAGACGTGACATATCTTGAGCGTGGTTTGGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt       <  1:3102636/91‑1 (MQ=255)
                                                                                    cATTCGATTGAATGATTAGAGGGAGGAATAATACGTGACATATCGTGAGCGTGGTTTGGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt       <  1:3100960/91‑1 (MQ=255)
                                                                                    cATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt       <  1:2833321/91‑1 (MQ=255)
                                                                                    cATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt       <  2:345066/91‑1 (MQ=255)
                                                                                     aTTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTg      >  1:3736829/1‑91 (MQ=255)
                                                                                     aTTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTg      <  2:2003007/91‑1 (MQ=255)
                                                                                     aTTCAATCGACTGATTAGAGGGAGGACTAAAACGTGACATATCCTGAGCGTGGTTTAGCATTAAGCGCAGAATGTAAAAGAACCGCATTTg      >  2:3634632/1‑91 (MQ=255)
                                                                                     aTTCAATCGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTg      >  2:3563264/1‑91 (MQ=255)
                                                                                      ttCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGa     <  1:3394171/91‑1 (MQ=255)
                                                                                       tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt    <  1:1163236/91‑1 (MQ=255)
                                                                                       tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt    >  2:1092386/1‑91 (MQ=255)
                                                                                       tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt    >  1:606716/1‑91 (MQ=255)
                                                                                       tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt    >  2:2548889/1‑91 (MQ=255)
                                                                                        cGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAtt   >  1:3764006/1‑91 (MQ=255)
                                                                                         gATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAttt  >  1:2693193/1‑91 (MQ=255)
                                                                                         |                                                                                          
TTTTATTTTTAAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGATTT  >  USA300TCH1516_ALE/2608320‑2608499

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