BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,362,261	A→T	8.5%	intergenic (‑193/+25)	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,261	0	A	T	8.5%	517.4 / 31.3	152	intergenic (‑193/+25)	USA300TCH1516_01271/lysP_1	hypothetical protein/Lysine‑specific permease
Reads supporting (aligned to +/- strand): ref base A (71/68); new base T (6/7); total (77/75)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 155 total aligned reads displayed.

ACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTT  >  USA300TCH1516_ALE/1362173‑1362351

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

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

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

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

  tAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCa                                                                                        >  1:1514240/1‑91 (MQ=255)

  tAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCa                                                                                        <  1:1482061/91‑1 (MQ=255)

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

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

   aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                                       >  1:3461631/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:698024/91‑1 (MQ=25)

            acTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTAt                                                                              <  1:3112457/91‑1 (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:1886963/91‑1 (MQ=255)

                         aaaaGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttattta                                                                 <  1:1152296/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                   aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                       <  1:2670208/91‑1 (MQ=255)

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

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

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

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

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

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

                                            cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                              <  1:1261353/91‑1 (MQ=255)

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

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

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

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

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

                                                  tACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGt                                        >  1:1603081/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                          tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                                >  1:4107269/1‑91 (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:2555431/1‑91 (MQ=255)

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

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

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

                                                                      aaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAg                    >  1:3802361/1‑91 (MQ=255)

                                                                       aTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc                   >  1:1367147/1‑91 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                   aTCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTAACAAAtttt       >  1:3724118/1‑91 (MQ=255)

                                                                                     cTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgt     <  1:3375850/91‑1 (MQ=255)

                                                                                      tAAAGCATCTGCCGTATGAGAATTTATTCATTTTTGATTGTCGTAATCGTGTGGTTTTAAATAAATTGTTTCTAGCTTTACAAAttttgtt    >  1:2003101/1‑91 (MQ=255)

                                                                                        aaGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgtttt  >  1:3676276/1‑91 (MQ=255)

                                                                                        aaGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgtttt  >  1:2308104/1‑91 (MQ=255)

ACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTT  >  USA300TCH1516_ALE/1362173‑1362351

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