BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,432,243	A→T	5.8%	intergenic (+590/+57)	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,243	0	A	T	5.8%	363.4 / 10.6	104	intergenic (+590/+57)	birA_2/USA300HOU_RS12525	Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/58); new base T (3/3); total (43/61)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 115 total aligned reads displayed.

GAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAA  >  USA300TCH1516_ALE/2432155‑2432332

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

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

     cTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAg                                                                                    <  1:251557/91‑1 (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:2027833/1‑91 (MQ=255)

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

           cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                              <  1:2025533/91‑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: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)

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

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

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

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

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

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

                              tttatttaTTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaa                                                           <  1:2053554/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: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:306457/91‑1 (MQ=255)

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

                                              cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                           >  1:1730488/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)

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

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

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

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

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

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

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

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

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

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

                                                            tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGc                             <  1:957826/91‑1 (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:928209/91‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                         gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt                >  1:3956852/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)

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

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

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

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

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

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

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

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

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

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

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

                                                                                      taattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACga   <  1:1990304/91‑1 (MQ=255)

                                                                                       aattaatCGAAATCAGCTTTTTTCATTACCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa  <  1:2073739/91‑1 (MQ=255)

                                                                                       aattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa  <  1:665851/91‑1 (MQ=255)

                                                                                       aattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa  <  1:3458912/91‑1 (MQ=255)

                                                                                       aattaatAGAAATCAGCTTTTTTCATTACCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa  <  1:2073451/91‑1 (MQ=255)

                                                                                       aattaatAGAAATCAGCTTTTTTCATTACCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa  <  1:2073539/91‑1 (MQ=255)

GAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAA  >  USA300TCH1516_ALE/2432155‑2432332

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