Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,088 T→A 14.4% 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,0880TA14.4% 677.0 / 92.0 236intergenic (+82/‑280)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (95/107);  new base A (20/14);  total (115/121)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

Only 100 of 245 total aligned reads displayed.

AAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACA  >  USA300TCH1516_ALE/949997‑950176
                                                                                           |                                                                                        
aaaCGACTAAAGGAACAGATAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAt                                                                                          <  1:2412035/91‑1 (MQ=255)
  aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAtt                                                                                         <  2:838021/91‑1 (MQ=255)
  aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAtt                                                                                         <  2:3069720/91‑1 (MQ=255)
  aCGACCAAAGGAACAGATAAATCGTTCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAtt                                                                                         <  1:3616381/91‑1 (MQ=255)
    gACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCtaataa                                                                                       <  1:3253270/91‑1 (MQ=255)
     aCTAAAGGAGCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTaaa                                                                                      <  1:1307843/91‑1 (MQ=255)
     aCTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGCTTACTCTCTATTaaa                                                                                      >  1:3519687/1‑91 (MQ=255)
     aCTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTaaa                                                                                      >  1:2441250/1‑91 (MQ=255)
     aCTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTaaa                                                                                      >  2:682438/1‑91 (MQ=255)
      cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                     <  1:832818/91‑1 (MQ=255)
      cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                     <  1:3189086/91‑1 (MQ=255)
      cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                     <  2:2299998/91‑1 (MQ=255)
       tAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtt                                                                                    >  1:100958/1‑91 (MQ=255)
        aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                   >  1:3324885/1‑91 (MQ=255)
         aaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtttt                                                                                  <  1:447318/91‑1 (MQ=255)
         aaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtttt                                                                                  >  1:1206032/1‑91 (MQ=255)
         aaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtttt                                                                                  <  2:2116600/91‑1 (MQ=255)
          aGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTtta                                                                                 >  1:3657035/1‑91 (MQ=255)
               cAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGa                                                                            <  2:1367145/91‑1 (MQ=255)
                aGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGAt                                                                           >  2:1921939/1‑91 (MQ=255)
                  aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                         >  1:1695555/1‑91 (MQ=255)
                  aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                         >  2:1888609/1‑91 (MQ=255)
                   tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                        <  2:3752662/91‑1 (MQ=255)
                   tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                        <  1:2477506/91‑1 (MQ=255)
                    aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTa                                                                       <  2:1667039/91‑1 (MQ=37)
                     aaacGTACAGAAAATCATACAAATGATCTGAAATGAAAGTTTCTTAAATTATAATAGAGATTAGTATTTATTAAATTTTATTAGATACTaa                                                                      <  2:1215804/88‑1 (MQ=255)
                     aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                      <  2:1958711/91‑1 (MQ=255)
                     aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                      >  1:3858939/1‑91 (MQ=255)
                     aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                      >  2:2737981/1‑91 (MQ=255)
                      aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                     >  2:1778084/1‑91 (MQ=255)
                       tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                    <  1:1452566/91‑1 (MQ=255)
                       tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                    <  1:3447213/91‑1 (MQ=255)
                        cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGCCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTTAtct                                                                   <  1:3070212/91‑1 (MQ=255)
                        cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtct                                                                   <  1:1992915/91‑1 (MQ=255)
                        cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtct                                                                   <  2:1415998/91‑1 (MQ=25)
                          tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                 <  1:1070822/91‑1 (MQ=255)
                          tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                 <  1:2995360/91‑1 (MQ=255)
                           aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                >  2:3276213/1‑91 (MQ=255)
                             cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                              >  1:2746067/1‑91 (MQ=255)
                              gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                             <  1:255874/91‑1 (MQ=255)
                               aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTg                                                            >  2:2487160/1‑91 (MQ=255)
                                      tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                     >  1:3430889/1‑91 (MQ=255)
                                      tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                     >  1:3756644/1‑91 (MQ=255)
                                        cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                   >  2:3473744/1‑91 (MQ=255)
                                        cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                   >  1:903049/1‑91 (MQ=255)
                                        cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                   >  2:935531/1‑91 (MQ=255)
                                         aaaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATaa                                                  >  1:1909151/1‑91 (MQ=255)
                                            gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                               <  1:819527/91‑1 (MQ=255)
                                            gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                               >  2:3525067/1‑91 (MQ=255)
                                            gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTACTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                               >  2:3194610/1‑91 (MQ=255)
                                              gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                             <  2:2444637/91‑1 (MQ=255)
                                              gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                             >  1:3884510/1‑91 (MQ=255)
                                               tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                            >  2:1403897/1‑91 (MQ=255)
                                                cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                           >  2:1465851/1‑91 (MQ=255)
                                                cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                           <  1:3337268/91‑1 (MQ=255)
                                                 tgaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAatat                                          >  2:2783870/1‑91 (MQ=255)
                                                  gaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtg                                         >  2:3391594/1‑91 (MQ=255)
                                                    aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                       >  1:2373669/1‑91 (MQ=255)
                                                     atgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgatt                                      >  1:689678/1‑91 (MQ=255)
                                                      tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                     <  1:1831514/91‑1 (MQ=255)
                                                       gaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattga                                    <  1:377702/91‑1 (MQ=255)
                                                        aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                   >  1:582449/1‑91 (MQ=255)
                                                        aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                   >  2:3861848/1‑91 (MQ=255)
                                                        aaaGTTTCTTAGACTATAAAAGAGATAAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                   >  2:2211787/1‑91 (MQ=255)
                                                         aaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATgattgatt                                  <  1:2649181/91‑1 (MQ=255)
                                                          aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                 >  2:1144353/1‑91 (MQ=255)
                                                           gTTTCTTAGACTATAGAATAGATTAGTATCTATAAAATTTTATTAGCTACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCt                                <  1:1108127/91‑1 (MQ=38)
                                                            tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                               <  1:2210/91‑1 (MQ=255)
                                                               cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                            >  2:378626/1‑91 (MQ=255)
                                                               cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                            >  2:569462/1‑91 (MQ=255)
                                                                ttAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTa                           >  2:1730076/1‑91 (MQ=255)
                                                                 tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                          >  1:2810474/1‑91 (MQ=255)
                                                                  aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                         <  2:587527/91‑1 (MQ=255)
                                                                   gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                        <  2:903049/91‑1 (MQ=255)
                                                                    aCTATAAAAGATATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                       <  1:3523625/91‑1 (MQ=255)
                                                                     cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                      <  2:2388503/91‑1 (MQ=255)
                                                                      tataAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTa                     <  2:3868633/91‑1 (MQ=255)
                                                                        taAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                   <  2:2059162/91‑1 (MQ=255)
                                                                        taAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                   <  2:3649684/91‑1 (MQ=255)
                                                                        aaaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                   <  1:2559084/90‑1 (MQ=255)
                                                                          aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                 <  1:1723956/91‑1 (MQ=255)
                                                                           aaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAAt                >  2:3580939/1‑91 (MQ=255)
                                                                            agagATTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg               >  1:2363663/1‑91 (MQ=25)
                                                                              agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt             <  2:2444233/91‑1 (MQ=255)
                                                                              agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt             <  2:2296629/91‑1 (MQ=255)
                                                                              agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt             <  2:1518035/91‑1 (MQ=255)
                                                                              agaTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt             >  1:504123/1‑91 (MQ=25)
                                                                               gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt            >  2:2989785/1‑91 (MQ=255)
                                                                               gaTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt            >  2:3063439/1‑91 (MQ=25)
                                                                                aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt           >  2:259194/1‑91 (MQ=25)
                                                                                aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt           >  2:1301229/1‑91 (MQ=25)
                                                                                 ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa          <  1:322940/91‑1 (MQ=255)
                                                                                  tAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTaa         >  1:954078/1‑91 (MQ=255)
                                                                                   aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg        <  2:1046604/91‑1 (MQ=255)
                                                                                   aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTTTATTTTTACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg        >  2:282512/1‑91 (MQ=21)
                                                                                   aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg        >  1:2261645/1‑91 (MQ=37)
                                                                                      aTCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg     >  1:499666/1‑91 (MQ=37)
                                                                                       tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa    <  1:1412663/91‑1 (MQ=255)
                                                                                       tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa    >  1:2025186/1‑91 (MQ=255)
                                                                                         tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa  >  2:271810/1‑91 (MQ=255)
                                                                                           |                                                                                        
AAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACA  >  USA300TCH1516_ALE/949997‑950176

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: