Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,093 T→A 17.0% intergenic (+87/‑275) yugI_2 → / → USA300HOU_RS04740 General stress protein 13/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE950,0930TA17.0% 749.6 / 124.9 274intergenic (+87/‑275)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (118/107);  new base A (20/26);  total (140/134)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 291 total aligned reads displayed.

CTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAT  >  USA300TCH1516_ALE/950003‑950183
                                                                                          |                                                                                          
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            >  2:2815833/1‑91 (MQ=255)
 tAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtt                                                                                           <  2:356549/91‑1 (MQ=255)
  aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                          >  1:1391497/1‑91 (MQ=255)
  aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                          >  1:1842287/1‑91 (MQ=255)
   aaGGGACAGATCGATCGTACCGAAGATCATACAATGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGGCTTGAGTATCTATTAAAtttt                                                                                         <  2:2269735/91‑1 (MQ=255)
    aGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTtta                                                                                        <  1:2272331/91‑1 (MQ=255)
        aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                    <  2:940310/91‑1 (MQ=38)
         cAGATGAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGa                                                                                   <  2:3546731/91‑1 (MQ=37)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                >  2:2351349/1‑91 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                                <  1:2795342/91‑1 (MQ=37)
             tAAGTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                               <  1:3191180/91‑1 (MQ=255)
             tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACt                                                                               <  2:3743800/91‑1 (MQ=37)
             tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACt                                                                               <  1:186141/91‑1 (MQ=37)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  1:1030183/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             >  2:1843314/1‑91 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  2:2297357/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  1:1637243/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  2:606262/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTaa                                                                             <  2:3604519/91‑1 (MQ=25)
                aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                            <  2:3061701/91‑1 (MQ=255)
                 tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                           <  1:1356330/91‑1 (MQ=255)
                   gTACCGAAAATCCTGCCAAGGTTCTGAATTGAAAGTTTCTTAGACTATAAAAGCGAGTAGTTTCTAATAAAATTTAATAGATACTAAtctc                                                                         <  2:2176738/91‑1 (MQ=16)
                   gTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtctc                                                                         >  1:373983/1‑91 (MQ=255)
                   gTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtctc                                                                         >  2:3656431/1‑91 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        >  1:47940/1‑91 (MQ=255)
                     aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                       >  2:686534/1‑91 (MQ=255)
                     aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                       >  1:3257363/1‑91 (MQ=255)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                      >  2:271342/1‑91 (MQ=255)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                      <  2:2069608/91‑1 (MQ=255)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATAAAATTTTATTAGATACTAATCTCttt                                                                      >  1:1824891/1‑91 (MQ=255)
                       cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                     <  1:2702623/91‑1 (MQ=255)
                       cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                     <  2:2032340/91‑1 (MQ=255)
                        gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                                    <  1:3103521/91‑1 (MQ=255)
                         aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTg                                                                   >  2:2615126/1‑91 (MQ=255)
                         aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTg                                                                   >  2:1596935/1‑91 (MQ=255)
                          aaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGt                                                                  <  1:1246255/91‑1 (MQ=255)
                             tCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTa                                                               <  2:1659980/91‑1 (MQ=255)
                             tCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTa                                                               >  1:1498625/1‑91 (MQ=255)
                                tACAAATGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                            <  2:3716714/91‑1 (MQ=255)
                                tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                            >  1:1692253/1‑91 (MQ=255)
                                 aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                           <  2:3041516/91‑1 (MQ=255)
                                    aaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                        <  1:249575/91‑1 (MQ=255)
                                    aaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                        >  1:3486911/1‑91 (MQ=255)
                                     aGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACg                                                       <  2:3019694/91‑1 (MQ=255)
                                      gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                      <  1:2298608/91‑1 (MQ=255)
                                        gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                                    >  1:2392177/1‑91 (MQ=255)
                                         tCTGATATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                   <  2:3646942/91‑1 (MQ=255)
                                         tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                   >  1:2192226/1‑91 (MQ=255)
                                         tCTGAAACGAAAGTTTCTTAGACTATAAAAGAGACTACCACCTATCAAACTTCACTACACACTAATCTCTTCCTGTCTACGATAACGTAat                                                   >  2:3944264/1‑91 (MQ=255)
                                          cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                  <  1:1171668/91‑1 (MQ=255)
                                               atgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgatt                                             <  2:1714906/91‑1 (MQ=255)
                                                  aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                          <  1:763839/91‑1 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        <  1:3303313/91‑1 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        >  1:153393/1‑91 (MQ=255)
                                                      tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                      >  1:2250994/1‑91 (MQ=255)
                                                        tCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAtt                                    >  1:2755647/1‑91 (MQ=255)
                                                         cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                                   >  1:1441711/1‑91 (MQ=255)
                                                            aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                >  2:2897913/1‑91 (MQ=255)
                                                            aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                >  1:3027398/1‑91 (MQ=255)
                                                            aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                >  1:2029914/1‑91 (MQ=255)
                                                               cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                             >  1:105146/1‑91 (MQ=255)
                                                               cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                             <  1:2721650/91‑1 (MQ=255)
                                                               cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGACTACGATAACGTAATATGATTGATTCTATTTACACGt                             <  1:3343800/91‑1 (MQ=255)
                                                                 ataAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTAc                           >  1:1426304/1‑91 (MQ=255)
                                                                  taAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                          >  1:337988/1‑91 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         <  2:3130369/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  2:2896448/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCAACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  2:1653293/91‑1 (MQ=255)
                                                                      agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      >  1:1331168/1‑91 (MQ=255)
                                                                      agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      <  2:1717973/91‑1 (MQ=255)
                                                                      agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      <  2:2956574/91‑1 (MQ=255)
                                                                      agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      <  2:1530098/91‑1 (MQ=255)
                                                                        agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                    <  1:960964/91‑1 (MQ=255)
                                                                         gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   >  2:1300134/1‑91 (MQ=255)
                                                                          aTTAGTATCGAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAGATGGttt                  >  1:1969253/1‑91 (MQ=25)
                                                                           ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                 >  2:10303/1‑91 (MQ=255)
                                                                             aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:2513135/1‑91 (MQ=255)
                                                                             aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:360522/1‑91 (MQ=255)
                                                                             aGTATCTAATAAAATTTAATAAATACTAATCTCTTTTATTCTACGATCACGTAATATGATTGATCCTATTTACACCTACAAATGGTTTAAg               >  2:3422951/1‑91 (MQ=18)
                                                                              gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg              >  2:2530888/1‑91 (MQ=255)
                                                                              gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg              >  1:2989510/1‑91 (MQ=255)
                                                                              gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg              <  2:3260708/91‑1 (MQ=255)
                                                                                aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            >  2:788314/1‑91 (MQ=255)
                                                                                aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            >  2:1347718/1‑91 (MQ=255)
                                                                                aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            >  2:35870/1‑91 (MQ=255)
                                                                                aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATACTATTTACACGTACAAATGGTTTAAGGTg            >  1:3167410/1‑91 (MQ=255)
                                                                                aTCTAGTAATTTTTAGTAGATACTGATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            <  1:2339595/91‑1 (MQ=255)
                                                                                 tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           >  2:3455040/1‑91 (MQ=255)
                                                                                 tCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAATGGTTTAAGGTGAc          >  1:55292/1‑91 (MQ=37)
                                                                                 ccTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           >  2:3128212/2‑91 (MQ=255)
                                                                                   taatAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGATGACa         >  1:3685180/4‑91 (MQ=38)
                                                                                   tATTAAATTTTATTAGATACTAATCTCTTTTTTTTTATTTTTACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa         >  2:2404500/1‑91 (MQ=255)
                                                                                     ttAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata       <  1:3675858/91‑1 (MQ=255)
                                                                                     ttAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata       <  1:1778700/91‑1 (MQ=255)
                                                                                       aaaCTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTCACATATc     >  2:3507838/1‑91 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  2:857638/91‑1 (MQ=255)
                                                                                         aTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCa   >  1:2592781/1‑91 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  2:2132592/91‑1 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  2:2357367/91‑1 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  1:617625/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAT  >  USA300TCH1516_ALE/950003‑950183

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: