BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,934,780	A→T	6.0%	N156K (AAT→AAA)	USA300HOU_RS09605 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,934,780	0	A	T	6.0%	713.1 / 25.9	201	N156K (AAT→AAA)	USA300HOU_RS09605	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (73/116); new base T (4/8); total (77/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 209 total aligned reads displayed.

ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAAT  >  USA300TCH1516_ALE/1934690‑1934869

aCTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGaaa                                                                                           >  1:807838/1‑91 (MQ=255)

aCTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGaaa                                                                                           >  1:3430628/1‑91 (MQ=255)

     aGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaa                                                                                      >  1:2160130/1‑91 (MQ=255)

     aGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaa                                                                                      >  1:1249833/1‑91 (MQ=255)

          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaa                                                                                 <  1:760669/91‑1 (MQ=255)

          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaa                                                                                 <  1:759985/91‑1 (MQ=255)

          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaa                                                                                 >  1:721222/1‑91 (MQ=255)

          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaa                                                                                 <  1:2215531/91‑1 (MQ=255)

            tttcttttctTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAaataaata                                                                               <  1:1641969/91‑1 (MQ=255)

            tttcttttctTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAaataaata                                                                               <  1:3488572/91‑1 (MQ=255)

             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                              <  1:813286/91‑1 (MQ=255)

             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTCACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                              >  1:3595575/1‑91 (MQ=255)

              tcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATAtt                                                                             >  1:2312313/1‑91 (MQ=255)

               cttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTc                                                                            >  1:212892/1‑91 (MQ=255)

               cttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTc                                                                            >  1:4109880/1‑91 (MQ=255)

                ttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCt                                                                           >  1:1046188/1‑91 (MQ=255)

                ttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTATTTGAATTTCAAAAAAATATTCt                                                                           <  1:1575495/91‑1 (MQ=37)

                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTATTTGAATTTCAAAAAAATATTCtc                                                                          <  1:2183895/91‑2 (MQ=37)

                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTATTTGAATTTCAAAAAAATATTCtc                                                                          <  1:3155025/91‑2 (MQ=37)

                    cttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGc                                                                       <  1:1134956/91‑1 (MQ=255)

                     ttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCg                                                                      <  1:3998237/91‑1 (MQ=255)

                        aaaCTCTTCGTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAt                                                                   <  1:3973304/91‑1 (MQ=255)

                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                  <  1:1403949/91‑1 (MQ=255)

                          aCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAttt                                                                 >  1:3171266/1‑91 (MQ=255)

                          aCTCTTCGTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAttt                                                                 <  1:2067692/91‑1 (MQ=255)

                          aCTCTTCGTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAttt                                                                 <  1:2068478/91‑1 (MQ=255)

                           ctctTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtttt                                                                <  1:1871426/91‑1 (MQ=255)

                             cttcttCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGa                                                              >  1:2707728/1‑91 (MQ=255)

                               tcttcCAGTTCCAAATAGTTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAtt                                                            <  1:3549135/91‑1 (MQ=255)

                               tcttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAtt                                                            >  1:3577245/1‑91 (MQ=255)

                               tcttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAtt                                                            >  1:1895026/1‑91 (MQ=255)

                                         ccAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATc                                                  >  1:3126624/1‑91 (MQ=255)

                                         ccAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATc                                                  <  1:1225191/91‑1 (MQ=255)

                                            aataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaa                                               >  1:697232/1‑91 (MQ=255)

                                               aatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAat                                            >  1:1805659/1‑91 (MQ=255)

                                               aatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAat                                            >  1:1741300/1‑91 (MQ=255)

                                               aatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAat                                            <  1:1487119/91‑1 (MQ=255)

                                               aatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAat                                            <  1:3193053/91‑1 (MQ=255)

                                                attattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAata                                           <  1:459766/91‑1 (MQ=255)

                                                 ttattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAatat                                          <  1:2093987/91‑1 (MQ=255)

                                                 ttattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAatat                                          <  1:534290/91‑1 (MQ=255)

                                                 ttattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAatat                                          <  1:1472313/91‑1 (MQ=255)

                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATc                                         <  1:1620511/91‑1 (MQ=255)

                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                       <  1:3925003/91‑1 (MQ=255)

                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                       >  1:2776676/1‑91 (MQ=255)

                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                       >  1:4061949/1‑91 (MQ=255)

                                                     tGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtt                                      >  1:3750799/1‑91 (MQ=255)

                                                      gACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAttt                                     >  1:334389/1‑91 (MQ=255)

                                                       aCTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtttt                                    <  1:3728040/91‑1 (MQ=255)

                                                       aCTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtttt                                    <  1:2956102/91‑1 (MQ=255)

                                                        cTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTc                                   <  1:564350/91‑1 (MQ=255)

                                                          gATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGa                                 >  1:1007925/1‑91 (MQ=255)

                                                          gATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGa                                 <  1:307264/91‑1 (MQ=255)

                                                           aTTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGat                                <  1:3133244/91‑1 (MQ=255)

                                                           aTTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGat                                <  1:2653749/91‑1 (MQ=255)

                                                           aTTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGat                                <  1:3123353/91‑1 (MQ=255)

                                                            ttGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGata                               <  1:1282810/91‑1 (MQ=255)

                                                            ttGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGata                               <  1:259783/91‑1 (MQ=255)

                                                             tGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGatat                              >  1:2392029/1‑91 (MQ=255)

                                                             tGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGatat                              <  1:32586/91‑1 (MQ=255)

                                                               ttCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATaa                            >  1:461171/1‑91 (MQ=255)

                                                                tCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAc                           >  1:930242/1‑91 (MQ=255)

                                                                tCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAc                           >  1:1934902/1‑91 (MQ=255)

                                                                tCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAc                           >  1:67324/1‑91 (MQ=255)

                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                         <  1:1838/91‑1 (MQ=255)

                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                         <  1:1351948/91‑1 (MQ=255)

                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                         <  1:3701710/91‑1 (MQ=255)

                                                                   aaaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtt                        <  1:2577378/91‑1 (MQ=255)

                                                                   aaaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtt                        <  1:2254408/91‑1 (MQ=255)

                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                       <  1:1446275/91‑1 (MQ=255)

                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                       >  1:4142099/1‑91 (MQ=255)

                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                       <  1:2433988/91‑1 (MQ=255)

                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                       <  1:2260426/91‑1 (MQ=255)

                                                                     aTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtttt                      <  1:2533595/91‑1 (MQ=255)

                                                                     aTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtttt                      <  1:2153543/91‑1 (MQ=255)

                                                                         gagaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAAtt                  <  1:3040393/91‑1 (MQ=255)

                                                                          agaATATTTATTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTATTAATTc                 >  1:2080939/1‑91 (MQ=255)

                                                                          agaATATTTATTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTATTAATTc                 >  1:2080977/1‑91 (MQ=255)

                                                                            aaTATTTTTTTGCAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtt               >  1:3198286/1‑91 (MQ=255)

                                                                            aaTATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtt               <  1:2358174/91‑1 (MQ=255)

                                                                             atatTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCttt              <  1:2996681/91‑1 (MQ=255)

                                                                              tatTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtttt             <  1:109974/91‑1 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            >  1:879542/1‑91 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            <  1:3303138/91‑1 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            <  1:3273523/91‑1 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            <  1:1668054/91‑1 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            <  1:1655202/91‑1 (MQ=255)

                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa            <  1:2540631/91‑1 (MQ=255)

                                                                                tttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAg           >  1:586318/1‑91 (MQ=255)

                                                                                tttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAg           <  1:681625/91‑1 (MQ=255)

                                                                                 ttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGa          >  1:2903173/1‑91 (MQ=255)

                                                                                  tttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGat         <  1:2380406/91‑1 (MQ=255)

                                                                                   ttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGata        <  1:1413614/91‑1 (MQ=255)

                                                                                    tttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat       <  1:653361/91‑1 (MQ=255)

                                                                                    tttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat       <  1:248413/91‑1 (MQ=255)

                                                                                    tttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat       <  1:1170924/91‑1 (MQ=255)

                                                                                     ttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATg      <  1:1711989/91‑1 (MQ=255)

                                                                                      tGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGa     <  1:2014157/91‑1 (MQ=255)

                                                                                       gAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGaa    <  1:1331354/91‑1 (MQ=255)

                                                                                         aaTTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAat  >  1:306286/1‑91 (MQ=255)

ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAAT  >  USA300TCH1516_ALE/1934690‑1934869

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