Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,412 T→C 5.7% intergenic (‑263/‑26) 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,4120TC 5.7% 580.9 / 16.2 186intergenic (‑263/‑26)USA300HOU_RS13435/USA300TCH1516_02485putative oxidoreductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (75/100);  new base C (5/6);  total (80/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

TTATTTTTAAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGATTTTCG  >  USA300TCH1516_ALE/2608322‑2608502
                                                                                          |                                                                                          
ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                            >  2:1546872/1‑91 (MQ=255)
ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                            <  2:3759351/91‑1 (MQ=255)
ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                            >  1:3741403/1‑91 (MQ=255)
ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGAtt                                                                                            >  2:2793956/1‑91 (MQ=255)
ttatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATATATCATTCGAtt                                                                                            >  2:3232458/1‑91 (MQ=255)
 tatttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAGAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTg                                                                                           <  2:3415556/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)
  atttttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGa                                                                                          <  1:2970714/91‑1 (MQ=255)
     tttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATg                                                                                       >  2:3831697/1‑91 (MQ=255)
     tttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCAATCGAATg                                                                                       <  1:555661/91‑1 (MQ=255)
      ttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCAATCGAATGa                                                                                      <  1:3645985/91‑1 (MQ=255)
      ttaAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCAATCGAATGa                                                                                      <  1:177731/91‑1 (MQ=255)
       taAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGAt                                                                                     >  2:1203874/1‑91 (MQ=255)
        aaaaaTGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGAtt                                                                                    <  1:3306093/91‑1 (MQ=255)
        aaaaaTGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCAATCGAATGAta                                                                                    <  1:2054421/91‑2 (MQ=255)
         aaaaTGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTa                                                                                   <  2:3653438/91‑1 (MQ=255)
              gTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgaggg                                                                              <  2:1548677/91‑1 (MQ=255)
              gTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgaggg                                                                              >  2:1797586/1‑91 (MQ=255)
                 gAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAgagggagg                                                                           >  1:2600976/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)
                   aaaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGaa                                                                         >  1:3514927/1‑91 (MQ=255)
                    aaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAAt                                                                        <  2:47020/91‑1 (MQ=255)
                    aaaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAACTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAAt                                                                        <  1:3096876/91‑1 (MQ=255)
                     aaTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATa                                                                       <  2:1867539/91‑1 (MQ=255)
                      aTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaa                                                                      >  1:993264/1‑91 (MQ=255)
                      aTCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaa                                                                      <  1:713383/91‑1 (MQ=255)
                       tCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATaaa                                                                     >  1:1345834/1‑91 (MQ=255)
                          tCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACg                                                                  <  1:1399113/91‑1 (MQ=255)
                          tCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACg                                                                  >  1:3769493/1‑91 (MQ=255)
                            aCTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTg                                                                <  1:1080826/91‑1 (MQ=255)
                            aCTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTg                                                                >  1:732601/1‑91 (MQ=255)
                            aCTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTg                                                                <  1:1080806/91‑1 (MQ=255)
                                cttTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACat                                                            <  1:1120686/91‑1 (MQ=255)
                                cttTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACat                                                            <  1:904195/91‑1 (MQ=255)
                                cttTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACat                                                            <  2:1383287/91‑1 (MQ=255)
                                      tGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATg                                                      <  2:2727898/91‑1 (MQ=255)
                                       gTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGa                                                     >  2:2034252/1‑91 (MQ=255)
                                          aaaaCATTATATTAGTTATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                                  <  2:2975244/91‑1 (MQ=255)
                                          aaaaCATTATATTAGTGATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                                  <  2:2978163/91‑1 (MQ=255)
                                          aaaaCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCg                                                  >  1:292807/1‑91 (MQ=255)
                                             aCATTATATTGGGAATAAAGTTAATGCTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                               <  2:2468135/91‑1 (MQ=255)
                                             aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGGCATATCATGAGCgtgt                                               >  1:3511436/1‑91 (MQ=255)
                                             aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                               >  1:2218881/1‑91 (MQ=255)
                                             aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                               <  2:2319160/91‑1 (MQ=255)
                                             aCATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCgtgt                                               <  2:849456/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)
                                                 tataTTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTa                                           >  2:1492192/1‑91 (MQ=255)
                                                  atatTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGGGACATATCATGAGCGTGTTTTAg                                          >  2:933093/1‑91 (MQ=255)
                                                   tatTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGc                                         <  1:359312/91‑1 (MQ=255)
                                                   tatTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGc                                         >  2:1713973/1‑91 (MQ=255)
                                                    atTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTCACATATCATGAGCGTGTTTTAGCa                                        >  2:3617796/1‑91 (MQ=255)
                                                     ttAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCAt                                       >  1:1960113/1‑91 (MQ=255)
                                                       aGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTa                                     <  2:612342/91‑1 (MQ=255)
                                                       aGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTa                                     <  2:1164544/91‑1 (MQ=255)
                                                              aaGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAg                              <  2:3865325/91‑1 (MQ=255)
                                                               aGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGa                             >  1:527253/1‑91 (MQ=255)
                                                               aGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGa                             <  2:777676/91‑1 (MQ=255)
                                                                gTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaa                            >  1:2106397/1‑91 (MQ=255)
                                                                gTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaa                            <  2:2891891/91‑1 (MQ=255)
                                                                 ttAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGaaa                           >  2:1924491/1‑91 (MQ=255)
                                                                 ttAATACTGTGATTTATCATTCGATTGAATGATTAGACGGAGGAATAAAACCTGACATATCATGATCGTGTTTTAGCATTAAGAGTAGaaa                           >  1:3607908/1‑91 (MQ=255)
                                                                  tAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAg                          >  1:2413904/1‑91 (MQ=255)
                                                                  tAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAg                          <  2:1411028/91‑1 (MQ=255)
                                                                    aTACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTa                        >  2:51313/1‑91 (MQ=255)
                                                                    aTACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTa                        <  2:3068977/91‑1 (MQ=255)
                                                                    aTACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTa                        >  1:634220/1‑91 (MQ=255)
                                                                     tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                       <  2:1953920/91‑1 (MQ=255)
                                                                     tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                       <  1:1422165/91‑1 (MQ=255)
                                                                     tACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTaa                       <  2:3884222/91‑1 (MQ=255)
                                                                         gtgATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGa                   >  2:881478/1‑91 (MQ=255)
                                                                         gtgATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGa                   >  1:2662059/1‑91 (MQ=255)
                                                                          tgATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGaa                  >  2:990650/1‑91 (MQ=255)
                                                                           gATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAAc                 <  2:2746596/91‑1 (MQ=255)
                                                                            aTTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAccc                >  2:3546808/1‑91 (MQ=255)
                                                                             tttATCATTCGATTGAATGATTAGAGGGAGGAATAAAGCGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCg               >  1:3265118/1‑91 (MQ=255)
                                                                             tttATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCg               <  2:1310412/91‑1 (MQ=255)
                                                                              ttATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGc              >  1:3843756/1‑91 (MQ=255)
                                                                               tATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCa             >  1:2809957/1‑91 (MQ=255)
                                                                               tATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCa             <  2:2844978/91‑1 (MQ=255)
                                                                               aaTCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCa             <  2:3849074/90‑1 (MQ=255)
                                                                                 tCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACTGCAtt           >  2:3118908/1‑91 (MQ=255)
                                                                                 tCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAtt           >  1:3643472/1‑91 (MQ=255)
                                                                                 tCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAtt           >  1:454002/1‑91 (MQ=255)
                                                                                 tCATTCAATCGAATGATAAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAtt           >  2:1008039/1‑91 (MQ=255)
                                                                                  cATTCGATTGAATGATTAGAGGGAGGAATAAGACGTGACATATCGTGAGCGTGCTTTGGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt          <  1:3101064/91‑1 (MQ=255)
                                                                                  cATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCAttt          <  2:345066/91‑1 (MQ=255)
                                                                                   aTTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTg         <  2:2003007/91‑1 (MQ=255)
                                                                                   aTTCAATCGACTGATTAGAGGGAGGACTAAAACGTGACATATCCTGAGCGTGGTTTAGCATTAAGCGCAGAATGTAAAAGAACCGCATTTg         >  2:3634632/1‑91 (MQ=255)
                                                                                   aTTCAATCGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTg         >  2:3563264/1‑91 (MQ=255)
                                                                                     tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt       >  2:2548889/1‑91 (MQ=255)
                                                                                     tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt       >  1:2913355/1‑91 (MQ=255)
                                                                                     tCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAt       <  1:1163236/91‑1 (MQ=255)
                                                                                       gATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGAttt     >  1:3447145/1‑91 (MQ=255)
                                                                                         ttGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGATTttc   >  2:3371865/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTATTTTTAAAAATGTTGAAAAATCGTCACTTCTTTATTGTAAAAACATTATATTAGTAATAAAGTTAATACTGTGATTTATCATTCGATTGAATGATTAGAGGGAGGAATAAAACGTGACATATCATGAGCGTGTTTTAGCATTAAGAGCAGAAAGTAAAAGAACCGCATTTGATTTTCG  >  USA300TCH1516_ALE/2608322‑2608502

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: