BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,432,240	A→T	5.7%	intergenic (+587/+60)	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,240	0	A	T	5.7%	365.6 / 10.5	106	intergenic (+587/+60)	birA_2/USA300HOU_RS12525	Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/59); new base T (3/3); total (44/62)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 110 total aligned reads displayed.

AGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACG  >  USA300TCH1516_ALE/2432150‑2432330

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

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

     gAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatt                                                                                       <  1:796601/91‑1 (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)

              tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                              <  1:253990/91‑1 (MQ=25)

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

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

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

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

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

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

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

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

                    tGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                        >  1:595202/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)

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

                                aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                            <  1:3106896/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:3788951/91‑1 (MQ=255)

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

                                          attCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                                  >  1:3005567/1‑91 (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:860734/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                 tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGc                           <  1:3036382/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)

                                                                      aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaaa                      <  1:2745655/91‑1 (MQ=255)

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

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

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

                                                                        aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                    <  1:14291/91‑1 (MQ=255)

                                                                         aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                   <  1:1751038/91‑1 (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:1036613/1‑91 (MQ=25)

                                                                            cTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtta                >  1:1036612/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:696829/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                      tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat      <  1:2205068/91‑1 (MQ=255)

                                                                                        aaaTAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATaa    <  1:3649624/91‑1 (MQ=255)

                                                                                          ataatTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACg  <  1:62692/91‑1 (MQ=255)

AGAGTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACG  >  USA300TCH1516_ALE/2432150‑2432330

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