Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,148 C→T 5.4% intergenic (‑160/+60) 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,1480CT 5.4% 795.7 / 22.5 242intergenic (‑160/+60)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base C (119/109);  new base T (6/7);  total (126/116)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

Only 100 of 249 total aligned reads displayed.

AATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685058‑1685238
                                                                                          |                                                                                          
 atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                           <  2:1560892/91‑1 (MQ=255)
   atCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtg                                                                                         <  2:3320304/91‑1 (MQ=255)
    tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        <  2:3154956/91‑1 (MQ=255)
    tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        >  1:2229805/1‑91 (MQ=255)
     cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  1:3701864/91‑1 (MQ=255)
     cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  2:661214/91‑1 (MQ=255)
      aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  2:1610532/91‑1 (MQ=255)
      aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  1:1524104/91‑1 (MQ=255)
      aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  2:3606118/91‑1 (MQ=255)
       aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTAATATGTGAt                                                                                     <  1:168734/91‑1 (MQ=18)
       aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                                     <  1:311062/91‑1 (MQ=255)
         ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   >  1:1004152/1‑91 (MQ=255)
         ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   >  1:1969807/1‑91 (MQ=255)
         ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                                   <  1:3530640/91‑1 (MQ=21)
          tGATTATTTATAAACTATATAACTTTTTTTAATTAAAGATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                                                  <  1:3160817/91‑1 (MQ=255)
           gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                 >  1:2494446/1‑91 (MQ=255)
           gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                 <  1:482014/91‑1 (MQ=255)
           gATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgttt                                                                                 >  2:2632965/1‑91 (MQ=255)
              tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                              <  1:3037456/91‑1 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAATTTGACAACAAAAAAACATCACATATTATTGTCAGACACTAACTTGTGATgtttgttt                                                                             >  2:1465222/1‑91 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             <  1:1761485/91‑1 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             >  1:1471244/1‑91 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             <  1:2604457/91‑1 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             <  1:40549/91‑1 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             >  1:1336317/1‑91 (MQ=255)
                tttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGtttt                                                                            >  1:2846581/1‑91 (MQ=255)
                 ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                           >  1:3332636/1‑91 (MQ=255)
                 ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                           >  1:747679/1‑91 (MQ=255)
                  tataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtat                                                                          >  2:2160319/1‑91 (MQ=255)
                  tataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtat                                                                          <  2:1339601/91‑1 (MQ=255)
                   ataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtata                                                                         <  1:2295905/91‑1 (MQ=255)
                      aaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaa                                                                      <  2:2466009/91‑1 (MQ=255)
                       aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                                     <  1:2839950/91‑1 (MQ=255)
                        cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                    <  2:2078166/91‑1 (MQ=255)
                        cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                    >  2:155212/1‑91 (MQ=255)
                           tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                                 <  1:1078943/91‑1 (MQ=255)
                           tataACTTTTTGTAAGTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                                 <  2:632898/91‑1 (MQ=255)
                            ataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAata                                                                >  2:2701097/1‑91 (MQ=255)
                             taACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAatat                                                               <  1:970499/91‑1 (MQ=255)
                              aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                              >  2:1271821/1‑91 (MQ=255)
                               aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             >  2:587565/1‑91 (MQ=255)
                                  ttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtg                                                          <  1:1788346/91‑1 (MQ=255)
                                   tttttAATTAACAATTGACAACAGAAAAACATCACATATTAGTGTCAGACACTAACTTGGGATGTTTGTTTTATATAATTAATATTCgtgt                                                         <  2:3016687/91‑1 (MQ=255)
                                    ttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGtt                                                        >  2:3709984/1‑91 (MQ=255)
                                      ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                                      <  2:391664/91‑1 (MQ=255)
                                       taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                     >  1:2931277/1‑91 (MQ=255)
                                       taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                     >  1:1906523/1‑91 (MQ=255)
                                         atttaAATTTCACAACAAAAAAAAATCCCATATTAGTGTCAGTCACTAACTTGTTATGTTTGTCTTATATAATTAATATTCATCTTGaagt                                                   >  2:3804862/1‑88 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   >  1:539672/1‑91 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   >  1:539604/1‑91 (MQ=255)
                                          ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                                  >  2:3468807/1‑91 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 >  2:470086/1‑91 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  2:1077251/91‑1 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                               <  1:1855671/91‑1 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                               <  1:682660/91‑1 (MQ=255)
                                               aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACCAACTTGTGATGTTTGTTTTATATAATCAATATTCGTGTTGAGGGAGTCaa                                             >  2:3675961/1‑91 (MQ=255)
                                                aTTGACAACAAAAACACATCCCATATTAGTGTCAGAGAATAACTTGTGCTGGTTGTTTTATATAATTAATATTCGTGTTGAGGGCGTCacg                                            >  2:681563/1‑89 (MQ=255)
                                                   gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                         >  1:2329550/1‑91 (MQ=255)
                                                   gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                         >  1:48043/1‑91 (MQ=255)
                                                    acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                        >  1:1888831/1‑91 (MQ=255)
                                                    acaacaAAAAAACATCACATATTAGTGTCAGACACTAAATTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                        >  2:3849545/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCGACAGAgg                                       >  2:3503267/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                       >  2:1784211/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                       >  2:294027/1‑91 (MQ=255)
                                                       acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                                     <  2:3761578/91‑1 (MQ=255)
                                                        caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                    <  2:2094421/91‑1 (MQ=255)
                                                        caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                    >  1:2461953/1‑91 (MQ=255)
                                                         aaaaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatt                                   <  2:549991/91‑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)
                                                             aaaCATCACATATTAGTATCAGACACTAACTTATGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                               <  1:3384997/91‑1 (MQ=255)
                                                               aCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattat                             <  1:3287229/91‑1 (MQ=255)
                                                               aCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattat                             >  1:2835589/1‑91 (MQ=255)
                                                                cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                            <  1:782561/91‑1 (MQ=255)
                                                                cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                            <  1:3610472/91‑1 (MQ=255)
                                                                   cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                         >  1:2195007/1‑91 (MQ=255)
                                                                   cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                         >  1:2746810/1‑91 (MQ=255)
                                                                     caTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                       >  1:511298/1‑91 (MQ=255)
                                                                      atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      >  2:2495612/1‑91 (MQ=255)
                                                                      atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      <  2:1938432/91‑1 (MQ=255)
                                                                      aagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      >  1:863656/4‑91 (MQ=37)
                                                                        atTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttta                    <  2:3622237/91‑1 (MQ=255)
                                                                         ttAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttac                   <  2:3373723/91‑1 (MQ=255)
                                                                          tAGTGTCAGACACTAGCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGGGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  2:1699437/91‑1 (MQ=255)
                                                                           aGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgt                 <  2:277245/91‑1 (MQ=255)
                                                                           aGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgt                 >  1:2791222/1‑91 (MQ=255)
                                                                             tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt               >  2:2944562/1‑91 (MQ=255)
                                                                              gtctgACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttta              >  2:761111/5‑91 (MQ=255)
                                                                               tctgACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttac             >  1:640588/4‑91 (MQ=255)
                                                                                cAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttacg            <  1:2884426/91‑1 (MQ=255)
                                                                                 aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt           >  2:2254526/1‑91 (MQ=255)
                                                                                 aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt           <  2:2851990/91‑1 (MQ=255)
                                                                                 aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt           <  2:814373/91‑1 (MQ=255)
                                                                                     acTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAg       >  1:3887844/1‑91 (MQ=255)
                                                                                       tAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCt     >  1:1135499/1‑91 (MQ=255)
                                                                                         aCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTc   >  1:2126148/1‑91 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt  >  1:2919441/1‑91 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt  <  2:288318/91‑1 (MQ=255)
                                                                                          |                                                                                          
AATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685058‑1685238

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: