BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	293,940	0	C	T	21.4%	334.1 / 84.3	159	intergenic (+57/‑471)	gatD/tarI'	Galactitol‑1‑phosphate 5‑dehydrogenase/Ribitol‑5‑phosphate cytidylyltransferase 2
Reads supporting (aligned to +/- strand): ref base C (61/64); new base T (19/15); total (80/79)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 171 total aligned reads displayed.

ACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAA  >  USA300TCH1516_ALE/293850‑294029

aCGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGCGGACAAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGc                                                                                           >  1:2458502/1‑91 (MQ=255)

aCGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGc                                                                                           <  2:3031828/91‑1 (MQ=255)

 cGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGcc                                                                                          >  2:2960714/1‑91 (MQ=255)

 cGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGcc                                                                                          >  1:2079044/1‑91 (MQ=255)

   aTTTCATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGt                                                                                        <  1:2575936/91‑1 (MQ=255)

    tttGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATgtcgtc                                                                                       >  1:350787/1‑91 (MQ=255)

    tttGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATgtcgtc                                                                                       >  1:2425442/1‑91 (MQ=255)

    tttGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTc                                                                                       >  1:485506/1‑91 (MQ=255)

     ttGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCt                                                                                      >  2:1535049/1‑91 (MQ=255)

       gATAAGGTCATGTTTACGATTTATTGGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTgg                                                                                    <  1:348342/91‑1 (MQ=255)

       gATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAGTGCCGTCTgg                                                                                    <  2:945627/91‑1 (MQ=255)

              tCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAg                                                                             <  2:1360076/91‑1 (MQ=255)

              tCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAg                                                                             <  2:1694984/91‑1 (MQ=255)

               cATGTTTACGATTTATGAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATTGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGa                                                                            <  2:1928625/91‑1 (MQ=255)

                aTGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGAt                                                                           <  1:2327375/91‑1 (MQ=255)

                aTGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGAt                                                                           <  1:2177230/91‑1 (MQ=255)

                aTGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGAt                                                                           <  2:3094194/91‑1 (MQ=255)

                    ttACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTa                                                                       <  2:431391/91‑1 (MQ=255)

                     tACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTAt                                                                      >  2:2578944/1‑91 (MQ=255)

                      aCGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATc                                                                     >  2:2096804/1‑91 (MQ=255)

                       cGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCg                                                                    >  2:520206/1‑91 (MQ=255)

                       cGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCg                                                                    >  2:1966300/1‑91 (MQ=255)

                       cGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCg                                                                    >  1:1812170/1‑91 (MQ=255)

                         aTTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGaa                                                                  >  1:836479/1‑91 (MQ=255)

                         aTTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGaa                                                                  >  1:701833/1‑91 (MQ=255)

                           ttattaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAaat                                                                <  1:1265938/91‑1 (MQ=255)

                            tattaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaata                                                               >  1:2122133/1‑91 (MQ=255)

                            tattaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaata                                                               >  2:220008/1‑91 (MQ=255)

                            tattaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaata                                                               <  1:2962300/91‑1 (MQ=255)

                            tattaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaata                                                               >  1:2121896/1‑91 (MQ=255)

                             attaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaataa                                                              <  1:2565263/91‑1 (MQ=255)

                             attaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAaataa                                                              >  2:2640471/1‑91 (MQ=255)

                              ttaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaataat                                                             <  2:1915213/91‑1 (MQ=255)

                              ttaGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAaataat                                                             >  1:878006/1‑91 (MQ=255)

                                 gTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCa                                                          >  2:2433693/1‑91 (MQ=255)

                                 gTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCa                                                          >  1:426606/1‑91 (MQ=255)

                                   aTGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATc                                                        <  2:2169804/91‑1 (MQ=255)

                                     gCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCtt                                                      <  1:2699388/91‑1 (MQ=255)

                                      cACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTc                                                     >  2:396668/1‑91 (MQ=255)

                                       aCCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATACATCTTCg                                                    <  2:1081919/91‑1 (MQ=255)

                                           ttGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAata                                                <  1:356641/91‑1 (MQ=255)

                                               ggACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACg                                            >  2:2822157/1‑91 (MQ=255)

                                                gACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGt                                           >  1:2921937/1‑91 (MQ=255)

                                                gACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGt                                           >  2:814244/1‑91 (MQ=255)

                                                gACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGt                                           >  2:2671361/1‑91 (MQ=255)

                                                gACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGt                                           <  1:2447797/91‑1 (MQ=255)

                                                gACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGt                                           <  1:2447823/91‑1 (MQ=255)

                                                 aCGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGtt                                          >  2:228991/1‑91 (MQ=255)

                                                  cGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTg                                         <  2:2009357/91‑1 (MQ=255)

                                                     aaaCGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATa                                      <  2:2866439/91‑1 (MQ=255)

                                                      aaCGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATaa                                     <  2:1881493/91‑1 (MQ=255)

                                                      aaCGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATaa                                     <  2:11490/91‑1 (MQ=255)

                                                       aCGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATaaa                                    <  1:1575029/91‑1 (MQ=255)

                                                        cGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAAt                                   >  1:2013143/1‑91 (MQ=255)

                                                            ggTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGcc                               >  1:1056632/1‑91 (MQ=255)

                                                            ggTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGcc                               <  2:1931224/91‑1 (MQ=255)

                                                            ggTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGcc                               <  1:895561/91‑1 (MQ=255)

                                                             gTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCg                              >  1:2751638/1‑91 (MQ=255)

                                                              tataGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCgg                             <  2:2758420/91‑1 (MQ=255)

                                                              tataGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCgg                             <  2:2756604/91‑1 (MQ=255)

                                                                taGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGtt                           <  1:560736/91‑1 (MQ=255)

                                                                      taGCTATCAAAGTGCGAATGCCGTCTGGTCTATAGAAACTATCGAAAAAATTCATCTTCGAATATATATTGATAAACAGCCGGTTTACTtg                     >  2:2885335/1‑91 (MQ=255)

                                                                        gCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTtgtg                   >  1:251643/1‑91 (MQ=255)

                                                                        gCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTtgtg                   >  1:892768/1‑91 (MQ=255)

                                                                          tATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTtgtgtg                 <  1:1373728/91‑1 (MQ=255)

                                                                          caTGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTtgtgtg                 <  2:3072381/90‑1 (MQ=255)

                                                                           aTGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGa                >  2:1658247/1‑91 (MQ=255)

                                                                           aTGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGa                <  1:2313263/91‑1 (MQ=255)

                                                                           aTGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGa                >  2:1658254/1‑91 (MQ=255)

                                                                            tGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaa               >  2:484779/1‑91 (MQ=255)

                                                                            tGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaa               >  1:1189950/1‑91 (MQ=255)

                                                                             gAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaaa              <  1:2756284/91‑1 (MQ=255)

                                                                             gAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaaa              <  1:999318/91‑1 (MQ=255)

                                                                             gAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaaa              <  1:3035489/91‑1 (MQ=255)

                                                                             gAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGaaa              <  2:2292670/91‑1 (MQ=255)

                                                                               aaGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAata            <  2:511766/91‑1 (MQ=255)

                                                                                 gTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATg          <  2:1202513/91‑1 (MQ=255)

                                                                                 gTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATg          <  1:1041464/91‑1 (MQ=255)

                                                                                 gTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATg          <  2:1202395/91‑1 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  1:914713/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  2:1797035/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  1:1529712/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  1:719861/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  1:1742124/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  2:607611/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  2:609543/1‑91 (MQ=255)

                                                                                  tGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTCGATAAATAGCCGGTTTACTTGTGTGAAATATGc         >  2:454239/1‑91 (MQ=255)

                                                                                   gCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCt        >  1:187760/1‑91 (MQ=255)

                                                                                   gCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCt        <  1:1073970/91‑1 (MQ=255)

                                                                                    cGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCtt       <  2:1762637/91‑1 (MQ=255)

                                                                                    cGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCtt       >  1:924085/1‑91 (MQ=255)

                                                                                     gAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGACATATGCTtg      >  1:2964241/1‑91 (MQ=255)

                                                                                      aaTGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTtgt     <  2:602164/91‑1 (MQ=255)

                                                                                      aaTGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTtgt     <  2:1217501/91‑1 (MQ=255)

                                                                                       aTGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTtgtg    >  1:2824059/1‑91 (MQ=255)

                                                                                        tGCCGTCTGGTCTATAGATACTATCGAAATGATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGGGTGAAATATGCTTGTGa   <  2:2705053/91‑1 (MQ=255)

                                                                                        tGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGa   >  2:1387867/1‑91 (MQ=255)

                                                                                        tGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGa   >  1:2216996/1‑91 (MQ=255)

                                                                                         gCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGaa  <  2:398681/91‑1 (MQ=255)

                                                                                         gCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGaa  >  2:1417771/1‑91 (MQ=255)

ACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAA  >  USA300TCH1516_ALE/293850‑294029

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: `‑`