BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,512	A→G	27.4%	S194S (TCA→TCG)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,512	0	A	G	27.4%	197.8 / 80.1	123	S194S (TCA→TCG)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (45/44); new base G (12/22); total (57/66)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 124 total aligned reads displayed.

CTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111373‑111650

cTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAg                                                                                                                                            <  1:273277/141‑1 (MQ=255)

cTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAg                                                                                                                                            <  2:207847/141‑1 (MQ=255)

  aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATAATATAATCAAGAGGTGCCGAGTTATTCGGCt                                                                                                                                          >  2:395918/1‑141 (MQ=255)

  aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATAAATCATATAAACCAAAGGTACCGAGTTAATCaacc                                                                                                                                          >  1:54573/1‑138 (MQ=25)

   agaaGGAACTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCta                                                                                                                                         <  2:349900/138‑2 (MQ=255)

     gAGGAAGTAAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATAGAATCCGGAGGTGCCGAGTTATTCGGCTaaa                                                                                                                                       >  2:1014189/1‑141 (MQ=21)

      aGGAAGTTAAAGTAGAAAAAATAAAAAAAGAAATCGAATACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAATATGGAGGTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                                      <  2:204052/141‑1 (MQ=21)

      aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                                      <  1:608534/141‑1 (MQ=255)

       ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCtata                                                                                                                                     <  2:388075/141‑1 (MQ=255)

       ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTGTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCacaa                                                                                                                                     <  2:11471/141‑4 (MQ=255)

           gTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCGTATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                                                 <  2:788060/141‑1 (MQ=255)

            tcaaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTACGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                                >  2:486636/3‑141 (MQ=255)

             tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                               <  2:88411/141‑1 (MQ=255)

             tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                               <  1:203831/141‑1 (MQ=255)

             tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGGTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                               <  2:214075/141‑1 (MQ=255)

             tAAAGAAGAAAACATAAAAAAAGAAACCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                               <  1:573054/141‑1 (MQ=255)

              aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                              >  1:935279/1‑141 (MQ=255)

              aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                                                              <  1:1028286/141‑1 (MQ=255)

               aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                             <  2:876126/141‑1 (MQ=255)

               aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATAAAGAAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATACAATCAAGAGTTGCCGAGTTATTCAGCTCAATAACAATTaa                                                                                                                             >  2:207267/1‑141 (MQ=21)

                agaagtAAACATAACAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                                                            <  1:266582/141‑1 (MQ=255)

                agaagaAAACATAAAAAAAGAAATTGAAAACTTGAAGTTTTTTGCGCAATATAGAAATTTTAAAGATTTAATGAATTATAAAGATGTAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                                                            <  1:1044008/141‑1 (MQ=17)

                agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCCATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTGATTCGGCTAAATATCAATTAAc                                                                                                                            <  2:917896/141‑1 (MQ=255)

                 gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                           <  2:463258/141‑1 (MQ=255)

                 gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                           >  1:567985/1‑141 (MQ=255)

                 gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                           <  1:663908/141‑1 (MQ=255)

                    gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTCGCTATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGTTCCCGAGTTATTCATCTTAATATCAATTAACTaaa                                                                                                                        >  2:572349/1‑140 (MQ=18)

                    aaaaaaaaaaaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                        <  2:1038318/140‑1 (MQ=255)

                    aaaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                                        <  1:1020392/140‑1 (MQ=255)

                        aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                    >  1:703625/1‑141 (MQ=255)

                         cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                                                   >  2:291947/1‑141 (MQ=255)

                            aaaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                                                                >  1:760637/1‑141 (MQ=255)

                             aaaaaaGCAATGCAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGTTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAGTAACTAATGTTGATtag                                                                                                               <  2:298738/141‑2 (MQ=21)

                             aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCTATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCATAGGTACCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                               >  1:758514/1‑141 (MQ=255)

                             aaaaaaGAAAACAAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                               <  2:383367/141‑1 (MQ=255)

                              aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtata                                                                                                              <  2:130895/141‑1 (MQ=255)

                               aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                                                                             <  1:622077/141‑1 (MQ=255)

                                aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCCTATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                            >  2:761654/1‑141 (MQ=255)

                                aaaGAAATCGAAAACTTTAAGTTTTCTGCGTAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGCTTATAAt                                                                                                            >  2:440501/1‑141 (MQ=17)

                                  aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                          <  2:515134/141‑1 (MQ=255)

                                    aattCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATGCAGAGGTTCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                        <  1:461559/138‑1 (MQ=17)

                                    aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                        <  1:461343/141‑1 (MQ=255)

                                    aaaGCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCACTATCATTTAACTAATGATGATTATAATGTaa                                                                                                        <  2:486942/141‑1 (MQ=255)

                                      aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                      >  2:309011/1‑141 (MQ=255)

                                      aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                      >  2:142564/1‑141 (MQ=255)

                                       tCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACAAATGGTGATTATAATTTAAAAc                                                                                                     >  2:23623/1‑141 (MQ=18)

                                        cGAAAACTCTAAATTTTTTGCGCAATATAGCAATTTCAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACCAATGATGATTATAATGTAAAACa                                                                                                    <  1:883641/141‑1 (MQ=17)

                                          aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                  >  2:545933/1‑141 (MQ=255)

                                          aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                  <  1:141911/141‑1 (MQ=255)

                                                ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGCTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAAAAtaaatca                                                                                            >  1:85272/1‑135 (MQ=18)

                                                  aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCTATTAACTAATGATGATTATAATGTAAAACAATTACGCaa‑a                                                                                         <  1:158206/141‑1 (MQ=255)

                                                   aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa‑aa                                                                                        >  2:342013/1‑141 (MQ=255)

                                                   aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa‑aa                                                                                        >  2:341873/1‑141 (MQ=255)

                                                      ttttttCTCAATACAGCAACTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAA‑AAGat                                                                                     <  2:654391/141‑1 (MQ=14)

                                                      tttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTAAGAAGCTGGAGATATATCATATAATCCAGAAGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAA‑AAGat                                                                                     <  2:323694/141‑1 (MQ=17)

                                                       ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAA‑AAGata                                                                                    <  2:915762/141‑1 (MQ=255)

                                                            cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACa                                                                               <  2:798314/141‑1 (MQ=255)

                                                             gccatatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTATGTAA‑AAGATATGACAt                                                                              <  1:526655/138‑1 (MQ=255)

                                                             gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACAt                                                                              <  1:758401/141‑1 (MQ=255)

                                                              cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAA‑AAGATATGAtata                                                                             <  1:25600/141‑1 (MQ=14)

                                                              cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATa                                                                             <  2:62395/141‑1 (MQ=255)

                                                               aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAA‑AAGATATGATATAc                                                                            <  2:605969/141‑1 (MQ=14)

                                                                 tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGCGGTGCCATGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACca                                                                          <  1:1043017/141‑1 (MQ=255)

                                                                 tagaGCAATTGTAAAGACTTAATGAAGTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACca                                                                          <  2:680394/138‑1 (MQ=255)

                                                                      taaTTTTAAAGATTTAAGGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACaaa                                                                     <  2:80303/140‑1 (MQ=255)

                                                                          tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGAGATGACATACCAACAAATaaa                                                                 >  1:645987/1‑141 (MQ=255)

                                                                          tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAACTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATaaa                                                                 <  1:761654/141‑1 (MQ=255)

                                                                             aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGAAATACCAACAAATAAAgcg                                                              >  2:60217/1‑141 (MQ=255)

                                                                                gagttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCg                                                           <  1:159677/138‑1 (MQ=255)

                                                                                gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATTCCAACAAATAAAGCGCCg                                                           >  1:633993/1‑141 (MQ=255)

                                                                                     aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACACATAAAGCGCCGAAGCt                                                      <  1:89250/141‑1 (MQ=255)

                                                                                       tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                    >  2:104189/1‑141 (MQ=255)

                                                                                       tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                    <  2:8740/141‑1 (MQ=255)

                                                                                         aaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGctgctg                                                  >  1:169855/1‑141 (MQ=255)

                                                                                         aaTTAAAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                                  <  1:486636/141‑1 (MQ=255)

                                                                                          aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACTTAA‑AAGATATGAAATACCAACAATTAAAGAGACGAGGAtgatgt                                                 >  2:700453/1‑141 (MQ=17)

                                                                                               aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGaaa                                            >  2:41846/1‑141 (MQ=255)

                                                                                                aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAg                                           >  2:664874/1‑141 (MQ=255)

                                                                                                aaGATGGAGATATAGCATATAATCCAGAGGTGCCGAGTTATTCAAATCAATAACAAGTAACTAATGATGATTATAATGTACAAAAATTACGTAA‑AAGATATGACATAACTACAAAGAAAGCGCCGAAGCTGTTGTTGAAAg                                           >  2:433744/1‑141 (MQ=17)

                                                                                                    tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAc                                       <  2:886209/141‑1 (MQ=255)

                                                                                                      gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                     >  1:214207/1‑141 (MQ=255)

                                                                                                       agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    >  1:253722/1‑141 (MQ=255)

                                                                                                       agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    >  1:658539/1‑141 (MQ=255)

                                                                                                        aatatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAC‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTTAAAGGTACAggg                                   <  2:358382/140‑1 (MQ=255)

                                                                                                          gataaCATATAAGCCAGAGGTGCCGAGTTATTCAGCGCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGaa                                 <  2:434130/140‑1 (MQ=255)

                                                                                                           atatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGAGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                >  2:743695/1‑141 (MQ=255)

                                                                                                             atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATTATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                              <  2:366067/141‑1 (MQ=255)

                                                                                                               cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaa                            >  1:4877/1‑141 (MQ=255)

                                                                                                                     aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                      >  1:349579/1‑141 (MQ=255)

                                                                                                                     aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                      >  1:36867/1‑141 (MQ=255)

                                                                                                                      tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGCTGATTAAAATGAAAAACAATTACGTAATACTGAATCACGT‑GGGACGAGAAAATCTCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtc                     >  2:747909/1‑141 (MQ=9)

                                                                                                                       ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGCTATGACATACCAACAAATAAAGCGCCGAAGCTGTTATTGAAAGGCACAGGGAATTTAAAAGGTtca                    >  2:55609/1‑141 (MQ=255)

                                                                                                                        cAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcat                   <  1:86518/141‑1 (MQ=255)

                                                                                                                           agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAg                >  2:848749/1‑141 (MQ=255)

                                                                                                                             gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt              >  1:62889/1‑141 (MQ=255)

                                                                                                                             gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGATTCATCAGtg              >  1:102532/1‑140 (MQ=255)

                                                                                                                             gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATATAAAAGGTTCATCAGtt              >  2:820166/1‑141 (MQ=255)

                                                                                                                             gTGCAGAGTTATTCACGTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt              <  2:664541/141‑1 (MQ=255)

                                                                                                                                 cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGat          >  2:122270/1‑141 (MQ=255)

                                                                                                                                       aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaa    >  1:636774/1‑141 (MQ=255)

CTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAA‑AAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111373‑111650

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