BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	17.0%	I173I (ATT→ATA)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,660	0	T	A	17.0%	223.0 / 42.7	101	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (25/57); new base A (7/10); total (32/69)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

Only 100 of 106 total aligned reads displayed.

CCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTT  >  USA300TCH1516_ALE/110520‑110800

ccGGTTAAAATGGTCGCTAATTAAATCATTCACACAAAAGCAATCAAAGATGAACAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATTAGAAAGACGGAAATAtt                                                                                                                                              >  2:785657/1‑141 (MQ=17)

 cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                             <  2:308650/141‑1 (MQ=255)

             tCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                                 <  1:370304/141‑1 (MQ=255)

             tCGATAATAAAAACATTCCAACAAAAGAAAATAAAGATGAAACCATTCGCGAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                                 <  1:894193/141‑1 (MQ=17)

             tCGAAAATAAAATCAGTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                                 <  1:180889/141‑1 (MQ=255)

             gcGATAATAAAATCATACCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATCTAATCCag                                                                                                                                 <  2:82963/140‑1 (MQ=255)

              cGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                                                                                >  1:92632/1‑141 (MQ=255)

                ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCAAAAgg                                                                                                                              >  2:44860/1‑141 (MQ=255)

                 taataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGt                                                                                                                             >  2:269871/1‑141 (MQ=255)

                   acaaaaGCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                           <  1:847835/139‑1 (MQ=255)

                    taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                          >  1:653618/1‑141 (MQ=255)

                        aTCATTCCAACAAAAGAAATTAAAGAGGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                      <  2:187026/141‑1 (MQ=255)

                        aTCATTCAAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                      <  1:74959/141‑1 (MQ=255)

                           tcgacaacaaAAGATATTAAAGATGAAAAAAAAAAAAAAGAAATTCAAAACTTTACGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                                   <  2:49745/137‑1 (MQ=17)

                           aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                                   >  2:484109/1‑141 (MQ=255)

                               caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               <  1:494048/141‑1 (MQ=255)

                               caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               <  1:69173/141‑1 (MQ=255)

                               caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               <  1:45751/141‑1 (MQ=255)

                                  caaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                            <  1:105649/141‑1 (MQ=255)

                                  caaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAAACGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                            <  2:476542/141‑1 (MQ=255)

                                  aaaaatgaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                            <  2:187077/134‑1 (MQ=255)

                                  aaaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                            <  2:112033/140‑1 (MQ=255)

                                  aaaaaGAAATTAAAAATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATACTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                            <  1:432235/140‑1 (MQ=255)

                                   aaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAGAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                           <  2:280225/141‑1 (MQ=255)

                                   aaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                           <  1:755646/141‑1 (MQ=255)

                                   aaaaGAAATTAAAGACGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                           >  1:677284/1‑141 (MQ=255)

                                    aaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTAGGCTaga                                                                                                          >  2:751699/1‑139 (MQ=255)

                                       gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                       >  1:831131/1‑141 (MQ=255)

                                       aaaaCTAAAGATGAAACAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                       <  1:117151/140‑1 (MQ=255)

                                         aaTTAAAGATGAAAAAAAAAAAAAAAAAATCCAAAACTTTAAGTTACTTGTTCGATATAGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGTCTAAATATCa                                                                                                     <  2:323183/141‑1 (MQ=17)

                                                 aTGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATCTAATCCAGTGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                             <  1:867457/141‑1 (MQ=255)

                                                  tGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                                            <  1:15116/141‑1 (MQ=255)

                                                  tGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                                            <  1:808634/141‑1 (MQ=255)

                                                     aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATCACGCGGCAGCTACTTCATATAATCGAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                         <  2:192325/141‑1 (MQ=17)

                                                        aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCATTTAACTAAtgatga                                                                                      <  1:495687/141‑1 (MQ=255)

                                                        aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                      <  1:261123/141‑1 (MQ=255)

                                                         aTAAAAAAAGAAAGCGAGAACTTTAAGTTCGTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCAGATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                     <  1:54965/141‑1 (MQ=255)

                                                            aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                  <  2:455891/141‑1 (MQ=255)

                                                            aaaaaaGAAATCGAAAACCTTAAGGTATTGGTTCAAGATGGCGACTTAAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATACCAATTAACTAATGATGATtat                                                                                  <  1:344675/141‑1 (MQ=17)

                                                              aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTGTCTTAAGAGGATtatat                                                                                >  1:568321/1‑140 (MQ=17)

                                                                  gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                            <  2:16519/141‑1 (MQ=255)

                                                                    aaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAGCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                          <  1:380907/141‑1 (MQ=255)

                                                                             cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTACCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                 <  1:681464/141‑1 (MQ=255)

                                                                              tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                                >  2:386905/1‑141 (MQ=255)

                                                                                 aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGACTTTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTGACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                             >  1:856706/1‑141 (MQ=255)

                                                                                    ttCTTTGTTCAATATGGCGACTTTAAAACTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga                                                          <  1:701689/141‑1 (MQ=255)

                                                                                      cTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata                                                        <  1:406861/141‑1 (MQ=255)

                                                                                           ttCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                                   <  1:736317/141‑1 (MQ=255)

                                                                                           ttCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                                   <  2:100721/141‑1 (MQ=255)

                                                                                             cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAtata                                                 <  2:493469/141‑1 (MQ=255)

                                                                                                    gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                          >  1:726378/1‑141 (MQ=255)

                                                                                                    gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                          >  2:791974/1‑141 (MQ=255)

                                                                                                      gACTTAAAAAATTTGAAAAATTATAAAGACGGAGATTTCTAATATAATCCAGAGGGGCCGAGTTATTCGGCTAAATATCAATGAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGt                                        <  2:712236/141‑1 (MQ=18)

                                                                                                         tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                                     >  2:791217/1‑141 (MQ=255)

                                                                                                           tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAATACGCAAAAGATATGAGATACCGAAGAGTAAAGc                                   >  1:364532/1‑141 (MQ=255)

                                                                                                             aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                                 >  1:322206/1‑141 (MQ=255)

                                                                                                             aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCAACGAGTAAAGCTc                                 >  2:184333/1‑141 (MQ=255)

                                                                                                               aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCa                               <  1:484109/141‑1 (MQ=255)

                                                                                                                 tttaatgaaTTAAATAGATGGAGATATATAATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAACGATGATTATAATGTAAAACAACTACGCAAAAGATATGATATACCGACGAGTAAAGCTTCaaa                             <  1:651355/134‑1 (MQ=14)

                                                                                                                 tttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                             <  2:588022/141‑1 (MQ=255)

                                                                                                                 tttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                             <  1:269871/141‑1 (MQ=255)

                                                                                                                 tttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCAAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCACAAGATATGATATACCGACGAGTAAAGCTCCaaa                             <  1:270094/141‑1 (MQ=255)

                                                                                                                  ttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                            <  2:906979/135‑1 (MQ=255)

                                                                                                                   taatgtaTTAAAAAGATTGTGATATGTCAAATAATCCAGAGGTGCCCAGTTATTCGGCTAAATATCAATTAACTAATGATCATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                           <  2:802314/135‑1 (MQ=17)

                                                                                                                   taatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                           <  2:81192/136‑1 (MQ=255)

                                                                                                                       gaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTCTAGCTCCAAAGTTAtt                       <  1:453869/140‑1 (MQ=255)

                                                                                                                       aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAtt                       >  1:408809/1‑141 (MQ=255)

                                                                                                                           tataAAGACGGAGATCTTTCATATAACCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCTAAAGATATGATATACCGACGAGTCAAGCTCCAAAGTTATTGTTa                   <  1:552495/141‑1 (MQ=255)

                                                                                                                           tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa                   <  1:469819/141‑1 (MQ=255)

                                                                                                                            ataAAGACGGAGATATTTCATATAACCCAGAGGTGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa                  <  2:626959/141‑1 (MQ=255)

                                                                                                                             taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACAGACGAGTAAAGCTCCAAAGTTCTTGTTaaa                 >  2:123953/1‑141 (MQ=255)

                                                                                                                             aaaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTGCGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTCTTGTTaca                 <  2:672739/140‑3 (MQ=18)

                                                                                                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                >  2:210662/1‑141 (MQ=255)

                                                                                                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                >  1:889218/1‑141 (MQ=255)

                                                                                                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                >  2:599875/1‑141 (MQ=255)

                                                                                                                              aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATATCGACGAGTAAAGATCCAAAGTTATTGTTaaaa                >  2:459979/1‑141 (MQ=255)

                                                                                                                              aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                >  2:67797/1‑141 (MQ=255)

                                                                                                                              aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGATATTGTTaaaa                >  2:429911/1‑141 (MQ=255)

                                                                                                                               aatatggAGATATTTCATATAATTCAGAGGTGCCGAGTTAGTCGTCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACTCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  2:303820/136‑1 (MQ=18)

                                                                                                                               aaggtggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  1:524708/136‑1 (MQ=255)

                                                                                                                               aagaagGAGAAATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTCAATAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGTCTAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  2:402291/141‑1 (MQ=17)

                                                                                                                               aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  1:282578/141‑1 (MQ=255)

                                                                                                                               aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  1:492671/141‑1 (MQ=255)

                                                                                                                               aaGACGGAGATATTTCATATAATCCAGAGGTGCCGATTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  2:896035/141‑1 (MQ=255)

                                                                                                                               aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  1:5662/141‑1 (MQ=255)

                                                                                                                               aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               <  1:725951/141‑1 (MQ=255)

                                                                                                                                agatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg              >  1:767254/5‑141 (MQ=255)

                                                                                                                                  aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt            <  1:291879/141‑1 (MQ=255)

                                                                                                                                  aCGGAGAAATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt            <  1:291902/141‑1 (MQ=255)

                                                                                                                                   tggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc           <  1:516833/140‑1 (MQ=255)

                                                                                                                                     gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAg         >  1:201441/1‑141 (MQ=255)

                                                                                                                                     gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGATATTGTTAAAAGGTTCAg         >  2:415649/1‑141 (MQ=255)

                                                                                                                                       gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg       <  1:72457/141‑1 (MQ=255)

                                                                                                                                       gaTACTTCATACAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg       <  2:478546/141‑1 (MQ=255)

                                                                                                                                        atatatCATATAATCAAGAGGTGACGAGATATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGATATTGTTAAAAGGCTCAGGGa      >  2:882160/1‑141 (MQ=18)

                                                                                                                                         tataTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa     <  1:69438/141‑1 (MQ=255)

                                                                                                                                         tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa     <  1:690281/141‑1 (MQ=255)

                                                                                                                                         tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGATACCAAGTTATTGTTAAAAGGATCAGGGaa     >  2:43732/1‑141 (MQ=255)

                                                                                                                                          atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTCTTGTCAAAAGGTTCAGGGAAt    <  1:467890/141‑1 (MQ=255)

                                                                                                                                          atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAt    <  1:367392/141‑1 (MQ=255)

CCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTT  >  USA300TCH1516_ALE/110520‑110800

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