Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,909,134 T→A 7.6% intergenic (‑200/‑60) tal ← / → USA300HOU_RS09460 Transaldolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,909,1340TA 7.6% 756.7 / 40.8 233intergenic (‑200/‑60)tal/USA300HOU_RS09460Transaldolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (90/125);  new base A (8/10);  total (98/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 238 total aligned reads displayed.

ATATCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTA  >  USA300TCH1516_ALE/1909044‑1909223
                                                                                          |                                                                                         
atatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAACATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACt                                                                                           <  2:3200208/91‑1 (MQ=255)
atatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACt                                                                                           >  1:3821573/1‑91 (MQ=255)
 tatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTTTGTCTACAAGTTAATAAAACTa                                                                                          >  1:1595833/1‑91 (MQ=255)
  atCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTAt                                                                                         >  2:2627119/1‑91 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                      >  2:3828899/1‑91 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                      <  1:389813/91‑1 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                      >  2:1529449/1‑91 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                    <  1:3733544/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                    <  1:902723/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                    <  1:2369665/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                    <  1:903134/91‑1 (MQ=255)
        cTTTTATAATATATAGCTAGTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                   <  1:271345/91‑1 (MQ=255)
        cTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                   <  1:1895307/91‑1 (MQ=255)
          tttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTAt                                                                                 <  1:1765537/91‑1 (MQ=255)
          tttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGTTTTAt                                                                                 <  1:3853556/91‑1 (MQ=255)
           ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                <  2:3697598/91‑1 (MQ=255)
            tataATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGTAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGa                                                                               >  1:3211121/1‑91 (MQ=255)
            tataATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGa                                                                               <  1:2152507/91‑1 (MQ=255)
            tataATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGa                                                                               <  1:121590/91‑1 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAATCTATTGTTTTATGAt                                                                              <  1:3307555/91‑1 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                              <  1:341043/91‑1 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                              <  1:3505642/91‑1 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                              >  1:2480248/1‑91 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                             >  2:1954230/1‑91 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGTTTTAttaac                                                                             <  1:1081748/91‑5 (MQ=255)
               aataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTa                                                                                          <  2:3858324/77‑1 (MQ=255)
               aataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTg                                                                            <  2:2737250/91‑1 (MQ=255)
                atatatAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGt                                                                           <  1:2012487/91‑1 (MQ=255)
                  atataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGttt                                                                         >  1:1156670/1‑91 (MQ=255)
                   tataGCTAATTTTTTTCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGtttt                                                                        >  2:2888660/1‑90 (MQ=255)
                   tataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTa                                                                        <  1:3652729/91‑1 (MQ=255)
                    ataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTAt                                                                       <  2:2908323/91‑1 (MQ=255)
                       gagaaTTTTGTGCATATATGTTTTAAAAGTGACTTTAGCAAACATTGTGGCTACAAGTTAATAAAAGTATGGTTTTATGATGGTTTAGGAt                                                                    <  2:321721/88‑1 (MQ=255)
                         tAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATAc                                                                  <  2:22163/91‑1 (MQ=255)
                          aaTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACg                                                                 <  1:828234/91‑1 (MQ=255)
                          aaTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACg                                                                 <  2:875620/91‑1 (MQ=255)
                           aTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGa                                                                <  1:2780796/91‑1 (MQ=255)
                           aTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGa                                                                <  1:164552/91‑1 (MQ=255)
                            ttttttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGat                                                               <  2:1559761/91‑1 (MQ=255)
                             tttttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                      >  1:2659834/1‑67 (MQ=255)
                                 tGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTa                                                          >  2:2046173/1‑91 (MQ=255)
                                 tGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTa                                                          <  1:2110550/91‑1 (MQ=255)
                                  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaa                                                         >  2:2391103/1‑91 (MQ=255)
                                  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaa                                                         >  1:3697329/1‑91 (MQ=255)
                                    cTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAAt                                                       <  2:1562027/91‑1 (MQ=255)
                                     tataTTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATg                                                      <  1:2731816/91‑1 (MQ=255)
                                        atTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATa                                                   >  2:3403943/1‑91 (MQ=255)
                                          tttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAat                                                 <  2:1492550/91‑1 (MQ=255)
                                           ttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatc                                                >  1:1060881/1‑91 (MQ=255)
                                              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                             >  1:2539579/1‑91 (MQ=255)
                                              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                             >  1:2167285/1‑91 (MQ=255)
                                              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                             <  1:1176175/91‑1 (MQ=255)
                                                aaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGCTTGTTTATGATACGATATGTAAATGATAATCATCCa                                           >  1:3347424/1‑91 (MQ=255)
                                                aaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCa                                           <  2:3733407/91‑1 (MQ=255)
                                                  aaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCagg                                         <  1:2308554/91‑1 (MQ=255)
                                                  aaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCagg                                         <  1:636019/91‑1 (MQ=255)
                                                  aaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCagg                                         >  1:496761/1‑91 (MQ=255)
                                                    tttcTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                       >  2:2624099/4‑91 (MQ=255)
                                                       cTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                     >  2:252792/1‑90 (MQ=255)
                                                       cTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTc                                    <  2:133276/91‑1 (MQ=255)
                                                       cTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTc                                    <  2:2106/91‑1 (MQ=255)
                                                         ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGCAGGTCTa                                  >  1:1937384/1‑91 (MQ=255)
                                                          tAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAg                                 >  1:807828/1‑91 (MQ=255)
                                                          tAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAg                                 >  2:2697949/1‑91 (MQ=255)
                                                           aGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGc                                >  1:2228001/1‑91 (MQ=255)
                                                           aGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGc                                <  2:583333/91‑1 (MQ=255)
                                                           aGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGc                                >  1:582976/1‑91 (MQ=255)
                                                            gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                               <  2:3576441/91‑1 (MQ=255)
                                                            gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                               <  1:3038476/91‑1 (MQ=255)
                                                            gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                               >  2:1122700/1‑91 (MQ=255)
                                                                 cATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtct                          <  1:1298794/91‑1 (MQ=255)
                                                                 cATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtct                          >  2:651185/1‑91 (MQ=255)
                                                                    tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                       <  2:621038/91‑1 (MQ=255)
                                                                    tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                       <  1:2183788/91‑1 (MQ=255)
                                                                    tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                       <  1:948250/91‑1 (MQ=255)
                                                                    tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                       >  2:771054/1‑91 (MQ=255)
                                                                      tgtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTc                     >  2:1670597/1‑91 (MQ=255)
                                                                      tgtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTc                     <  2:1010282/91‑1 (MQ=255)
                                                                       gtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCa                    >  2:1075243/1‑91 (MQ=255)
                                                                         cTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaa                  <  1:2801345/91‑1 (MQ=255)
                                                                          tACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaa                 <  1:1036777/91‑1 (MQ=255)
                                                                          tACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaa                 <  1:670569/91‑1 (MQ=255)
                                                                           aCAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaa                >  2:3805620/1‑91 (MQ=255)
                                                                            cAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaaa               <  2:2258642/91‑1 (MQ=255)
                                                                             aaGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaaaa              <  2:1880500/91‑1 (MQ=255)
                                                                               gTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGTGGTCTAGCATGTCTCGTTCAAAAAAATa            <  1:1674232/91‑1 (MQ=255)
                                                                               gTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATa            >  2:1636818/1‑91 (MQ=255)
                                                                                 taataaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACt          >  2:2668872/1‑91 (MQ=255)
                                                                                  cataaAACAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtt         >  1:877257/2‑91 (MQ=255)
                                                                                  cataaAACAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtt         >  2:1142655/2‑91 (MQ=255)
                                                                                  cataaAACAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtt         >  1:3080608/2‑91 (MQ=255)
                                                                                  aataaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGATAAAAAAAATACtt         <  2:2972910/91‑1 (MQ=255)
                                                                                   ataaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACttt        <  1:2911235/91‑1 (MQ=255)
                                                                                    taaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt       >  2:1462356/1‑91 (MQ=255)
                                                                                      aaaCAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTCTTAc     >  2:685349/1‑91 (MQ=255)
                                                                                       aaCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACt    <  1:517678/91‑1 (MQ=255)
                                                                                        aCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACtt   >  2:3050660/1‑91 (MQ=255)
                                                                                        aCTATTGTATTAAAATTGTTTATTATACTATATGTAAATAATAATCATCCAGGTGGTCTATCAAATAAAATTAAAAAAAATACTTTTACtt   <  2:871736/91‑1 (MQ=255)
                                                                                         cTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTa  <  2:3320332/91‑1 (MQ=255)
                                                                                         cTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCAGCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTa  >  1:3086883/1‑91 (MQ=255)
                                                                                          |                                                                                         
ATATCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTA  >  USA300TCH1516_ALE/1909044‑1909223

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: