BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,362,255	A→T	9.5%	intergenic (‑187/+31)	USA300TCH1516_01271 ← / ← lysP_1	hypothetical protein/Lysine‑specific permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,362,255	0	A	T	9.5%	492.0 / 34.7	149	intergenic (‑187/+31)	USA300TCH1516_01271/lysP_1	hypothetical protein/Lysine‑specific permease
Reads supporting (aligned to +/- strand): ref base A (65/68); new base T (6/8); total (73/76)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 154 total aligned reads displayed.

AACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTT  >  USA300TCH1516_ALE/1362165‑1362345

aaCAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATa                                                                                            >  1:1387146/1‑91 (MQ=255)

aaCAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCAGTAAAAAGCGGCAAATGCTATAGATa                                                                                            >  1:3852792/1‑91 (MQ=255)

  cAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAAt                                                                                          <  1:2408778/91‑1 (MQ=255)

    gAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCt                                                                                        <  1:3122411/91‑1 (MQ=255)

        aCTAGTTAATACACTTTATACTAGAAAAAGATATTAGAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAg                                                                                    >  1:3062155/1‑91 (MQ=255)

        aCTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTatag                                                                                    <  1:1145719/91‑4 (MQ=38)

        aCTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTGTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAg                                                                                    <  1:3585469/91‑1 (MQ=255)

         cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                                   <  1:3789183/91‑1 (MQ=255)

         cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                                   >  1:732616/1‑91 (MQ=255)

         cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                                   >  1:677483/1‑91 (MQ=255)

           aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCAt                                                                                 <  1:3170388/91‑1 (MQ=37)

           aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAATGCAt                                                                                 >  1:3461809/1‑91 (MQ=255)

           aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                                 <  1:3911212/91‑1 (MQ=255)

              tAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTg                                                                              >  1:1389657/1‑91 (MQ=255)

                aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                            <  1:698011/91‑1 (MQ=25)

                aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                            <  1:697936/91‑1 (MQ=25)

                aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                            <  1:1378486/91‑1 (MQ=25)

                aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                            <  1:698024/91‑1 (MQ=25)

                  acacTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGt                                                                          <  1:459197/91‑1 (MQ=255)

                   cacTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTa                                                                         >  1:3252012/1‑91 (MQ=255)

                    acTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTAt                                                                        <  1:3112457/91‑1 (MQ=255)

                     cTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATg                                                                       <  1:1695878/91‑1 (MQ=255)

                         ataCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGaa                                                                   >  1:3467090/1‑91 (MQ=255)

                          taCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAat                                                                  >  1:2476968/1‑91 (MQ=255)

                            cTAGAAAAAGATATTAATGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                <  1:299874/91‑1 (MQ=255)

                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                <  1:1146286/91‑1 (MQ=255)

                              aGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttat                                                              <  1:3450152/91‑1 (MQ=255)

                              aGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttat                                                              <  1:1886963/91‑1 (MQ=255)

                                    aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttattt                                                        >  1:1506238/1‑91 (MQ=255)

                                     gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                       <  1:2451364/91‑1 (MQ=255)

                                     gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                       >  1:3717719/1‑91 (MQ=255)

                                     gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                       >  1:3720295/1‑91 (MQ=255)

                                      atatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAttttt                                                      >  1:3359013/1‑91 (MQ=255)

                                       tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGTAAATGCTATAGATGATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                                     >  1:3855242/1‑91 (MQ=255)

                                       tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                                     <  1:2337361/91‑1 (MQ=255)

                                       tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGAAAATGCTATAGATCATCTAAAGCATCTGCCGTATGAGAATTTATTTACTTTTg                                                     >  1:3857455/1‑91 (MQ=255)

                                          tAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGAtt                                                  <  1:3915879/91‑1 (MQ=255)

                                           aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                 >  1:2880274/1‑91 (MQ=255)

                                           aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                 >  1:1022567/1‑91 (MQ=255)

                                           aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                 >  1:4104132/1‑91 (MQ=255)

                                           aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                 >  1:528871/1‑91 (MQ=255)

                                            aaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCCGCCGTATGAGAATTTATTTATTTTTGATTGt                                                <  1:2970533/91‑1 (MQ=255)

                                               tataTCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCAt                                             <  1:3647788/91‑1 (MQ=255)

                                                 tatCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATaa                                           <  1:1558187/91‑1 (MQ=255)

                                                   tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                         <  1:1552157/91‑1 (MQ=255)

                                                   tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                         >  1:1874018/1‑91 (MQ=255)

                                                    cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                        >  1:314658/1‑91 (MQ=255)

                                                    cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                        >  1:314723/1‑91 (MQ=255)

                                                      gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                      <  1:178789/91‑1 (MQ=255)

                                                      gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                      >  1:623584/1‑91 (MQ=255)

                                                      gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                      <  1:2854467/91‑1 (MQ=255)

                                                       cTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtgt                                     <  1:3246334/91‑1 (MQ=255)

                                                        tttACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGtgtg                                    <  1:3936408/91‑1 (MQ=255)

                                                           acacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtt                                 <  1:3190475/91‑1 (MQ=255)

                                                            cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                                >  1:2470728/1‑91 (MQ=255)

                                                             acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                               <  1:2844424/91‑1 (MQ=255)

                                                              ccATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTa                              >  1:400118/1‑91 (MQ=255)

                                                               cATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaa                             >  1:1088112/1‑91 (MQ=255)

                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            <  1:73690/91‑1 (MQ=255)

                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            <  1:3897220/91‑1 (MQ=255)

                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            <  1:2174503/91‑1 (MQ=255)

                                                                  tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                          >  1:4107269/1‑91 (MQ=255)

                                                                  tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                          <  1:2524636/91‑1 (MQ=255)

                                                                    aaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaa                        >  1:587802/1‑91 (MQ=255)

                                                                     catacGGCAGATGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                       >  1:2687099/5‑91 (MQ=21)

                                                                     aaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                       <  1:924729/91‑1 (MQ=255)

                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      <  1:2803065/91‑1 (MQ=255)

                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      >  1:2024162/1‑91 (MQ=255)

                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      <  1:2803066/91‑1 (MQ=255)

                                                                        gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt                    >  1:1245515/1‑91 (MQ=255)

                                                                        gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt                    >  1:3810867/1‑91 (MQ=255)

                                                                         cGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGtt                   >  1:1363194/1‑91 (MQ=255)

                                                                          ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                  >  1:3421450/1‑91 (MQ=255)

                                                                          ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                  <  1:2460102/91‑1 (MQ=255)

                                                                          ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                  >  1:2555431/1‑91 (MQ=255)

                                                                           gCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTc                 >  1:3889437/1‑91 (MQ=255)

                                                                            cAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCt                >  1:3072076/1‑91 (MQ=255)

                                                                             aaaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa               >  1:1176849/1‑91 (MQ=255)

                                                                               aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTGTTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc             <  1:3288553/91‑1 (MQ=255)

                                                                               aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc             >  1:2594066/1‑91 (MQ=255)

                                                                               aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc             <  1:2181886/91‑1 (MQ=255)

                                                                                  cTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt          <  1:1606926/91‑1 (MQ=255)

                                                                                  cTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt          <  1:1981511/91‑1 (MQ=255)

                                                                                   tttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         >  1:2120107/3‑91 (MQ=38)

                                                                                   tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         >  1:2876222/1‑91 (MQ=255)

                                                                                   tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         >  1:1397982/1‑91 (MQ=255)

                                                                                   tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         >  1:1486388/1‑91 (MQ=255)

                                                                                    ttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        >  1:1731862/2‑91 (MQ=39)

                                                                                    ttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        >  1:3063147/2‑91 (MQ=39)

                                                                                    ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATACTCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        >  1:3388732/1‑91 (MQ=255)

                                                                                    ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        >  1:2558864/1‑91 (MQ=255)

                                                                                    ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        <  1:3859469/91‑1 (MQ=255)

                                                                                    ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc        <  1:3584722/91‑1 (MQ=255)

                                                                                     taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa       <  1:2682773/91‑1 (MQ=255)

                                                                                     taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa       >  1:3563691/1‑91 (MQ=255)

                                                                                      cgATAAGCTAACGCATCTGCCGTATGAGAATGTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa      <  1:1509776/90‑1 (MQ=255)

                                                                                      aGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa      >  1:2565052/1‑91 (MQ=255)

                                                                                      aGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa      <  1:3023419/91‑1 (MQ=255)

                                                                                        aTAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAt    <  1:1303811/91‑1 (MQ=255)

                                                                                         tAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAtt   <  1:3053396/91‑1 (MQ=255)

AACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTT  >  USA300TCH1516_ALE/1362165‑1362345

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