Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,494 A→T 7.5% intergenic (+72/+37) 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,4940AT 7.5% 785.9 / 37.9 234intergenic (+72/+37)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base A (107/109);  new base T (3/15);  total (110/124)
Fisher's exact test for biased strand distribution p-value = 7.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

AAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGG  >  USA300TCH1516_ALE/2243404‑2243584
                                                                                          |                                                                                          
aaGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCaa                                                                                            >  2:2343242/1‑91 (MQ=255)
  gagaCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATg                                                                                          >  1:1449355/1‑91 (MQ=255)
    gaCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGtt                                                                                        <  1:449944/91‑1 (MQ=255)
      cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTAGTTTCAATGTTTa                                                                                      <  1:2972620/91‑1 (MQ=255)
      cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTAGTTTCAATGTTTa                                                                                      <  1:2971819/91‑1 (MQ=255)
      cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGttaa                                                                                      <  1:2323497/91‑3 (MQ=255)
      cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                      <  1:1238577/91‑1 (MQ=255)
       ttAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTaa                                                                                     >  2:129135/1‑91 (MQ=255)
         aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACa                                                                                   >  2:1380064/1‑91 (MQ=255)
         aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACa                                                                                   >  1:2202884/1‑91 (MQ=255)
          gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                  <  2:63490/91‑2 (MQ=255)
          gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                  <  1:1027359/91‑1 (MQ=255)
          gAATTATAAAAGTCACTATTGTTATTTCATAACTTGTCACTAATATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                  <  2:2680816/91‑2 (MQ=38)
           aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACatt                                                                                 <  2:3243759/91‑3 (MQ=38)
           aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                 >  1:1644570/1‑91 (MQ=255)
           aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                 >  1:2787863/1‑91 (MQ=255)
            aTTATAAAGGGCACTATTGTTATTTCATACCTTGTCACTAACATAATATTTAAATTATACGAAAAACCGGTACTTTCATTGTTAAACattg                                                                                <  2:3594438/91‑4 (MQ=37)
            aTTATAAAAGTCACTATTGTTATTTCATACCCTGTCACTAACATAACATTAACATTATACGCAAAACCGATACTTTCAATGCCTAACAAtc                                                                                >  2:463948/1‑90 (MQ=255)
             ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                               <  1:794075/91‑1 (MQ=38)
             ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                               <  1:12082/91‑1 (MQ=38)
               ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                             <  1:385945/91‑1 (MQ=255)
                 aaaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGt                                                                           >  1:2079566/1‑91 (MQ=255)
                   aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                         <  1:638821/91‑1 (MQ=255)
                   aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                         <  2:3580947/91‑1 (MQ=255)
                    aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                        <  2:1986963/91‑1 (MQ=255)
                    aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                        <  1:2835994/91‑1 (MQ=255)
                     gTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                       <  1:2735119/91‑1 (MQ=255)
                      tCACTATTGTTATTTCATACCTTGTCACTAACGTAATATTAAAATTATACGAAAAGCCGGTACGTTCAATGTTTAACAATGAAAGTACCgg                                                                      <  1:611693/91‑1 (MQ=255)
                            ttGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttat                                                                >  1:783074/1‑91 (MQ=255)
                              gTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatta                                                              <  1:1304690/91‑1 (MQ=255)
                                tATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGc                                                            <  2:1660045/91‑1 (MQ=255)
                                 aTTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCa                                                           >  1:1956392/1‑91 (MQ=255)
                                   ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTCCCGGTTTTATTAGCaaa                                                         <  1:2532845/91‑1 (MQ=255)
                                   ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                         >  1:1277638/1‑91 (MQ=255)
                                   ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                         >  1:1892892/1‑91 (MQ=255)
                                     cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                       <  1:2377888/91‑1 (MQ=255)
                                     cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                       >  2:544314/1‑91 (MQ=255)
                                        aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                    <  2:2684491/91‑1 (MQ=255)
                                        aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTCCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                    >  1:3875263/1‑91 (MQ=255)
                                          cTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattt                                                  >  2:3176599/1‑91 (MQ=255)
                                           ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                 >  2:2280563/1‑91 (MQ=255)
                                           ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                 >  1:899304/1‑91 (MQ=255)
                                            tGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                <  2:3821199/91‑1 (MQ=255)
                                                aCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttg                                            >  1:2772794/1‑91 (MQ=255)
                                                aCTAACATAATATTAAAAATATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAGATTTATTTATTttg                                            <  2:1057477/91‑1 (MQ=255)
                                                 cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                           >  2:474902/1‑91 (MQ=255)
                                                   aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                         >  1:1323610/1‑91 (MQ=255)
                                                   aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                         <  1:2605151/91‑1 (MQ=255)
                                                   aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                         <  1:2119366/91‑1 (MQ=255)
                                                   aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                         <  1:23124/91‑1 (MQ=255)
                                                        aataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtt                                    <  1:463841/91‑1 (MQ=255)
                                                        aataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtt                                    <  1:1908640/91‑1 (MQ=255)
                                                         atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                   >  1:3655814/1‑91 (MQ=255)
                                                         atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                   >  2:2646916/1‑91 (MQ=255)
                                                         atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTACTTATTTTGTTGATAttt                                   >  2:1796292/1‑91 (MQ=255)
                                                          tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                  >  1:796028/1‑91 (MQ=255)
                                                          tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                  >  2:60143/1‑91 (MQ=255)
                                                           atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                 >  2:811701/1‑91 (MQ=255)
                                                             tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                               >  1:2424354/1‑91 (MQ=255)
                                                              aaaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                              >  2:3663054/1‑91 (MQ=255)
                                                              aaaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                              >  2:151896/1‑91 (MQ=255)
                                                                aaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTGTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGaa                            <  2:3186863/91‑1 (MQ=255)
                                                                     taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                       >  2:21576/1‑91 (MQ=255)
                                                                     taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                       >  1:3332155/1‑91 (MQ=255)
                                                                     taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                       <  2:2255515/91‑1 (MQ=255)
                                                                      aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                      >  2:3687302/1‑91 (MQ=255)
                                                                       cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                     >  1:65896/1‑91 (MQ=255)
                                                                       cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                     >  1:2227876/1‑91 (MQ=255)
                                                                       cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                     <  1:983055/91‑1 (MQ=255)
                                                                        gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                    >  1:1789531/1‑91 (MQ=255)
                                                                        gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                    >  1:2262889/1‑91 (MQ=255)
                                                                        gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                    >  1:336939/1‑91 (MQ=255)
                                                                         aaaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATa                   >  2:1892727/1‑91 (MQ=255)
                                                                           aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                 >  2:2456888/1‑91 (MQ=255)
                                                                           aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                 >  1:1700177/1‑91 (MQ=255)
                                                                           aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                 <  2:3332155/91‑1 (MQ=255)
                                                                            aaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAAt                <  1:3724313/91‑1 (MQ=255)
                                                                             aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg               <  1:320394/91‑1 (MQ=255)
                                                                              ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa              <  1:1580992/91‑1 (MQ=255)
                                                                              ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATCTATTTATCTTGTTGATATATTAATGAAGCTTCAATAAATGa              >  2:2729692/1‑91 (MQ=255)
                                                                               cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt             >  2:117799/1‑91 (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           >  2:3199040/1‑91 (MQ=255)
                                                                                 gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt           <  2:659314/91‑1 (MQ=255)
                                                                                 gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt           <  1:2634698/91‑1 (MQ=255)
                                                                                 gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt           <  1:668892/91‑1 (MQ=255)
                                                                                  tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa          >  2:3486890/1‑91 (MQ=255)
                                                                                  tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa          >  1:1799275/1‑91 (MQ=255)
                                                                                    cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa        <  1:235871/91‑1 (MQ=255)
                                                                                    cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa        >  1:1093464/1‑91 (MQ=255)
                                                                                    cTTTCAACGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa        <  1:1166274/91‑1 (MQ=255)
                                                                                     tttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaa       <  1:2257721/91‑1 (MQ=255)
                                                                                      ttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaaa      >  2:782992/1‑91 (MQ=255)
                                                                                       tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt     >  2:2051355/1‑91 (MQ=255)
                                                                                       tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt     <  2:1222272/91‑1 (MQ=255)
                                                                                          aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg  >  1:3448301/1‑91 (MQ=255)
                                                                                          aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg  >  2:2521992/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGG  >  USA300TCH1516_ALE/2243404‑2243584

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: