Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,139 A→T 7.5% intergenic (‑151/+69) 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,1390AT 7.5% 840.7 / 43.8 253intergenic (‑151/+69)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base A (115/119);  new base T (8/11);  total (123/130)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 258 total aligned reads displayed.

AACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGT  >  USA300TCH1516_ALE/1685049‑1685229
                                                                                          |                                                                                          
  ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                          >  1:2651931/1‑91 (MQ=255)
     ataGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGacac                                                                                       <  1:3752763/91‑1 (MQ=255)
     ataGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGacac                                                                                       <  1:3659104/91‑1 (MQ=255)
      taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      <  1:3275869/91‑1 (MQ=255)
      taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      <  1:1327741/91‑1 (MQ=255)
      taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      <  1:142659/91‑1 (MQ=255)
      taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      <  1:2274567/91‑1 (MQ=255)
      taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      >  1:3533270/1‑91 (MQ=255)
       aGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTa                                                                                     <  2:48674/91‑1 (MQ=255)
       aGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTa                                                                                     <  1:246040/91‑1 (MQ=38)
       aGAATATCAATTTGATTATTTATAAACTATATAACTTATTTTAATTAAAAATTGACAACAAAAAAACATCACATTTTAGTGTCAGACACTa                                                                                     >  1:3319473/1‑91 (MQ=255)
        gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                                    >  2:1092529/1‑91 (MQ=255)
         aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                                   >  1:518468/1‑91 (MQ=255)
         aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTaat                                                                                   <  2:2947948/91‑2 (MQ=38)
           tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                 <  2:1840458/91‑1 (MQ=255)
           tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                 >  2:1041306/1‑91 (MQ=255)
            atCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtg                                                                                <  2:3320304/91‑1 (MQ=255)
            atCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtg                                                                                >  1:687765/1‑91 (MQ=255)
             tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTCAAAATTCACAACAAAAAAACATCACTCATTAGTGTCAGACACTAACTtgt                                                                               >  2:3673707/1‑91 (MQ=255)
              cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                              <  2:2819920/91‑1 (MQ=255)
              cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                              >  1:1288328/1‑91 (MQ=255)
               aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAGCTTGTGa                                                                             >  1:1962146/1‑91 (MQ=255)
               aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                             <  1:2256369/91‑1 (MQ=255)
               aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                             <  2:3606118/91‑1 (MQ=255)
               aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                             <  1:901884/91‑1 (MQ=255)
               aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                             <  2:1610532/91‑1 (MQ=255)
                aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                            <  1:2732059/91‑1 (MQ=255)
                aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAt                                                                            <  1:28056/91‑1 (MQ=25)
                aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAt                                                                            <  1:28411/91‑1 (MQ=25)
                 tttGATTATTTATAAACTATATAACTATATTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtg                                                                           <  2:2682245/91‑1 (MQ=255)
                  ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTTACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTGTTGAtgt                                                                          >  2:3505548/1‑91 (MQ=255)
                  ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                          >  1:1004152/1‑91 (MQ=255)
                  ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                          <  2:2113593/91‑1 (MQ=255)
                  ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                          <  1:3530640/91‑1 (MQ=21)
                    gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                        <  1:482014/91‑1 (MQ=255)
                    gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                        >  2:897119/1‑91 (MQ=255)
                    gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                        <  1:75825/91‑1 (MQ=255)
                     attattTATAAACTATATAACTTTTTTTAATTAAAAATAGAAAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttg                                                                       <  2:725178/91‑1 (MQ=255)
                        atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                    <  1:1761485/91‑1 (MQ=255)
                          ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                  <  1:219511/91‑1 (MQ=255)
                          ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                  >  1:3332636/1‑91 (MQ=255)
                           tataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtat                                                                 <  2:1339601/91‑1 (MQ=255)
                              aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                              <  1:1804624/91‑1 (MQ=255)
                               aaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaa                                                             <  2:2466009/91‑1 (MQ=255)
                                 cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                           <  1:3906328/91‑1 (MQ=255)
                                 cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                           >  2:155212/1‑91 (MQ=255)
                                   atataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaa                                                         <  1:21597/91‑1 (MQ=255)
                                    tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                        <  1:1078943/91‑1 (MQ=255)
                                    tataACTTTTTGTAAGTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                        <  2:632898/91‑1 (MQ=255)
                                       aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                     >  2:1271821/1‑91 (MQ=255)
                                        aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                    <  2:2924458/91‑1 (MQ=255)
                                            tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                                <  2:3141031/91‑1 (MQ=255)
                                             ttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGtt                                               >  2:3709984/1‑91 (MQ=255)
                                               ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                             >  2:3440303/1‑91 (MQ=255)
                                               ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                             <  2:1593787/91‑1 (MQ=255)
                                                taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                            >  2:2871742/1‑91 (MQ=255)
                                                taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                            >  1:246498/1‑91 (MQ=255)
                                                  attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTATCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                          <  1:257916/91‑1 (MQ=255)
                                                  attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                          >  1:539604/1‑91 (MQ=255)
                                                  attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                          <  1:3071319/91‑1 (MQ=255)
                                                   ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                         >  2:3468807/1‑91 (MQ=255)
                                                    taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                        <  1:542580/91‑1 (MQ=255)
                                                      aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                      <  1:1857036/91‑1 (MQ=255)
                                                        aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCaa                                    <  2:1480311/91‑1 (MQ=255)
                                                         aTTGACAACAAAAACACATCCCATATTAGTGTCAGAGAATAACTTGTGCTGGTTGTTTTATATAATTAATATTCGTGTTGAGGGCGTCacg                                   >  2:681563/1‑89 (MQ=255)
                                                            gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                >  1:2901546/1‑91 (MQ=255)
                                                            gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                >  1:2329550/1‑91 (MQ=255)
                                                            gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                >  1:48043/1‑91 (MQ=255)
                                                             acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                               >  2:3162165/1‑91 (MQ=255)
                                                             acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                               >  1:1888842/1‑91 (MQ=255)
                                                             acaacaAAAAAACATCACATATTAGTGTCAGACACTAAATTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                               >  2:3849545/1‑91 (MQ=255)
                                                              caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                              >  2:294056/1‑91 (MQ=255)
                                                              caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                              >  2:2230962/1‑91 (MQ=255)
                                                               aacaaAAAAACATCACATATTAGTGTCAGACACGAACATGTGATGTCTGTTTTATATAACTAATATTCGTGTTGAGGGAGTCAACAGAGGt                             >  1:3153515/1‑91 (MQ=255)
                                                                acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                            <  2:3761578/91‑1 (MQ=255)
                                                                 caaAAAAACATCTCATATTAGTGTCATACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTCAGGGAGTCAACAGAGGTat                           <  2:381141/91‑1 (MQ=255)
                                                                 caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                           <  2:2094343/91‑1 (MQ=255)
                                                                    aaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaa                        <  2:644186/91‑1 (MQ=255)
                                                                     caaaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                       >  1:2029650/2‑91 (MQ=21)
                                                                     aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                       <  2:1384913/91‑1 (MQ=255)
                                                                      aaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                      >  2:2147193/1‑91 (MQ=255)
                                                                      aaaCATCACATATTAGTATCAGACACTAACTTATGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                      <  1:3384997/91‑1 (MQ=255)
                                                                       aaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatta                     >  2:2058158/1‑91 (MQ=255)
                                                                        aCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattat                    <  1:3287229/91‑1 (MQ=255)
                                                                         cATCACAAGTTAGTGTCTGACACTAACTTGTGCTGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                   >  1:1671168/1‑91 (MQ=38)
                                                                          aTCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTAttt                  >  1:3018578/1‑91 (MQ=255)
                                                                           tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  1:2838408/1‑91 (MQ=255)
                                                                           tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  1:861240/1‑91 (MQ=255)
                                                                            cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                <  2:3958851/91‑1 (MQ=255)
                                                                            cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                <  1:150050/91‑1 (MQ=255)
                                                                            cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                >  1:2746810/1‑91 (MQ=255)
                                                                               atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt             >  2:2495612/1‑91 (MQ=255)
                                                                                  ttAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttac          <  2:2707266/91‑1 (MQ=255)
                                                                                   tAGTGTCAGACACTAGCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGGGTCAACAGAGGTATTAATTATTTGAAtttacg         <  2:1699437/91‑1 (MQ=255)
                                                                                    aGTGTCAGACACTATCTTGTGATGTTTGTATTATATAATTAATATTCGGGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttgcat        <  2:570701/91‑5 (MQ=255)
                                                                                     gtgtCTGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtt       >  1:1994206/1‑91 (MQ=255)
                                                                                        tCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttac    >  2:2562941/1‑91 (MQ=255)
                                                                                         cAGACGCTAACTTGTGATGGTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttacg   <  2:1721516/91‑1 (MQ=255)
                                                                                         cAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttacg   <  1:2884426/91‑1 (MQ=255)
                                                                                          aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt  >  2:2385439/1‑91 (MQ=255)
                                                                                          |                                                                                          
AACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGT  >  USA300TCH1516_ALE/1685049‑1685229

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: