BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,432,246	A→G	5.7%	intergenic (+593/+54)	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,246	0	A	G	5.7%	348.1 / 8.2	105	intergenic (+593/+54)	birA_2/USA300HOU_RS12525	Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/59); new base G (3/3); total (43/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 = 8.70e-01

Only 100 of 115 total aligned reads displayed.

ATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAAGAA  >  USA300TCH1516_ALE/2432157‑2432336

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)

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

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

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

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

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

                          tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                                 <  1:1214920/91‑1 (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:3788951/91‑1 (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)

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

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

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

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

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

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

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

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

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

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

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

                                                  aaaaTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCttt                                         >  1:469053/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)

                                                      tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                                     >  1:3404433/1‑91 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                 aaaTAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATaa          <  1:3649624/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:665834/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)

                                                                                         aatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaaagaa  >  1:2876857/1‑91 (MQ=255)

ATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAAGAA  >  USA300TCH1516_ALE/2432157‑2432336

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