Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,504 A→T 8.5% intergenic (+82/+27) USA300HOU_RS11475 → / ← pyrG hypothetical protein/CTP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,243,5040AT 8.5% 774.5 / 47.5 235intergenic (+82/+27)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base A (112/103);  new base T (16/4);  total (128/107)
Fisher's exact test for biased strand distribution p-value = 1.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 245 total aligned reads displayed.

GAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGA  >  USA300TCH1516_ALE/2243414‑2243594
                                                                                          |                                                                                          
gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                            <  1:1027742/91‑1 (MQ=255)
gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                            >  2:1352064/1‑91 (MQ=255)
 aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                           >  1:2787863/1‑91 (MQ=255)
  aTTATAAAAGTCACTATTGTTATTTCATACCCTGTCACTAACATAACATTAACATTATACGCAAAACCGATACTTTCAATGCCTAACAAtc                                                                                          >  2:463948/1‑90 (MQ=255)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  1:12082/91‑1 (MQ=38)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  1:794075/91‑1 (MQ=38)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  1:507828/91‑1 (MQ=38)
    tataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaa                                                                                        >  1:1768636/1‑91 (MQ=255)
     ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                       <  1:385992/91‑1 (MQ=255)
     ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                       >  1:2367030/1‑91 (MQ=255)
        aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAACTTATACGAAAAACCGGTACTTTCAATGTTGAACAATGAAAGTa                                                                                    >  1:3878224/1‑91 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   <  1:638821/91‑1 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   >  2:3096707/1‑91 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   <  1:2902706/91‑1 (MQ=255)
          aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                                  <  1:2261808/91‑1 (MQ=255)
           gTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                                 <  1:2735119/91‑1 (MQ=255)
            tCACTATTGTTATTTCATACCTTGTCACTAACGTAATATTAAAATTATACGAAAAGCCGGTACGTTCAATGTTTAACAATGAAAGTACCgg                                                                                <  1:611693/91‑1 (MQ=255)
              aCTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtt                                                                              <  2:3784937/91‑1 (MQ=255)
                tATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtttt                                                                            <  2:2252058/91‑1 (MQ=255)
                 aTTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAATCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTtta                                                                           <  2:3594338/91‑1 (MQ=255)
                    gTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatta                                                                        <  1:1304691/91‑1 (MQ=255)
                     ttATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAg                                                                       <  2:2006044/91‑1 (MQ=255)
                     ttATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAg                                                                       <  1:1610172/91‑1 (MQ=255)
                     ttATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAg                                                                       >  1:223274/1‑91 (MQ=255)
                      tATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGc                                                                      <  2:1660045/91‑1 (MQ=255)
                         ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTCCCGGTTTTATTAGCaaa                                                                   <  1:2532845/91‑1 (MQ=255)
                         ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                                   >  1:1895058/1‑91 (MQ=255)
                          tCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAat                                                                  <  2:2057665/91‑1 (MQ=255)
                           cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                                 <  1:2377888/91‑1 (MQ=255)
                             tACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAattta                                                               <  1:1690869/91‑1 (MQ=255)
                              aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                              <  2:2684491/91‑1 (MQ=255)
                               ccTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                             <  1:1387370/91‑1 (MQ=255)
                               ccTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                             <  1:644933/91‑1 (MQ=255)
                                cTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattt                                                            >  2:3176990/1‑91 (MQ=255)
                                 ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                           >  1:899304/1‑91 (MQ=255)
                                  tGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                          <  2:3821199/91‑1 (MQ=255)
                                  tGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                          <  2:140337/91‑1 (MQ=255)
                                     cACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTAtttt                                                       <  1:1515833/91‑1 (MQ=255)
                                      aCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttg                                                      <  1:129135/91‑1 (MQ=255)
                                       cTAACATAATATTAAAATTCTACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                                     >  2:474843/1‑91 (MQ=255)
                                       cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                                     >  2:668892/1‑91 (MQ=255)
                                        tAACATAATATTAAAATTATACTAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAGTTTATTTATTttgtt                                                    <  1:2104087/91‑1 (MQ=255)
                                         aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                                   <  1:23151/91‑1 (MQ=255)
                                          aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                                  >  1:2000152/1‑91 (MQ=255)
                                              aataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtt                                              <  1:1908640/91‑1 (MQ=255)
                                              aataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtt                                              <  1:463841/91‑1 (MQ=255)
                                               atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                             >  2:267503/1‑91 (MQ=255)
                                               atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                             >  1:3655814/1‑91 (MQ=255)
                                               atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTACTTATTTTGTTGATAttt                                             >  2:1796292/1‑91 (MQ=255)
                                                tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                            >  2:60143/1‑91 (MQ=255)
                                                tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                            >  1:796028/1‑91 (MQ=255)
                                                 atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                           >  2:811701/1‑91 (MQ=255)
                                                 atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                           <  2:2274598/91‑1 (MQ=255)
                                                 atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                           >  1:277623/1‑91 (MQ=255)
                                                   tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                                         >  2:654487/1‑91 (MQ=255)
                                                    aaaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                                        >  2:151896/1‑91 (MQ=255)
                                                       aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAg                                     >  2:3220928/1‑91 (MQ=255)
                                                             cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                               >  2:2784125/1‑91 (MQ=255)
                                                             cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                               <  1:983055/91‑1 (MQ=255)
                                                             cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                               >  1:659314/1‑91 (MQ=255)
                                                              gAAAAACCGGTACTTTCTATGTTTAACAATGAGAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                              <  1:3323795/91‑1 (MQ=255)
                                                              gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                              >  1:336939/1‑91 (MQ=255)
                                                               aaaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATa                             >  2:1892727/1‑91 (MQ=255)
                                                                 aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                           <  2:3332155/91‑1 (MQ=255)
                                                                 aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                           >  2:2456888/1‑91 (MQ=255)
                                                                  aaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAAt                          <  1:3724313/91‑1 (MQ=255)
                                                                   aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg                         <  1:320394/91‑1 (MQ=255)
                                                                    ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                        >  1:751189/1‑91 (MQ=255)
                                                                     cGGTACTTTCATTGTTAAACATTGAACGTACCGGTTTTATTAGCTAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAATTGAt                       >  2:5630/1‑91 (MQ=25)
                                                                     cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                       >  2:117799/1‑91 (MQ=255)
                                                                     cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                       <  2:951741/91‑1 (MQ=255)
                                                                     cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                       <  1:2160077/91‑1 (MQ=255)
                                                                     cGGGACTTTCAATGTTCAACAATGAAAGTACCGGTTCTAGTGGCAAATTTATTTATTTTGCTGCTATTTTAATGAAGATTCAATAAATGAt                       >  2:2885455/1‑91 (MQ=255)
                                                                      ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                      <  2:223274/91‑1 (MQ=255)
                                                                      ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTCATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                      <  2:3065692/91‑1 (MQ=255)
                                                                       gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  1:2511167/1‑91 (MQ=255)
                                                                       gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  2:3199040/1‑91 (MQ=255)
                                                                       gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     <  1:2634698/91‑1 (MQ=255)
                                                                        tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa                    >  2:3486890/1‑91 (MQ=255)
                                                                        tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa                    >  1:1799275/1‑91 (MQ=255)
                                                                          cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa                  <  2:2335825/91‑1 (MQ=255)
                                                                          cTTTCAACGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa                  <  1:1166393/91‑1 (MQ=255)
                                                                             tcattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt               >  1:182438/5‑91 (MQ=38)
                                                                             tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt               <  2:1222272/91‑1 (MQ=255)
                                                                             tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt               <  1:1117740/91‑1 (MQ=255)
                                                                              cattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc              >  1:2187226/4‑91 (MQ=38)
                                                                               attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             >  1:2536678/3‑91 (MQ=38)
                                                                               aaTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAttg             >  2:1944721/1‑89 (MQ=255)
                                                                               aaTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             <  1:22970/91‑1 (MQ=255)
                                                                                aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg            >  1:3448301/1‑91 (MQ=255)
                                                                                aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg            >  2:2521992/1‑91 (MQ=255)
                                                                                  gTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGt          >  2:3674305/1‑91 (MQ=255)
                                                                                     tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg       >  1:3709928/1‑91 (MQ=255)
                                                                                     tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg       <  1:224484/91‑1 (MQ=255)
                                                                                      aaCAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCgg      <  2:3644401/91‑1 (MQ=255)
                                                                                       aCAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGa     >  1:3025460/1‑91 (MQ=255)
                                                                                        cAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGAc    >  1:102461/1‑91 (MQ=255)
                                                                                         aaTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACg   <  1:2569549/91‑1 (MQ=255)
                                                                                         aaTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACg   >  1:2609169/1‑91 (MQ=255)
                                                                                          aTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa  <  2:2567729/91‑1 (MQ=255)
                                                                                          |                                                                                          
GAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGA  >  USA300TCH1516_ALE/2243414‑2243594

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: