BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	667,064	T→G	9.0%	D18E (GAT→GAG)	USA300TCH1516_00600 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	667,064	0	T	G	9.0%	732.4 / 52.7	231	D18E (GAT→GAG)	USA300TCH1516_00600	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (96/114); new base G (17/4); total (113/118)
Fisher's exact test for biased strand distribution p-value = 2.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTGAT  >  USA300TCH1516_ALE/666974‑667154

gCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGAt                                                                                            >  2:2496572/1‑91 (MQ=255)

  tAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATgg                                                                                          <  1:1032213/91‑1 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGt                                                                                         >  1:3438107/1‑91 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGt                                                                                         <  2:1822933/91‑1 (MQ=255)

    aaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTc                                                                                        <  2:683733/91‑1 (MQ=255)

       ttGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATa                                                                                     <  1:1704895/91‑1 (MQ=255)

          ttATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAg                                                                                  >  1:1446463/1‑91 (MQ=255)

            aTCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGt                                                                                <  1:1130901/91‑1 (MQ=255)

               gTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACa                                                                             >  1:1499305/1‑91 (MQ=255)

               gTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACa                                                                             >  1:1699290/1‑91 (MQ=255)

                  aaaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAat                                                                          <  1:146346/91‑1 (MQ=255)

                  aaaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAat                                                                          <  1:2884312/91‑1 (MQ=255)

                   aaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAata                                                                         <  1:732117/91‑1 (MQ=255)

                    atatGGAGGTTATATCAAAGACTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAatat                                                                        <  2:3495894/91‑1 (MQ=255)

                     tatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATAtt                                                                       >  2:1709115/1‑91 (MQ=255)

                      atGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAACAATGACCCTCTAGATGGTCATAAAGGTACAAATATTg                                                                      <  1:3345930/91‑1 (MQ=255)

                      atGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTg                                                                      <  2:2447016/91‑1 (MQ=255)

                      atGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTg                                                                      >  1:390342/1‑91 (MQ=255)

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTATATGGTCATAAAGGTACAAATATTgg                                                                     >  1:3671824/1‑91 (MQ=255)

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTAAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTgg                                                                     <  2:1005034/91‑1 (MQ=255)

                        ggAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGt                                                                    <  2:3357036/91‑1 (MQ=255)

                         gAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGtt                                                                   >  1:2514467/1‑91 (MQ=255)

                         gAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGtt                                                                   >  2:754472/1‑91 (MQ=255)

                         gAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGtt                                                                   <  1:3230895/91‑1 (MQ=255)

                          aGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTc                                                                  <  1:2547759/91‑1 (MQ=255)

                           ggTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCa                                                                 <  1:963624/91‑1 (MQ=255)

                           ggTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCa                                                                 <  2:1144255/91‑1 (MQ=255)

                            gTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCaa                                                                >  1:2209341/1‑91 (MQ=255)

                             ttATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAAt                                                               >  2:3853817/1‑91 (MQ=255)

                                tatCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAg                                                            <  1:2349509/91‑1 (MQ=255)

                                taacATTTTATGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAg                                                            <  2:3508163/88‑1 (MQ=255)

                                  tCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCt                                                          >  2:3754176/1‑91 (MQ=255)

                                   cATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCtt                                                         <  1:511268/91‑1 (MQ=255)

                                   cATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCtt                                                         <  2:2557150/91‑1 (MQ=255)

                                     ttGTGTGAGAAAAAAGGCTTTAATTTTAATATCATAATAAATGACCCTCTGGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       <  2:2282120/91‑1 (MQ=255)

                                     ttGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       <  1:3572216/91‑1 (MQ=255)

                                     ttGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       <  2:2336235/91‑1 (MQ=255)

                                     ttGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       <  2:2869112/91‑1 (MQ=255)

                                      tGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGa                                                      <  2:211702/91‑1 (MQ=255)

                                       gTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGAc                                                     <  2:2358521/91‑1 (MQ=255)

                                        tCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACa                                                    >  1:3591630/1‑91 (MQ=255)

                                        tCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACa                                                    <  1:1975961/91‑1 (MQ=255)

                                           gagaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAata                                                 >  2:1629100/1‑91 (MQ=255)

                                            agaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAatat                                                >  1:780048/1‑91 (MQ=255)

                                             gaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATAtt                                               >  2:371799/1‑91 (MQ=255)

                                             gaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATAtt                                               <  1:3831747/91‑1 (MQ=255)

                                              aaaaaaGGCTTTAATTTTAATATCATAAAAAGTGACCCTCTAGAGGGTCATAAGGGTGCAAATATTGGTTCAATTAGCTTAGACAATATTg                                              <  1:3160525/91‑1 (MQ=255)

                                                 aaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAc                                           >  1:821821/1‑91 (MQ=255)

                                                  aaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGTTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAcc                                          <  2:2786349/91‑1 (MQ=255)

                                                  aaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAcc                                          >  1:1228622/1‑91 (MQ=255)

                                                  aaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAcc                                          >  1:2174570/1‑91 (MQ=255)

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        >  2:85417/1‑91 (MQ=255)

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        >  2:525490/1‑91 (MQ=255)

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        >  2:455781/1‑91 (MQ=255)

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        <  2:307052/91‑1 (MQ=255)

                                                     gCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGt                                       <  1:2118797/91‑1 (MQ=255)

                                                      cTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtt                                      >  1:1022134/1‑91 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAGTATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    <  2:854817/91‑1 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGCTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    <  2:3572366/91‑1 (MQ=255)

                                                          aattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTa                                  <  1:1791048/91‑1 (MQ=255)

                                                          aattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTa                                  >  1:988983/1‑91 (MQ=255)

                                                           atttaaagATAAAAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTAt                                 <  2:3519361/91‑1 (MQ=255)

                                                            ttataataaCATAAAAAATGACCCTCTAGATGGTCATAAAAGTACAAATATTGGTTCAATTAACTTAGACAATATTGCACCAGTTTTTATc                                <  1:1180002/88‑1 (MQ=255)

                                                                aatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATg                            <  2:3248203/91‑1 (MQ=255)

                                                                 atatCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGt                           <  1:2759092/91‑1 (MQ=255)

                                                                  tatCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGtt                          >  1:429744/1‑91 (MQ=255)

                                                                  tatCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGtt                          <  2:2920815/91‑1 (MQ=255)

                                                                    tCATAAGAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGc                        <  2:1424472/91‑1 (MQ=255)

                                                                    tCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGc                        >  2:1412216/1‑91 (MQ=255)

                                                                     cATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCt                       >  1:1408153/1‑91 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                      >  1:1833229/1‑91 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                      <  2:3279946/91‑1 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                      <  2:2148103/91‑1 (MQ=255)

                                                                        aaaaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAAc                    <  2:1655350/91‑1 (MQ=255)

                                                                        aaaaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAAc                    <  1:226111/91‑1 (MQ=255)

                                                                         aaaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACa                   >  1:2656600/1‑91 (MQ=255)

                                                                          aaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACaa                  <  2:412209/91‑1 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTCGACAATATTGCACCCGTTCTTATCGATGTTGCTAACaca                 >  1:2475802/1‑89 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACaaa                 >  1:645662/1‑91 (MQ=255)

                                                                            aaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGATAATATTGCACCAGTTTTTATCGATGTTGCTAACAaag                <  1:3344614/91‑1 (MQ=255)

                                                                              tGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaa              <  2:2415984/91‑1 (MQ=255)

                                                                                aCCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGc            <  2:572795/91‑1 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACATTATTGCACCAGTTTTTATTGATGTTGCTAACAAAGAAGCa           >  1:393918/1‑91 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa           >  1:2013020/1‑91 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa           <  2:734523/91‑1 (MQ=255)

                                                                                  catctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt          >  2:2837667/3‑91 (MQ=255)

                                                                                  catctAGAGGGTCATAAAGGTACAAATATAGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt          >  2:3451186/3‑91 (MQ=255)

                                                                                   ctctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:3703583/1‑91 (MQ=255)

                                                                                   ctctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  2:525677/1‑91 (MQ=255)

                                                                                   atctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:3284946/2‑91 (MQ=255)

                                                                                   atctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:240254/2‑91 (MQ=255)

                                                                                   atctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:744107/2‑91 (MQ=255)

                                                                                    tctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAttt        <  2:2660376/91‑1 (MQ=255)

                                                                                    tctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAttt        >  1:1630771/1‑91 (MQ=255)

                                                                                     ctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTa       <  2:1187045/91‑1 (MQ=255)

                                                                                      tAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGAATTTAt      <  1:3075678/91‑1 (MQ=255)

                                                                                       aGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAtt     <  1:2494191/91‑1 (MQ=255)

                                                                                       aGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAtt     <  2:1987294/91‑1 (MQ=255)

                                                                                       aGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAtt     >  1:564566/1‑91 (MQ=255)

                                                                                        gATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTg    <  1:692447/91‑1 (MQ=255)

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTGAT  >  USA300TCH1516_ALE/666974‑667154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`