Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,149 T→A 6.8% intergenic (‑161/+59) USA300HOU_RS08395 ← / ← rpsU hypothetical protein/30S ribosomal protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,685,1490TA 6.8% 788.4 / 32.1 242intergenic (‑161/+59)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base T (117/107);  new base A (11/7);  total (128/114)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTG  >  USA300TCH1516_ALE/1685059‑1685239
                                                                                           |                                                                                          
atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                            <  2:1976693/91‑1 (MQ=255)
   tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                         <  2:703643/91‑1 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                        <  2:2819920/91‑1 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                        <  2:114373/91‑1 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                        <  2:662420/91‑1 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTGTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                        <  2:602284/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  1:901884/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  2:1610532/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  1:2256369/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  1:1524104/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACGTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  2:2197442/91‑1 (MQ=255)
      aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                                      <  2:2479799/91‑1 (MQ=255)
       tttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtg                                                                                     >  2:3324859/1‑91 (MQ=255)
       tttAATTATTTATAAACTATATAACTTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAt                                                                                      <  1:501250/91‑1 (MQ=18)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTTACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTGTTGAtgt                                                                                    >  2:3505548/1‑91 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                    >  1:1004152/1‑91 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                                    <  1:3530824/91‑1 (MQ=21)
         tGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                                                   >  1:389670/1‑91 (MQ=255)
          gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                  >  2:2633012/1‑91 (MQ=255)
          gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                  <  1:2341764/91‑1 (MQ=255)
          gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                  >  1:2494446/1‑91 (MQ=255)
              atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                              <  1:1761485/91‑1 (MQ=255)
              atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                              >  1:1336317/1‑91 (MQ=255)
                 tataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtat                                                                           <  2:1339710/91‑1 (MQ=255)
                  ataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtata                                                                          <  1:2295905/91‑1 (MQ=255)
                  ataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtata                                                                          >  2:521662/1‑91 (MQ=255)
                   taAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatat                                                                         >  1:3257931/1‑91 (MQ=255)
                    aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                                        <  1:1804624/91‑1 (MQ=255)
                     acgtatataACTTATTTTAATAAAAAATTGACAACAAAAAAACCTCACATATTAGTGTCAGACCCGAACTTGTGATGTTTGTTTTATAtaa                                                                       <  2:3855009/88‑1 (MQ=255)
                      aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                                      <  1:2839950/91‑1 (MQ=255)
                      aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                                      >  1:534901/1‑91 (MQ=255)
                       cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                     >  2:155212/1‑91 (MQ=255)
                       cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                     <  2:2078166/91‑1 (MQ=255)
                          tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                                  <  1:1078943/91‑1 (MQ=255)
                           ataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAata                                                                 >  2:2701097/1‑91 (MQ=255)
                            taACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACATATTAGTGTCAGACACTAACACGTGATGTTTGTCTTATATAATTAatat                                                                >  2:3468845/1‑91 (MQ=255)
                            taACTTTTTTTAATTAAAAATTGACAACAAAATAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAatat                                                                <  1:3355120/91‑1 (MQ=255)
                            taACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAatat                                                                <  1:970499/91‑1 (MQ=255)
                             aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                               >  2:1271821/1‑91 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                              <  2:2924458/91‑1 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                              <  1:2129854/91‑1 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                              <  1:2845596/91‑1 (MQ=255)
                                   tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                                          <  2:3141031/91‑1 (MQ=255)
                                      ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                                       >  2:3440303/1‑91 (MQ=255)
                                      ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                                       <  2:1593787/91‑1 (MQ=255)
                                       taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                      >  1:2833302/1‑91 (MQ=255)
                                       taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                      >  1:1906523/1‑91 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTATCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    <  1:257916/91‑1 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    >  1:539672/1‑91 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    >  2:3804745/1‑91 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    <  1:3071319/91‑1 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    >  1:3800535/1‑91 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                  <  1:2975196/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTAGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                  <  2:2997299/91‑1 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                                <  1:682660/91‑1 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                                >  1:3111810/1‑91 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                                <  1:1855671/91‑1 (MQ=255)
                                               aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCaa                                              <  2:1480311/91‑1 (MQ=255)
                                               aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACAAACTTGCGATGTTTGTTTTATATAACCAATATTCGTGTTGAGGGAGTCaa                                              >  2:3676601/1‑91 (MQ=255)
                                                   gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                          >  1:48096/1‑91 (MQ=255)
                                                    acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                         >  1:1888842/1‑91 (MQ=255)
                                                    acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                         >  2:3162165/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                        >  2:1784211/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                        >  1:2738266/1‑91 (MQ=255)
                                                       acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                                      >  2:283182/1‑91 (MQ=255)
                                                        caaAAAAACATCTCATATTAGTGTCATACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTCAGGGAGTCAACAGAGGTat                                     <  2:381141/91‑1 (MQ=255)
                                                        caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                     >  1:331138/1‑91 (MQ=255)
                                                          aacaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatta                                   <  1:2899370/88‑1 (MQ=255)
                                                            caaaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                                 >  1:2029650/2‑91 (MQ=21)
                                                            aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                                 <  2:1384913/91‑1 (MQ=255)
                                                             aaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                                >  2:2147173/1‑91 (MQ=255)
                                                             aaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                                >  1:2711709/1‑91 (MQ=255)
                                                             aaaCATCACATATTAGTCTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                                <  1:3385078/91‑1 (MQ=255)
                                                                cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             <  1:782561/91‑1 (MQ=255)
                                                                cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             <  1:3610514/91‑1 (MQ=255)
                                                                   cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                          <  1:150050/91‑1 (MQ=255)
                                                                   cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                          >  1:2195007/1‑91 (MQ=255)
                                                                    acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                         >  2:2106914/1‑91 (MQ=255)
                                                                      atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                       >  2:2495612/1‑91 (MQ=255)
                                                                      atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTGGAAtt                       >  2:3658561/1‑91 (MQ=255)
                                                                      aagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                       >  2:635800/4‑91 (MQ=37)
                                                                           aGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgt                  >  1:2791222/1‑91 (MQ=255)
                                                                            gtgtCTGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtt                 >  1:1994206/1‑91 (MQ=255)
                                                                             tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt                <  1:71455/91‑1 (MQ=255)
                                                                             tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt                >  2:2944562/1‑91 (MQ=255)
                                                                                cAGACGCTAACTTGTGATGGTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttacg             <  2:1721516/91‑1 (MQ=255)
                                                                                 aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            <  2:814373/91‑1 (MQ=255)
                                                                                 aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            >  2:2254526/1‑91 (MQ=255)
                                                                                   acacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGc          >  1:631745/1‑91 (MQ=255)
                                                                                    cacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCa         <  1:2928345/91‑1 (MQ=255)
                                                                                    cacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCa         >  1:117547/1‑91 (MQ=255)
                                                                                     acTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATAAATTATTTGAATTTACGTTTACGTGCAg        <  2:3621244/91‑1 (MQ=255)
                                                                                      cTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAgg       >  2:3927178/1‑90 (MQ=255)
                                                                                      cTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGc       <  1:414559/91‑1 (MQ=255)
                                                                                        aaCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCtt     >  2:3727041/1‑91 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt   <  2:288318/91‑1 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt   >  1:1695209/1‑91 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt   <  2:566155/91‑1 (MQ=255)
                                                                                           ttGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTg  >  1:2433608/1‑91 (MQ=255)
                                                                                           ttGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTg  >  1:287243/1‑91 (MQ=255)
                                                                                           |                                                                                          
ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTG  >  USA300TCH1516_ALE/1685059‑1685239

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