Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 372,106 T→A 14.7% intergenic (‑166/‑204) ylbJ ← / → lip2 Sporulation integral membrane protein YlbJ/Lipase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE372,1060TA14.7% 376.4 / 48.2 128intergenic (‑166/‑204)ylbJ/lip2Sporulation integral membrane protein YlbJ/Lipase 2
Reads supporting (aligned to +/- strand):  ref base T (53/56);  new base A (14/5);  total (67/61)
Fisher's exact test for biased strand distribution p-value = 4.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

ATTTTCATACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGTATTATTAA  >  USA300TCH1516_ALE/372017‑372194
                                                                                         |                                                                                        
attTTCATACCACTACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGAtt                                                                                         >  1:3179227/1‑91 (MQ=255)
 ttttCATACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTc                                                                                        <  1:3249359/91‑1 (MQ=255)
 ttttCATACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTc                                                                                        >  1:1658180/1‑91 (MQ=255)
      aTACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGAATCAacac                                                                                   <  1:287726/91‑1 (MQ=255)
       tACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAAcacc                                                                                  >  1:2971367/1‑91 (MQ=255)
       tACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAAcacc                                                                                  >  1:2971662/1‑91 (MQ=255)
       tACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAAcacc                                                                                  <  1:1442770/91‑1 (MQ=255)
        aCCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAAcacca                                                                                 <  1:933789/91‑1 (MQ=255)
         ccAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTGGTGTTGAATCAacacaac                                                                                <  1:882065/91‑4 (MQ=39)
         ccAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTGGTGTTGAATCAacacaac                                                                                <  1:2693776/91‑4 (MQ=39)
         ccAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTGGTGTTGAATCAacacaac                                                                                <  1:1531476/91‑4 (MQ=39)
          cAAGACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTGGTGTTGAATcaacacaaca                                                                               <  1:3341815/91‑1 (MQ=38)
           aaTACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACaa                                                                              >  1:4117898/1‑91 (MQ=255)
             tACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACt                                                                            <  1:2032002/91‑1 (MQ=255)
             tACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACt                                                                            <  1:1075650/91‑1 (MQ=255)
                aaaTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTaaa                                                                         <  1:2355471/91‑1 (MQ=255)
                aaaTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTaaa                                                                         <  1:2355404/91‑1 (MQ=255)
                 aaTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTaaaa                                                                        >  1:2318183/1‑91 (MQ=255)
                 aaTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTaaaa                                                                        >  1:3096907/1‑91 (MQ=255)
                  aTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAg                                                                       >  1:2429436/1‑91 (MQ=255)
                  aTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAg                                                                       <  1:3456632/91‑1 (MQ=255)
                  aTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAg                                                                       <  1:2866245/91‑1 (MQ=255)
                  aTTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAg                                                                       <  1:2866518/91‑1 (MQ=255)
                   ttAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAgt                                                                      >  1:2486828/1‑91 (MQ=255)
                   ttAACTAATAATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAAGTAAAAgt                                                                      >  1:2079738/1‑91 (MQ=255)
                    tAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAgtg                                                                     <  1:2063043/91‑1 (MQ=255)
                    tAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAgtg                                                                     <  1:518114/91‑1 (MQ=255)
                     aaCTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAgtgt                                                                    >  1:1354290/1‑91 (MQ=255)
                     aaCTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAgtgt                                                                    >  1:1200133/1‑91 (MQ=255)
                        tAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTta                                                                 <  1:2053219/91‑1 (MQ=255)
                        tAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTta                                                                 <  1:1872201/91‑1 (MQ=255)
                        tAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTta                                                                 <  1:1872150/91‑1 (MQ=255)
                         aaTTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTtat                                                                >  1:2744264/1‑91 (MQ=255)
                           ttATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATaa                                                              <  1:4056366/91‑1 (MQ=255)
                              tatataGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAatt                                                           <  1:236981/91‑1 (MQ=255)
                                 ataGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAattatt                                                        >  1:2582307/1‑91 (MQ=255)
                                 ataGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAattatt                                                        >  1:2582330/1‑91 (MQ=255)
                                   aGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCGACACCACAACTAAAAGTGTTTATAAATTATTTg                                                      >  1:3685493/1‑91 (MQ=255)
                                   aGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAGTTATTTg                                                      >  1:2067361/1‑91 (MQ=255)
                                   aGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTg                                                      >  1:2067362/1‑91 (MQ=255)
                                         aaaCTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATa                                                <  1:1001619/91‑1 (MQ=255)
                                          aaCTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATac                                               <  1:141397/91‑1 (MQ=255)
                                           aCTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATaca                                              <  1:1875900/91‑1 (MQ=255)
                                           aCTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATaca                                              >  1:221901/1‑91 (MQ=255)
                                           aCTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATaca                                              <  1:1974202/91‑1 (MQ=255)
                                           aCTATATTACTTAATAAAATATTTATCTTAAAAGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATaca                                              <  1:1189405/91‑1 (MQ=255)
                                              atatTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACat                                           >  1:3887687/1‑91 (MQ=255)
                                               tatTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACata                                          >  1:34395/1‑91 (MQ=255)
                                               tatTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACata                                          <  1:3698820/91‑1 (MQ=255)
                                               tatTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACata                                          >  1:4025646/1‑91 (MQ=255)
                                                 ttACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATAtt                                        <  1:1621066/91‑1 (MQ=255)
                                                 ttACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATAtt                                        <  1:1714271/91‑1 (MQ=255)
                                                 ttACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATAtt                                        >  1:783374/1‑91 (MQ=255)
                                                  tACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATAttt                                       <  1:1329089/91‑1 (MQ=255)
                                                  tACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATAttt                                       >  1:2489102/1‑91 (MQ=255)
                                                   aCTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTg                                      <  1:102752/91‑1 (MQ=255)
                                                   aCTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTg                                      <  1:2249044/91‑1 (MQ=255)
                                                        ataaAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAAt                                 <  1:2644675/91‑1 (MQ=255)
                                                        ataaAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAAt                                 <  1:2644313/91‑1 (MQ=255)
                                                        ataaAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAAt                                 >  1:3170654/1‑91 (MQ=255)
                                                         taaAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATg                                >  1:3893081/1‑91 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTTTGTTGATTCAACACCACAACTAAAAGTCTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            <  1:3597666/91‑1 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            >  1:1050942/1‑91 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            >  1:1621421/1‑91 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            >  1:3275650/1‑91 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            <  1:1957937/91‑1 (MQ=255)
                                                             atatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            >  1:2005935/1‑91 (MQ=255)
                                                             agatTTATCTTAAAGGTTGTTTTGTTGATTCAACTCCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTa                            <  1:3743274/89‑1 (MQ=255)
                                                              tatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAg                           >  1:128049/1‑91 (MQ=255)
                                                              tatTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAg                           >  1:2452080/1‑91 (MQ=255)
                                                               atTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGt                          >  1:1538953/1‑91 (MQ=255)
                                                               atTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGt                          >  1:4097743/1‑91 (MQ=255)
                                                                tttATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTa                         >  1:3901619/1‑91 (MQ=255)
                                                                   aTCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCg                      >  1:153892/1‑91 (MQ=255)
                                                                       tAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGAttt                  <  1:1710967/91‑1 (MQ=255)
                                                                        tttagttgtggtGTTGAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTa                 >  1:859826/5‑91 (MQ=37)
                                                                        aaaTGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTa                 <  1:1636107/91‑1 (MQ=255)
                                                                          aTGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAat               <  1:2746487/91‑1 (MQ=255)
                                                                          aTGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAat               <  1:1103042/91‑1 (MQ=255)
                                                                           agttgtggtGTTGAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAata              >  1:2795793/2‑91 (MQ=38)
                                                                           agttgtggtGTTGAAACAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATCTAata              >  1:2795881/2‑91 (MQ=38)
                                                                             ttgttgtGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATc            >  1:1305221/1‑91 (MQ=255)
                                                                             ttgttgtGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATc            >  1:1305725/1‑91 (MQ=255)
                                                                               gttgtgttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGt          >  1:1403367/1‑91 (MQ=255)
                                                                               gttgtgttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGt          >  1:3562888/1‑91 (MQ=255)
                                                                               gttgtgttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGt          >  1:3562199/1‑91 (MQ=255)
                                                                               gttgtgttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGt          <  1:2099791/91‑1 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAACGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGta         >  1:2063279/3‑91 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATGTCGta         >  1:133646/3‑91 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGta         >  1:907046/3‑91 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGta         >  1:133647/3‑91 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGta         >  1:3548663/3‑91 (MQ=255)
                                                                                tggtgttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGta         >  1:1704866/3‑91 (MQ=255)
                                                                                    gttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGtattat     <  1:2566538/91‑1 (MQ=255)
                                                                                    gttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGtattat     <  1:803768/91‑1 (MQ=255)
                                                                                    gttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGAGTAGTATCGATTTAATATCGtattat     >  1:3875724/1‑91 (MQ=255)
                                                                                    gttgATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGACTAGTATCGATTTAATATCGtattat     >  1:3876020/1‑91 (MQ=255)
                                                                                    gttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGtattat     >  1:1302547/1‑91 (MQ=255)
                                                                                     ttgAATCAACACAACAACTAAAAGTGTTTATAAATTATTTGGAAATAGACATATTTGTAAATGATTAGTATCGATTTGATATCGTattatt    >  1:4150728/1‑91 (MQ=255)
                                                                                       gATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGTATTATTaa  >  1:2015685/1‑91 (MQ=255)
                                                                                         |                                                                                        
ATTTTCATACCAATACAAATTAACTAATTATATATAGATTGAAACTATATTACTTAATAAAATATTTATCTTAAATGTTGTTGTGTTGATTCAACACCACAACTAAAAGTGTTTATAAATTATTTGGAAATACACATATTTGTAAATGATTAGTATCGATTTAATATCGTATTATTAA  >  USA300TCH1516_ALE/372017‑372194

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