BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,912,184	T→G	32.7%	N167K (AAT→AAG)	USA300HOU_RS15290 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,912,184	0	T	G	32.7%	62.1 / 57.9	65	N167K (AAT→AAG)	USA300HOU_RS15290	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (33/10); new base G (10/11); total (43/22)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

AAAGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGT  >  USA300TCH1516_ALE/1912045‑1912323

aaaGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAtt                                                                                                                                             <  1:1040713/141‑1 (MQ=255)

  aGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACCATGCATGAAACCAAATTCAAAATCCAATAGGAAAATCATGAAGAATtat                                                                                                                                           >  1:299601/1‑141 (MQ=255)

  aGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAAGGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAtaca                                                                                                                                           >  2:35578/1‑138 (MQ=255)

   gATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATtata                                                                                                                                          >  2:618516/1‑141 (MQ=255)

   gATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGtata                                                                                                                                          >  2:1045666/1‑141 (MQ=255)

      atGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATaagt                                                                                                                                       >  1:87847/1‑139 (MQ=255)

           tACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGaa                                                                                                                                  <  1:497489/141‑1 (MQ=255)

           aaCCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGaa                                                                                                                                  <  1:662137/140‑1 (MQ=255)

             ccGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGAACAAATACATGAATCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATa                                                                                                                                >  1:274205/1‑141 (MQ=14)

              cGTAATGATAGATATGGTTTTAAACGGGACTTCAAACTATATGAATGCGATGACCGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATAc                                                                                                                               <  2:641291/141‑1 (MQ=255)

                tAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAACTTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACtt                                                                                                                             >  1:83185/1‑141 (MQ=255)

                  aTGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTa                                                                                                                           >  1:958280/1‑141 (MQ=255)

                             ggTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAATCCAAAtt                                                                                                                >  1:1032117/1‑141 (MQ=255)

                                  tAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCa                                                                                                           >  1:853168/1‑141 (MQ=255)

                                  tAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCa                                                                                                           <  1:1192838/141‑1 (MQ=255)

                                       gagaCTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑c                                                                                                     >  2:973206/1‑141 (MQ=255)

                                        cgaCTTCAAACTATATGAATGCGCTGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGACTCTTTTCATGTCCAAATTAATCAAAAG‑Ct                                                                                                    >  2:243642/2‑141 (MQ=14)

                                         gaCTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑Ctt                                                                                                   >  1:743038/1‑141 (MQ=255)

                                             tCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAACTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTg                                                                                               <  2:429913/141‑1 (MQ=255)

                                             tCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTATAAGGCACCAATGAATGAAACAAAATACAAAATCCAATAAGAAAAACAGCAAAAATTAGAACCCGGAAAACATTAAAGTCCAAATAAATCAAAAG‑CTTTCta                                                                                               >  2:280988/1‑140 (MQ=255)

                                              cAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGa                                                                                              <  1:182517/141‑1 (MQ=255)

                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGaa                                                                                       >  1:503322/1‑141 (MQ=255)

                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGaa                                                                                       >  1:1106420/1‑141 (MQ=255)

                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTATAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGaa                                                                                       <  1:209885/141‑1 (MQ=255)

                                                       tGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATAAAAAG‑CTTTCTGAACAAGaaaa                                                                                     >  2:774313/1‑140 (MQ=255)

                                                            gCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGaaaa                                                                                <  1:293496/141‑1 (MQ=255)

                                                             cGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGaaaaa                                                                               >  2:998919/1‑141 (MQ=255)

                                                               aTGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGaaaaaaa                                                                             <  1:1185927/141‑1 (MQ=255)

                                                               aTGACTCGTCAGCATGTTCTTTGAGGCAACAATGCATGAAACAAAATTCAAAATCCAATAAGAAAATCATGAAAAATTATAATTGGGAATACTTTGAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAACCGaaaaaaa                                                                             >  2:305962/1‑141 (MQ=14)

                                                                 gACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCACTAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTACAAATAAAACAAAAG‑CTTTCTGAACCAGAAAAGaaaaagaac                                                                           >  2:841993/1‑137 (MQ=14)

                                                                 gACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATc                                                                           >  1:634779/1‑141 (MQ=255)

                                                                 gACTGTTCACCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCGAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATc                                                                           >  2:770177/1‑141 (MQ=255)

                                                                  aCTGTTCAGCATGTTCTTTGAGGCAACACTGCTCGAAACCAAATTCAAAATCCAATAAGAAAATGATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCt                                                                          >  2:10585/1‑141 (MQ=14)

                                                                  aCTGTTCAGCATGTGCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGCGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCt                                                                          <  1:334648/141‑1 (MQ=255)

                                                                      ttCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAg                                                                      >  1:714988/1‑141 (MQ=255)

                                                                          gCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCaa                                                                  <  1:430092/141‑1 (MQ=255)

                                                                             tGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAT‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGa                                                               >  1:923410/1‑141 (MQ=255)

                                                                             tGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATGGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTAAAAGa                                                               >  1:754390/1‑141 (MQ=255)

                                                                              gTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaa                                                              >  2:1034240/1‑141 (MQ=39)

                                                                               ttCTTTGAGGCAACAATGCGTGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGGATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaaa                                                             >  1:718472/1‑141 (MQ=255)

                                                                               ttCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaaa                                                             >  1:110947/1‑141 (MQ=38)

                                                                                   ttGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAAtt                                                         >  1:971663/1‑141 (MQ=39)

                                                                                   ttGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAAtt                                                         >  2:564372/1‑141 (MQ=39)

                                                                                   ttGAGGCAACAATGCATGAAAACAAATTCAAAATACAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACTAAAAAAATCTATAGTCAAAGAAAAAtt                                                         <  2:1108149/141‑1 (MQ=255)

                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCCTGAAGACTTATAATTGGGAATACTTTCAAGTCC‑AAGCACTCCAGAGTCGTTCTATAGCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGAta                                                     >  2:843177/1‑140 (MQ=11)

                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATg                                                     >  1:422967/1‑141 (MQ=37)

                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATg                                                     >  2:957534/1‑141 (MQ=37)

                                                                                        gCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGt                                                    <  1:1095157/141‑1 (MQ=25)

                                                                                           acaaTGCATGAAACCAAATTCAAAATCCGATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAATTATCTATAGTCAAAGAAAAATTGATGTaga                                                 <  2:383047/141‑1 (MQ=12)

                                                                                             aaTGCATGAAACCAAATTCAAAATCAAAGAAGAAAATCATGAAGAAGTATAAGTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACAAAAAAAAAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGc                                               <  2:87847/141‑1 (MQ=11)

                                                                                               tGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAAAGGCt                                             >  2:1141559/1‑141 (MQ=17)

                                                                                                gCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTg                                            <  1:1045666/141‑1 (MQ=17)

                                                                                                 cATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGt                                           <  1:956947/141‑1 (MQ=32)

                                                                                                    gAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCACAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGtttt                                        <  1:398903/141‑1 (MQ=17)

                                                                                                    gAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGtttt                                        <  1:1101268/141‑1 (MQ=17)

                                                                                                      aaCCCACTTCAAAACCCAATAAGGAAATCATGAAGAAGGATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTg                                      <  1:1169909/141‑1 (MQ=11)

                                                                                                       aCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTgg                                     >  1:644331/1‑141 (MQ=32)

                                                                                                       aCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAAAGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTGTTTgg                                     >  2:150913/1‑141 (MQ=17)

                                                                                                         cAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGAt                                   >  1:1023118/1‑141 (MQ=32)

                                                                                                            aTTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAACCTATAGTCAAAGAAAAACGGATGTAGAGCCTGTTTTGGGATTTa                                >  2:393580/1‑141 (MQ=255)

                                                                                                              tCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATg                              >  2:602978/1‑141 (MQ=17)

                                                                                                                                aTCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑ATATCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAGATTGGTGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTTGGTTTa            >  1:996914/1‑141 (MQ=11)

                                                                                                                                  cATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGGTTAAc          >  2:1061251/1‑141 (MQ=17)

                                                                                                                                   aTGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAAGCTATTTTTGGTTTAACt         >  1:400355/1‑141 (MQ=12)

                                                                                                                                     gaagaaTTATAATTGGAAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCg       >  1:592895/1‑141 (MQ=17)

                                                                                                                                       agaaGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGaa     <  1:108894/141‑1 (MQ=17)

                                                                                                                                          aTTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGt  <  2:379662/141‑1 (MQ=35)

AAAGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAG‑CTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGT  >  USA300TCH1516_ALE/1912045‑1912323

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