Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 803,539 A→G 17.2% 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,5390AG17.2% 564.5 / 111.3 250intergenic (‑239/+31)USA300HOU_RS03920/dtpTPutative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (106/101);  new base G (30/13);  total (136/114)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 253 total aligned reads displayed.

TATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803451‑803629
                                                                                        |                                                                                          
tataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAAAAAAATATAAAAAAACAGTAAGCCATTTAAATGACt                                                                                          <  2:55164/91‑1 (MQ=255)
 ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                         <  1:416639/91‑1 (MQ=255)
 ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                         <  1:21532/91‑1 (MQ=255)
 ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACtt                                                                                         >  2:2838095/1‑91 (MQ=255)
  taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                                        <  2:2172995/91‑1 (MQ=255)
  taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                                        >  2:1178212/1‑91 (MQ=255)
   aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                                       <  1:1091071/91‑1 (MQ=255)
     tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                                     >  2:2931408/1‑91 (MQ=255)
      gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                                    >  2:1958068/1‑91 (MQ=255)
      gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                                    <  2:1139360/91‑1 (MQ=255)
       aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTACTAGAAACCTATAAAAATCCAGTAAGCCATGTAAATGACTCACTGtt                                                                                   >  2:2439442/1‑91 (MQ=255)
       aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGATACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  1:3358979/91‑1 (MQ=37)
       aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  2:3171025/91‑1 (MQ=255)
       aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTAACTGtt                                                                                   <  1:3022406/91‑1 (MQ=38)
           tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                               <  2:1924861/91‑1 (MQ=255)
           tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                               <  2:239644/91‑1 (MQ=255)
           tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                               >  1:1503271/1‑91 (MQ=255)
             ttAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGa                                                                             >  2:3592008/1‑91 (MQ=255)
              tAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                            <  1:3155505/91‑1 (MQ=255)
              tAAGTTTTTATTTATACTGATATCTTGAATTCAATTAATAGAAACCTATAAAAAAACCGCACGCCATTTAAAGGACTTACTGTTTTTTGaa                                                                            >  2:1243280/1‑91 (MQ=255)
               aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAt                                                                           >  1:260046/1‑91 (MQ=255)
               aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAt                                                                           >  1:3847547/1‑91 (MQ=255)
                 gTTTTTATTTATACTGATATCTTGAGTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTa                                                                         >  1:2003538/1‑91 (MQ=255)
                      tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCca                                                                    >  1:3260438/1‑91 (MQ=255)
                       atttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaa                                                                   >  1:1389768/1‑91 (MQ=255)
                         ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                                 <  2:811478/91‑1 (MQ=255)
                         ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                                 <  1:3208763/91‑1 (MQ=255)
                          tataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaacaa                                                                >  1:635739/1‑91 (MQ=255)
                            taCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAata                                                              <  1:478237/91‑1 (MQ=255)
                            taCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAata                                                              >  2:1157964/1‑91 (MQ=255)
                               tGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTa                                                           >  1:1267205/1‑91 (MQ=255)
                                gATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTaa                                                          >  1:1579539/1‑91 (MQ=255)
                                 atatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAAc                                                         <  1:355649/91‑1 (MQ=255)
                                    tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                                      <  2:1484300/91‑1 (MQ=255)
                                    tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                                      >  1:2218229/1‑91 (MQ=255)
                                      ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                                    >  2:3423943/1‑91 (MQ=255)
                                      ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                                    >  2:3412754/1‑91 (MQ=255)
                                      ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                                    <  1:2518869/91‑1 (MQ=255)
                                       tGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATAc                                                   >  2:400908/1‑91 (MQ=255)
                                         aaTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCt                                                 <  1:2245848/91‑1 (MQ=255)
                                          aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:2094057/1‑91 (MQ=255)
                                           tttAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCttt                                               <  2:3236057/91‑1 (MQ=255)
                                              aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                            >  1:375168/1‑91 (MQ=255)
                                              aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                            >  1:3789529/1‑91 (MQ=255)
                                               attaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATc                                           <  1:3583714/91‑1 (MQ=255)
                                               attaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATc                                           >  1:695518/1‑91 (MQ=255)
                                                 taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                                         <  2:1736153/91‑1 (MQ=255)
                                                  aatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCt                                        >  2:1999173/1‑91 (MQ=255)
                                                      gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                                    <  2:2984807/91‑1 (MQ=255)
                                                      gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                                    <  1:2809710/91‑1 (MQ=255)
                                                        aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTCACGTATACCTTTCATTGCTTtgctg                                  >  2:822130/1‑91 (MQ=255)
                                                         aCCTATAAAAAAACAGTAAGCCATTTAAATGAGTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTCTCATCGCTTtgatga                                 <  1:3314093/91‑1 (MQ=255)
                                                         aCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatga                                 <  1:953035/91‑1 (MQ=255)
                                                            tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  2:3321419/1‑91 (MQ=255)
                                                            tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  2:140080/1‑91 (MQ=255)
                                                            tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  1:2645198/1‑91 (MQ=255)
                                                            tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  2:2614995/1‑91 (MQ=255)
                                                             ataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaa                             >  2:3772097/1‑91 (MQ=255)
                                                             ataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaa                             >  1:1705216/1‑91 (MQ=255)
                                                                 aaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt                         >  2:1042295/1‑91 (MQ=255)
                                                                   aaaCCGTAAGAAATTTACATGACTTACTGTTTTTTGAATTAGGCGAACCATATTAACGTATACATTTCATAGCTTTGATGATTAAAGGTga                       <  1:3084051/91‑1 (MQ=255)
                                                                   aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga                       <  1:882446/91‑1 (MQ=255)
                                                                   aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga                       <  1:1776142/91‑1 (MQ=255)
                                                                   aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga                       <  2:2351113/91‑1 (MQ=255)
                                                                   aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga                       >  2:713848/1‑91 (MQ=255)
                                                                     aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGGCAACAATATTAACGTATACCTTTCATCGCTTTGTTGATTAAAGGTgaga                     >  2:2101982/1‑91 (MQ=255)
                                                                     aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGAATAAAGGTgcga                     >  2:2702267/1‑91 (MQ=255)
                                                                        gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  >  1:1161782/1‑91 (MQ=255)
                                                                        gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  >  2:2305584/1‑91 (MQ=255)
                                                                        gTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  <  2:894369/91‑1 (MQ=255)
                                                                         tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  2:637616/1‑91 (MQ=255)
                                                                         tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  1:1482441/1‑91 (MQ=255)
                                                                         tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  2:3044332/1‑91 (MQ=255)
                                                                         tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  1:3451669/1‑91 (MQ=255)
                                                                         tAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  1:649061/1‑91 (MQ=255)
                                                                          aagtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                >  1:132245/5‑91 (MQ=255)
                                                                          aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                >  1:66571/1‑91 (MQ=255)
                                                                          aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                <  2:654694/91‑1 (MQ=255)
                                                                           agtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGCGCGATTGCTa               >  2:1606553/4‑91 (MQ=38)
                                                                           agtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGCATACCTTTCATCGCTTTGATGATTAAAGGAGAGATTGCTa               >  1:1402876/4‑91 (MQ=38)
                                                                            gCCATTTATTTGACTTACTGTTTTTTGAATTTGTCCAACTCTATTAACGTATACCTTTTGTCGCTTTGGTGATTCGAGGTGAGATTGCTaa              <  2:2413839/91‑1 (MQ=255)
                                                                            gCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa              <  2:3471778/91‑1 (MQ=255)
                                                                            gCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa              >  1:2079482/1‑91 (MQ=255)
                                                                             tcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt             >  2:2263481/2‑91 (MQ=255)
                                                                             ccATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt             <  2:47169/91‑1 (MQ=255)
                                                                              cATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt             >  1:2384854/1‑90 (MQ=255)
                                                                              cATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATa            >  2:2648622/1‑91 (MQ=255)
                                                                               aTTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           <  2:1768045/91‑1 (MQ=255)
                                                                               aTTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           <  1:287036/91‑1 (MQ=255)
                                                                               aTTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           >  1:524777/1‑91 (MQ=255)
                                                                                tttAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATATCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  2:3399590/1‑91 (MQ=255)
                                                                                 ttAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         >  2:438126/1‑91 (MQ=255)
                                                                                 ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:61571/91‑1 (MQ=255)
                                                                                  tAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAt        >  1:2443112/1‑91 (MQ=255)
                                                                                  tAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAt        >  2:1000702/1‑91 (MQ=255)
                                                                                    aaTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttg      >  1:3040775/1‑91 (MQ=255)
                                                                                     aTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttgt     <  1:2350793/91‑1 (MQ=255)
                                                                                       gACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttgttg   >  1:3023384/1‑91 (MQ=255)
                                                                                        aCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATtgttgt  >  1:747344/1‑91 (MQ=255)
                                                                                        aCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAACACAATtgttgt  >  1:3520682/1‑91 (MQ=255)
                                                                                        |                                                                                          
TATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803451‑803629

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