Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,088 T→A 16.1% intergenic (+82/‑280) 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,0880TA16.1% 768.4 / 117.0 273intergenic (+82/‑280)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (118/110);  new base A (18/26);  total (136/137)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 282 total aligned reads displayed.

AACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATA  >  USA300TCH1516_ALE/949998‑950178
                                                                                          |                                                                                          
gaCGACTAAAGGAACAGATAAATCGTACCGAAAATCCTACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAt                                                                                            <  2:3920689/90‑1 (MQ=255)
aaCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAt                                                                                            <  1:3628318/91‑1 (MQ=255)
 aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAtt                                                                                           <  2:2362940/91‑1 (MQ=255)
  cGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTa                                                                                          <  2:291203/91‑1 (MQ=255)
  cGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTa                                                                                          >  1:2751727/1‑91 (MQ=255)
    aCTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTaaa                                                                                        >  1:2096158/1‑91 (MQ=255)
      tAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtt                                                                                      <  2:356549/91‑1 (MQ=255)
       aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATCAAAGAGATTAGTCTCTATTAAAttt                                                                                     >  2:3538631/1‑91 (MQ=255)
         aGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTtta                                                                                   <  1:2272006/91‑1 (MQ=255)
             aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATTAAATTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                               <  2:1334939/91‑1 (MQ=37)
             aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAg                                                                               <  1:2174147/91‑1 (MQ=255)
             aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                               <  2:940310/91‑1 (MQ=38)
                 aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                           >  1:2286233/1‑91 (MQ=255)
                 aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                           <  2:2096859/91‑1 (MQ=37)
                  tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACt                                                                          <  1:3092630/91‑1 (MQ=37)
                  tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACt                                                                          <  1:186193/91‑1 (MQ=37)
                   aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTa                                                                         <  2:3650399/91‑1 (MQ=37)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        >  1:1832377/1‑91 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        >  2:1843314/1‑91 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        >  1:3269417/1‑91 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        <  2:2297346/91‑1 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        <  2:606262/91‑1 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        <  1:1030183/91‑1 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                        <  1:1637191/91‑1 (MQ=255)
                    aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTaa                                                                        <  2:3604519/91‑1 (MQ=25)
                     aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                       <  2:3061701/91‑1 (MQ=255)
                     aTCGTACCGAAAATCATACAAAGGGTCTGAAACGAAAGTTTCTTAGACTATCAACGAGATTAGTATCTATTACATTTTATTAGATACTAAt                                                                       >  2:2855690/1‑91 (MQ=255)
                      tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                      <  1:1356330/91‑1 (MQ=255)
                      acGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtc                                                                      <  2:3706625/90‑1 (MQ=25)
                        gTACCGAAAATCCTGCCAAGGTTCTGAATTGAAAGTTTCTTAGACTATAAAAGCGAGTAGTTTCTAATAAAATTTAATAGATACTAAtctc                                                                    <  2:2176738/91‑1 (MQ=16)
                          aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                  >  2:686534/1‑91 (MQ=255)
                          aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATCTCtt                                                                  <  2:2957245/91‑1 (MQ=25)
                           ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                 <  2:2067928/91‑1 (MQ=255)
                            cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                <  1:2702623/91‑1 (MQ=255)
                            cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtgtt                                                                <  2:2032164/91‑1 (MQ=255)
                             gAAAATCATACAAAGGGTCTGTAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                               <  2:1666030/91‑1 (MQ=255)
                              aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTg                                                              >  2:1596935/1‑91 (MQ=255)
                              aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCGTTTTg                                                              >  2:3638572/1‑91 (MQ=255)
                                aaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTc                                                            <  2:1192539/91‑1 (MQ=255)
                                  tCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTa                                                          <  2:1659339/91‑1 (MQ=255)
                                  tCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTACAAAAGAGACTAGTATCTATTATATTTTATTAGATACTAATCTCTTTTTGTCTa                                                          >  1:3152275/1‑91 (MQ=38)
                                  tCATACAAAGGGTCTGAAATGAAAGTTTCATAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTGTTTGTCTa                                                          >  2:3441607/1‑91 (MQ=255)
                                    aTACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACg                                                        <  1:1424401/91‑1 (MQ=255)
                                      aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                      <  2:3041516/91‑1 (MQ=255)
                                      aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                      >  1:1823298/1‑91 (MQ=255)
                                      aCAAAGGGTCTGAAATGAAAGTTTCTGAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                      >  2:2251745/1‑91 (MQ=255)
                                       cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                     <  2:265692/91‑1 (MQ=255)
                                         aaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                   <  1:249575/91‑1 (MQ=255)
                                          aGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACg                                                  >  1:1748127/1‑91 (MQ=255)
                                             gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                               >  1:2392177/1‑91 (MQ=255)
                                             gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                               >  1:3405755/1‑91 (MQ=255)
                                              tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                              <  1:3240945/91‑1 (MQ=255)
                                              tCTGAAACGAAAGTTTCTTAGACTATAAAAGAGACTACCACCTATCAAACTTCACTACACACTAATCTCTTCCTGTCTACGATAACGTAat                                              >  2:3944264/1‑91 (MQ=255)
                                               cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                             <  1:3768654/91‑1 (MQ=255)
                                               cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                             <  1:1171668/91‑1 (MQ=255)
                                               cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                             <  2:1401984/91‑1 (MQ=255)
                                                tgaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAatat                                            >  1:318951/1‑91 (MQ=255)
                                                    atgaaaGTTGCTTAGACTATAAAAGATATTAGTATCTATTGAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgatt                                        <  2:1714919/91‑1 (MQ=255)
                                                       aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTTTATTTTAATTTAATAtgattgat                                     >  2:748015/1‑91 (MQ=25)
                                                       aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                     <  1:763839/91‑1 (MQ=255)
                                                         aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                   >  1:153044/1‑91 (MQ=255)
                                                           tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                 >  1:2250980/1‑91 (MQ=255)
                                                            ttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAt                                >  1:2556781/1‑91 (MQ=255)
                                                             tCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAtt                               >  1:1666836/1‑91 (MQ=255)
                                                              cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                              <  1:1465616/91‑1 (MQ=255)
                                                              cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                              >  2:2565433/1‑91 (MQ=255)
                                                                 aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTTTCTACGATAACGTAATATGATTGATTCTATTTaca                           >  2:2897658/1‑91 (MQ=255)
                                                                 aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                           <  1:3730438/91‑1 (MQ=255)
                                                                  gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                          <  2:3400522/91‑1 (MQ=255)
                                                                  gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                          <  1:3308768/91‑1 (MQ=255)
                                                                   aCTATAAAAGAGATTAGTATCTATTAAATTTTTTTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                         <  1:3376108/91‑1 (MQ=255)
                                                                   aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                         <  1:197140/91‑1 (MQ=255)
                                                                    cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                        >  1:105146/1‑91 (MQ=255)
                                                                    cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                        <  1:2721650/91‑1 (MQ=255)
                                                                        aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                    >  1:2961785/1‑91 (MQ=255)
                                                                        aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                    <  2:3130369/91‑1 (MQ=255)
                                                                           agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                 <  2:272436/91‑1 (MQ=255)
                                                                           agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                 <  2:1530098/91‑1 (MQ=255)
                                                                           agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                 >  1:1331168/1‑91 (MQ=255)
                                                                           agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                 <  2:2956574/91‑1 (MQ=255)
                                                                             agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt               <  1:961390/91‑1 (MQ=255)
                                                                              gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt              >  1:1852739/1‑91 (MQ=255)
                                                                               aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt             >  1:559210/1‑91 (MQ=25)
                                                                                ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa            >  1:143816/1‑91 (MQ=255)
                                                                                ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa            >  1:550178/1‑91 (MQ=255)
                                                                                 tAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTaa           >  2:1744822/1‑91 (MQ=255)
                                                                                  aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg          >  1:1259301/1‑91 (MQ=255)
                                                                                  aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg          >  1:2513135/1‑91 (MQ=255)
                                                                                  aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg          >  2:3742257/1‑91 (MQ=255)
                                                                                  aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg          >  1:2777147/1‑91 (MQ=37)
                                                                                  aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg          >  1:848635/1‑91 (MQ=37)
                                                                                   gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg         <  2:3260708/91‑1 (MQ=255)
                                                                                   gTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg         >  1:3593820/1‑91 (MQ=255)
                                                                                    tATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGt        >  1:1369686/1‑91 (MQ=37)
                                                                                     aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg       >  2:35870/1‑91 (MQ=255)
                                                                                     aTCTAGTAATTTTTAGTAGATACTGATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg       <  1:2339595/91‑1 (MQ=255)
                                                                                      tCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa      >  1:2998895/1‑91 (MQ=37)
                                                                                      ccTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa      >  2:3128212/2‑91 (MQ=255)
                                                                                        tATTAAATTTTATTAGATACTAATCTCTTTTTTTTTATTTTTACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa    >  2:2404500/1‑91 (MQ=255)
                                                                                          ttAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata  <  1:3675858/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATA  >  USA300TCH1516_ALE/949998‑950178

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: