BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,432,237	C→T	5.3%	intergenic (+584/+63)	birA_2 → / ← USA300HOU_RS12525	Bifunctional ligase/repressor BirA/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,432,237	0	C	T	5.3%	378.7 / 7.6	112	intergenic (+584/+63)	birA_2/USA300HOU_RS12525	Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (44/62); new base T (3/3); total (47/65)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 114 total aligned reads displayed.

TTTAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATAT  >  USA300TCH1516_ALE/2432147‑2432326

tttAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATc                                                                                           >  1:3473343/1‑91 (MQ=255)

tttAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATc                                                                                           >  1:231325/1‑91 (MQ=255)

  tAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTGAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCaa                                                                                         >  1:3297306/1‑91 (MQ=255)

  tAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCaa                                                                                         <  1:2540058/91‑1 (MQ=255)

  tAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCaa                                                                                         <  1:2539916/91‑1 (MQ=255)

   agagTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCaaa                                                                                        <  1:910922/91‑1 (MQ=255)

     agTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAaata                                                                                      >  1:793776/1‑91 (MQ=255)

        gAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatt                                                                                   >  1:2637535/1‑91 (MQ=255)

          aTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                 <  1:3651805/91‑1 (MQ=255)

             cTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAg                                                                              <  1:251557/91‑1 (MQ=255)

              tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                             >  1:1364995/1‑91 (MQ=255)

               aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                            <  1:2407192/91‑1 (MQ=255)

                aTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTATTAGaaa                                                                           <  1:3229575/91‑1 (MQ=255)

                  aCTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                         <  1:2920223/91‑1 (MQ=255)

                  aCTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                         >  1:2027833/1‑91 (MQ=255)

                  aCTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                         >  1:366573/1‑91 (MQ=255)

                   cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                        <  1:2025533/91‑1 (MQ=255)

                     gaagcgcTAATTGAATCTTTATTTATTCACTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                      <  1:3326348/88‑1 (MQ=255)

                      aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                     <  1:3658222/91‑1 (MQ=255)

                       tGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCtt                                                                    <  1:684235/91‑1 (MQ=21)

                       tGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                    >  1:457077/1‑91 (MQ=255)

                       tGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                    <  1:1086693/91‑1 (MQ=255)

                        gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCttt                                                                   <  1:1828636/91‑1 (MQ=21)

                               ttGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCAtt                                                            <  1:2577869/91‑1 (MQ=255)

                                  aaTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGcc                                                         <  1:1542185/91‑1 (MQ=255)

                                   aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                        <  1:3106896/91‑1 (MQ=255)

                                   aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                        <  1:3106912/91‑1 (MQ=255)

                                    tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                       <  1:1214920/91‑1 (MQ=255)

                                    tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                       >  1:3717186/1‑91 (MQ=255)

                                    tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                       >  1:3415278/1‑91 (MQ=255)

                                     cTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaa                                                      >  1:1982017/1‑91 (MQ=255)

                                      tttatttaTTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaa                                                     <  1:2053554/91‑1 (MQ=255)

                                      tttatttaTTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaa                                                     <  1:2058002/91‑1 (MQ=255)

                                        tatttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaaa                                                   <  1:2008517/91‑1 (MQ=255)

                                         atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGGTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                  <  1:317469/91‑1 (MQ=255)

                                         atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                  <  1:317257/91‑1 (MQ=255)

                                         atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                  <  1:1534627/91‑1 (MQ=255)

                                          tttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACt                                                 <  1:1741682/91‑1 (MQ=255)

                                             attCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                              >  1:3005567/1‑91 (MQ=255)

                                             attCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                              <  1:3788951/91‑1 (MQ=255)

                                               tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAACTAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                            >  1:1224038/1‑91 (MQ=255)

                                               tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                            >  1:1224089/1‑91 (MQ=255)

                                                  gTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTcc                                         >  1:2893706/1‑91 (MQ=255)

                                                   tttCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTccc                                        <  1:3621262/91‑1 (MQ=255)

                                                    ttCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCg                                       <  1:2016344/91‑1 (MQ=255)

                                                    ttCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCg                                       <  1:860734/91‑1 (MQ=255)

                                                     tCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGa                                      >  1:3520323/1‑91 (MQ=255)

                                                      cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                     <  1:306457/91‑1 (MQ=255)

                                                      cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                     >  1:1730488/1‑91 (MQ=255)

                                                      cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                     >  1:1767723/1‑91 (MQ=255)

                                                       ataataAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAcc                                    >  1:2656893/1‑91 (MQ=255)

                                                       ataataAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAcc                                    <  1:161367/91‑1 (MQ=255)

                                                       ataataAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAcc                                    <  1:441591/91‑1 (MQ=255)

                                                          ataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctct                                 <  1:308733/91‑1 (MQ=255)

                                                          ataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctct                                 >  1:1548827/1‑91 (MQ=255)

                                                          ataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctct                                 >  1:3947106/1‑91 (MQ=255)

                                                          ataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctct                                 >  1:449219/1‑91 (MQ=255)

                                                           taaAATGCATATCAAAAAAGCTGCTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAGACGTAATGTCCCGACCTCtt                                <  1:2660442/91‑1 (MQ=255)

                                                           taaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCtt                                <  1:3659627/91‑1 (MQ=255)

                                                            aaaaTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCttt                               >  1:468105/1‑91 (MQ=255)

                                                            aaaaTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCttt                               >  1:469053/1‑91 (MQ=255)

                                                            aaaaTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCttt                               >  1:468746/1‑91 (MQ=255)

                                                              aaTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTAt                             <  1:3230551/91‑1 (MQ=255)

                                                                tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                           >  1:3404433/1‑91 (MQ=255)

                                                                tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                           <  1:2895883/91‑1 (MQ=255)

                                                                   atatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACg                        <  1:2371038/91‑1 (MQ=255)

                                                                    tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGc                       <  1:957826/91‑1 (MQ=255)

                                                                    tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGc                       >  1:1749198/1‑91 (MQ=255)

                                                                     atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCa                      >  1:490645/1‑91 (MQ=255)

                                                                       cAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATa                    <  1:960119/91‑1 (MQ=255)

                                                                          aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                 <  1:928209/91‑1 (MQ=255)

                                                                          aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                 <  1:928082/91‑1 (MQ=255)

                                                                           aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                >  1:1823623/1‑91 (MQ=255)

                                                                           aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                >  1:775949/1‑91 (MQ=255)

                                                                            aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc               <  1:3415125/91‑1 (MQ=255)

                                                                             aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACt              >  1:2320176/1‑91 (MQ=255)

                                                                              gCTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtt             >  1:3718958/1‑91 (MQ=25)

                                                                              gCTGATTTCGATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtt             <  1:3979193/91‑1 (MQ=255)

                                                                               cTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtta            >  1:1036612/1‑91 (MQ=25)

                                                                               cTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtta            >  1:1036613/1‑91 (MQ=25)

                                                                                tGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttat           <  1:1954620/91‑1 (MQ=255)

                                                                                tGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttat           >  1:78488/1‑91 (MQ=255)

                                                                                tGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttat           >  1:2954827/1‑91 (MQ=255)

                                                                                 gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt          <  1:3095573/91‑1 (MQ=255)

                                                                                 gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt          >  1:696829/1‑91 (MQ=255)

                                                                                 gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt          >  1:2853493/1‑91 (MQ=255)

                                                                                 gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt          >  1:3956852/1‑91 (MQ=255)

                                                                                  aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta         >  1:2070037/1‑91 (MQ=255)

                                                                                   tttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTAc        <  1:2709313/91‑1 (MQ=255)

                                                                                      cTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGa     <  1:2486434/91‑1 (MQ=255)

                                                                                      cTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGa     <  1:17860/91‑1 (MQ=255)

                                                                                       tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat    <  1:2244079/91‑1 (MQ=255)

                                                                                       tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat    <  1:3247912/91‑1 (MQ=255)

                                                                                       tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat    <  1:1873589/91‑1 (MQ=255)

                                                                                        aTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGata   <  1:1214794/91‑1 (MQ=255)

                                                                                        aTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGata   >  1:2906650/1‑91 (MQ=255)

                                                                                         tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat  >  1:803158/1‑91 (MQ=255)

                                                                                         tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat  >  1:4060618/1‑91 (MQ=255)

                                                                                         tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat  >  1:683411/1‑91 (MQ=255)

                                                                                         tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat  >  1:683266/1‑91 (MQ=255)

TTTAGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATAT  >  USA300TCH1516_ALE/2432147‑2432326

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