Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 803,539 A→G 6.5% intergenic (‑239/+31) USA300HOU_RS03920 ← / ← dtpT Putative lipid kinase/Di‑/tripeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE803,5390AG 6.5% 1432.8 / 59.7 423intergenic (‑239/+31)USA300HOU_RS03920/dtpTPutative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (186/209);  new base G (20/8);  total (206/217)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 433 total aligned reads displayed.

TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629
                                                                                          |                                                                                          
tCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGa                                                                                            >  1:1143031/1‑91 (MQ=255)
  tataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACt                                                                                          >  1:4917284/1‑91 (MQ=255)
   ataATGAGCATAGTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACtt                                                                                         <  2:4628718/91‑1 (MQ=255)
    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTa                                                                                        <  1:4317220/91‑1 (MQ=255)
    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                                        <  2:4220468/91‑1 (MQ=255)
     aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                                       <  1:5890008/91‑1 (MQ=255)
     aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                                       <  2:4571288/91‑1 (MQ=255)
      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACt                                                                                      >  2:6825247/1‑91 (MQ=255)
       tGAGCATATTAAGTTTGTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                                     <  2:3216521/91‑1 (MQ=255)
          gCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttt                                                                                  <  1:6783684/91‑1 (MQ=255)
           cacatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 <  2:6186657/88‑1 (MQ=255)
           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 >  2:6590087/1‑91 (MQ=255)
           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 <  1:3952039/91‑1 (MQ=255)
           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 >  2:1906494/1‑91 (MQ=255)
            atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttttt                                                                                <  1:4625179/91‑1 (MQ=255)
              atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTg                                                                              >  2:4102432/1‑91 (MQ=255)
                tAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                            <  2:2037992/91‑1 (MQ=255)
                tAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                            <  1:3789007/91‑1 (MQ=255)
                tAAGTTTTTATTTATACTAATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                            <  2:1370650/91‑1 (MQ=255)
                 aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAt                                                                           <  2:5422251/91‑1 (MQ=255)
                 aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAt                                                                           >  1:5311192/1‑91 (MQ=255)
                  aGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAtt                                                                          >  2:5253301/1‑91 (MQ=255)
                    tttttATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAg                                                                        >  1:1323584/1‑91 (MQ=255)
                      tttatttaTACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGc                                                                      >  2:3906829/1‑91 (MQ=255)
                       ttatttatACTTATATCTTGAATTTAATTAATAGAAACCTACAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGcc                                                                     <  2:6977816/91‑1 (MQ=255)
                        tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCca                                                                    <  1:687746/91‑1 (MQ=255)
                        tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCca                                                                    >  1:4783630/1‑91 (MQ=255)
                        tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCca                                                                    >  2:5654775/1‑91 (MQ=255)
                          tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                                  <  2:3415006/91‑1 (MQ=255)
                           ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                                 >  1:2789578/1‑91 (MQ=255)
                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaacaaa                                                               >  2:499722/1‑90 (MQ=255)
                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  1:6580223/1‑91 (MQ=255)
                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  1:2624244/1‑91 (MQ=255)
                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  1:506141/1‑91 (MQ=255)
                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  2:3332115/1‑91 (MQ=255)
                               aCTGATATCTTGCATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAatat                                                             <  2:4476388/91‑1 (MQ=255)
                               aCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAatat                                                             >  1:1021347/1‑91 (MQ=255)
                                cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATAtt                                                            <  1:5577407/91‑1 (MQ=255)
                                cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATAtt                                                            <  2:3806530/91‑1 (MQ=255)
                                cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATAtt                                                            >  2:2864418/1‑91 (MQ=255)
                                   atatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAAc                                                         <  2:3561883/91‑1 (MQ=255)
                                    tatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACg                                                        >  2:6793306/1‑91 (MQ=255)
                                      tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGCTTTTTGAATTAGGCCAACAATATTAACGta                                                      >  1:5397833/1‑91 (MQ=255)
                                       cTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtat                                                     <  2:2545199/91‑1 (MQ=255)
                                       cTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtat                                                     <  2:1862473/91‑1 (MQ=255)
                                       cTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtat                                                     <  2:1145395/91‑1 (MQ=255)
                                        ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                                    >  1:5962355/1‑91 (MQ=255)
                                          gAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATAcc                                                  >  2:6236180/1‑91 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:2872915/1‑91 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                <  1:253681/91‑1 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  2:1239497/1‑91 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                <  1:824925/91‑1 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:5031464/1‑91 (MQ=255)
                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  2:6443115/1‑91 (MQ=255)
                                             tttAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCttt                                               <  2:730964/91‑1 (MQ=255)
                                             tttAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCttt                                               <  1:4533226/91‑1 (MQ=255)
                                              ttaattaaTAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTc                                              <  1:1509092/91‑1 (MQ=255)
                                              ttaattaaTAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTc                                              <  2:2746200/91‑1 (MQ=255)
                                                aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                            <  1:1828278/91‑1 (MQ=255)
                                                 attaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTCTTGAATTAGGCCAACAATATTAACGTATACCTTTCATc                                           >  1:5778684/1‑91 (MQ=255)
                                                  ttaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCg                                          >  1:1439360/1‑91 (MQ=255)
                                                     atAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATcaaag                                       <  1:62633/91‑5 (MQ=255)
                                                      tAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCttt                                      <  2:2528443/91‑1 (MQ=255)
                                                       aGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtg                                     >  1:191427/1‑91 (MQ=255)
                                                          aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                                  >  1:5172518/1‑91 (MQ=255)
                                                          aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                                  >  1:4186417/1‑91 (MQ=255)
                                                           aCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatga                                 >  2:5810180/1‑91 (MQ=255)
                                                           aCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatga                                 >  2:1538717/1‑91 (MQ=255)
                                                            ccTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTgatgat                                <  1:2148543/91‑1 (MQ=255)
                                                            ccTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTgatgat                                <  1:5689704/91‑1 (MQ=255)
                                                            ccTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTgatgat                                <  1:4102432/91‑1 (MQ=255)
                                                             cTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGAtt                               >  1:3177304/1‑91 (MQ=255)
                                                             cTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGAtt                               <  1:2061466/91‑1 (MQ=255)
                                                              tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  1:1221292/1‑91 (MQ=255)
                                                               ataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTCTGATGATTaa                             >  2:5540889/1‑91 (MQ=255)
                                                                taAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGCATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaaa                            >  2:4715016/1‑91 (MQ=255)
                                                                 aaaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                           >  1:2927709/1‑91 (MQ=255)
                                                                 aaaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                           <  2:590737/91‑1 (MQ=255)
                                                                 aaaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                           <  1:6362332/91‑1 (MQ=255)
                                                                  aaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg                          <  2:5995033/91‑1 (MQ=255)
                                                                  aaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg                          <  2:1668418/91‑1 (MQ=255)
                                                                    aaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTg                        <  1:239115/91‑1 (MQ=255)
                                                                       aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgaga                     >  1:3842155/1‑91 (MQ=255)
                                                                         aGTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAtt                   >  2:6888454/1‑91 (MQ=255)
                                                                          gTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  >  1:736999/1‑91 (MQ=255)
                                                                            aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                >  2:664056/1‑91 (MQ=255)
                                                                             agtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTa               >  2:4368947/4‑91 (MQ=255)
                                                                             aGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTa               >  1:3551924/1‑91 (MQ=255)
                                                                              gtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa              >  2:235537/3‑91 (MQ=255)
                                                                               tcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt             >  2:3068032/2‑91 (MQ=255)
                                                                                cATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATa            <  1:613268/91‑1 (MQ=255)
                                                                                 aTTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           <  1:2235748/91‑1 (MQ=255)
                                                                                  tttAAATGGCTTACTGTTTTTTGTTTTATGTCCACACTATTAACGTATACCTTTTTTCGCTTTGATCATTATAGGTGAGATTGCTTCTACa          >  2:3715207/1‑91 (MQ=255)
                                                                                  tttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          <  1:1071081/91‑1 (MQ=255)
                                                                                  tttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:1745576/1‑91 (MQ=255)
                                                                                   ttAAATGGCTTACTGTTTTTTTATTTTTGTTAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         >  2:1703031/1‑91 (MQ=37)
                                                                                    tAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAt        >  1:3395278/1‑91 (MQ=255)
                                                                                      aaTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttg      >  2:750963/1‑91 (MQ=255)
                                                                                       aTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttgt     >  1:5992706/1‑91 (MQ=255)
                                                                                          aCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATtgttgt  >  1:987588/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629

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