BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,921,777	Δ1 bp	intergenic (‑267/‑118)	llmg_1941 ← / → glnP	conserved hypothetical protein/glutamine ABC transporter permease and substrate binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM406671	1,921,776	0	T	.	99.6%	4572.1 / 17.1	976	intergenic (‑266/‑119)	llmg_1941/glnP	conserved hypothetical protein/glutamine ABC transporter permease and substrate binding protein
Reads supporting (aligned to +/- strand): ref base T (0/3); major base . (435/530); minor base C (3/1); total (439/537)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 1010 total aligned reads displayed.

GGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCA  >  AM406671/1921615‑1921921

ggACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCt                                                                                                                                                   <  1:784944/149‑1 (MQ=255)

 gACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtc                                                                                                                                                  <  2:5078509/149‑2 (MQ=255)

  aCTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtcc                                                                                                                                                 <  1:1303262/149‑3 (MQ=255)

     ttttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTaa                                                                                                                                             <  1:8683319/149‑1 (MQ=255)

     ttttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTaa                                                                                                                                             <  1:3371065/149‑1 (MQ=255)

       ttGACAATATAAAC‑TTGCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAat                                                                                                                                           <  1:6183767/149‑1 (MQ=255)

        tGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAata                                                                                                                                          <  2:887492/149‑1 (MQ=255)

        tGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAata                                                                                                                                          >  2:5329687/1‑149 (MQ=255)

        tGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAata                                                                                                                                          <  1:3026500/149‑1 (MQ=255)

             atatAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAgga                                                                                                                                     >  1:2430353/1‑148 (MQ=255)

              tataAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAggaa                                                                                                                                    <  2:6246029/149‑3 (MQ=255)

                taAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAGGaaaa                                                                                                                                  <  1:3204574/149‑1 (MQ=255)

                taAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACATGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAAtgtatgaacc                                                                                                                                  >  2:2458527/1‑143 (MQ=255)

                  aaactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACt                                                                                                                                <  2:8353399/145‑1 (MQ=255)

                    actttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTaa                                                                                                                              >  2:6506398/3‑149 (MQ=255)

                    actttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTaa                                                                                                                              <  1:1005215/147‑1 (MQ=255)

                       ttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAg                                                                                                                           <  1:881094/149‑1 (MQ=255)

                        tCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTA‑TACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGa                                                                                                                          >  2:3208482/1‑149 (MQ=255)

                         cTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGaa                                                                                                                         <  1:2868021/149‑1 (MQ=255)

                          ttATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAAt                                                                                                                        >  1:487597/1‑149 (MQ=255)

                          ttATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTCGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAAt                                                                                                                        >  2:1804066/1‑149 (MQ=255)

                            atatTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGt                                                                                                                      <  2:5209232/149‑1 (MQ=255)

                             tatTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGtt                                                                                                                     >  1:962609/1‑149 (MQ=255)

                             tatTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGtt                                                                                                                     >  1:2656005/1‑149 (MQ=255)

                                  tttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTAtt                                                                                                                  >  1:7302178/8‑149 (MQ=255)

                                  tttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTC‑CGAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATACAATAAACTGAGCTTTTC‑TCCTAAATATA‑GGAAACTAAAAGAATGCTAtt                                                                                                                  >  2:4287756/8‑149 (MQ=255)

                                     atcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAat                                                                                                                                           <  1:5457968/118‑1 (MQ=255)

                                     atcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCg                                                                                                               <  2:8728157/145‑1 (MQ=255)

                                     atcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCg                                                                                                               >  2:3440362/5‑149 (MQ=255)

                                       cagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAt                                                                                                             <  1:6054828/147‑1 (MQ=255)

                                           gaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAttttt                                                                                                         >  2:1964517/1‑149 (MQ=255)

                                              actctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGTTATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAt                                                                                                       >  1:5767921/5‑149 (MQ=255)

                                              actctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAt                                                                                                       <  2:2096957/145‑1 (MQ=255)

                                                   caaatttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCaga                                                                                                  <  1:2764198/149‑1 (MQ=255)

                                                    aaatttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaa                                                                                                 >  2:4795554/1‑149 (MQ=255)

                                                    aaatttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaa                                                                                                 <  2:7186137/149‑1 (MQ=255)

                                                         ttcacaaa‑ttcaaaATGTC‑CGAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑AGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAAc                                                                                             <  1:3827684/146‑1 (MQ=255)

                                                          tcacaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACa                                                                                            <  1:3629991/147‑1 (MQ=255)

                                                           cacaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAt                                                                                           <  2:734140/148‑1 (MQ=255)

                                                                aattcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAttttt                                                                                       >  2:3786927/3‑149 (MQ=255)

                                                                aattcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAttttt                                                                                       >  2:8517571/3‑149 (MQ=255)

                                                                 attcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAtttttt                                                                                      <  2:754135/148‑1 (MQ=255)

                                                                  ttcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAttttttt                                                                                     <  1:1145674/149‑1 (MQ=255)

                                                                  ttcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAttttttt                                                                                     <  1:5544201/149‑1 (MQ=255)

                                                                   tcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTa                                                                                    >  1:225756/1‑149 (MQ=255)

                                                                      aaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAAtt                                                                                 >  1:7622049/1‑149 (MQ=255)

                                                                      aaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAAtt                                                                                 >  1:8011432/1‑149 (MQ=255)

                                                                           gtcagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAt                                                                             >  2:360409/5‑149 (MQ=255)

                                                                             cagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttt                                                                           <  1:4125973/147‑1 (MQ=255)

                                                                              agAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAtttt                                                                          <  1:7211418/148‑1 (MQ=255)

                                                                              agAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCACAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAtttt                                                                          <  2:4661873/148‑1 (MQ=255)

                                                                               gAAAACCT‑AACA‑TCAACTTCA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttttt                                                                         <  1:4503715/149‑1 (MQ=255)

                                                                                 aaaCCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttttttt                                                                       <  1:2051726/149‑1 (MQ=255)

                                                                                         aacatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTa                                                                 <  2:1354109/145‑1 (MQ=255)

                                                                                          acatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTaa                                                                >  2:7802283/4‑149 (MQ=255)

                                                                                           catCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAg                                                               >  1:1712357/3‑149 (MQ=255)

                                                                                             tCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAggg                                                             >  2:2335215/1‑149 (MQ=255)

                                                                                              cAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGa                                                            >  2:2660533/1‑149 (MQ=255)

                                                                                                aCTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATa                                                          <  2:8278992/149‑1 (MQ=255)

                                                                                                        tCTAG‑ACAATTTTTTTTTTTTAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAATGGATAAGCTCaa                                                   >  2:6379648/1‑149 (MQ=255)

                                                                                                        tCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCaa                                                   <  1:6475278/149‑1 (MQ=255)

                                                                                                           tagaCAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATCCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCANTTTTTTATATTAAGGGATAAGCTCAAAt                                                 <  1:5058577/146‑1 (MQ=255)

                                                                                                            agaCAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAAtt                                                >  2:5622302/3‑148 (MQ=255)

                                                                                                             gaCAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAAtct                                               <  1:8974422/148‑1 (MQ=255)

                                                                                                                 aTTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCata                                           <  1:59267/149‑1 (MQ=255)

                                                                                                                   ttttGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATAtt                                         >  2:5758013/1‑149 (MQ=255)

                                                                                                                    tttGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTc                                        <  2:5444119/149‑1 (MQ=255)

                                                                                                                       gtgATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAg                                     <  2:4156989/149‑1 (MQ=255)

                                                                                                                       gtgATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAg                                     >  1:6911879/1‑149 (MQ=255)

                                                                                                                         gATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGga                                   >  1:185478/1‑149 (MQ=255)

                                                                                                                          atatAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGgag                                  >  1:8484807/1‑149 (MQ=255)

                                                                                                                          atatAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGgag                                  >  1:4623938/1‑149 (MQ=255)

                                                                                                                                attttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACaaa                             <  2:668347/148‑1 (MQ=255)

                                                                                                                                 ttttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACaaaa                            >  2:1614738/1‑149 (MQ=255)

                                                                                                                                  tttACCAG‑AATAAAATAAACTGTGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAAAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAAt                           <  2:1127103/149‑1 (MQ=255)

                                                                                                                                   ttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATg                          >  2:2979209/1‑149 (MQ=255)

                                                                                                                                     cccAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATgaa                        <  1:7092125/149‑1 (MQ=255)

                                                                                                                                       ccagaataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGCAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATgaag                       <  1:6684829/145‑1 (MQ=255)

                                                                                                                                        cagaataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATgaaga                      <  2:3733820/146‑1 (MQ=255)

                                                                                                                                         agaataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATCTTCAAGGAGACAAAATgaagaa                     >  1:5247868/3‑149 (MQ=255)

                                                                                                                                         agaataaaagaaACTGGGCTGGGC‑TCCTAAATATA‑GGAAACTGAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATGTTTGTATATGAAGGGATAAGCGCAAATCTCAGATGCAAGGAGACAAACTgaagaa                     <  2:4953077/147‑1 (MQ=255)

                                                                                                                                          gaataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGaaa                    <  2:705222/148‑1 (MQ=255)

                                                                                                                                           aataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAat                   >  1:913032/1‑149 (MQ=255)

                                                                                                                                           aataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAat                   >  1:8909408/1‑149 (MQ=255)

                                                                                                                                            ataaaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAATTGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAatt                  >  2:7669451/1‑149 (MQ=255)

                                                                                                                                            ataaaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAatt                  >  2:6863396/1‑149 (MQ=255)

                                                                                                                                            ataaaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAatt                  <  2:404076/149‑1 (MQ=255)

                                                                                                                                             taaaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAatta                 >  2:8130902/1‑149 (MQ=255)

                                                                                                                                               aaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAattatt               >  2:798047/1‑149 (MQ=255)

                                                                                                                                               aaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAattatt               >  1:4697595/1‑149 (MQ=255)

                                                                                                                                               aaataaaCTGAGCATAT‑ATCCTAAATA‑AATGAAACTAAAAGAATGTTATTCAGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGAGATAAGCTCAAATCTCATATTCAAAGAGACAAAATGAAGAAattatt               <  2:9051931/149‑1 (MQ=255)

                                                                                                                                                aataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAttt              >  1:626972/1‑149 (MQ=255)

                                                                                                                                                aataaaCTGAGCTTGTC‑TCCAAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAattatat              >  2:6087483/1‑147 (MQ=255)

                                                                                                                                                 ataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAtttt             >  1:3986012/1‑149 (MQ=255)

                                                                                                                                                  taaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCatat                                          >  1:8342612/1‑119 (MQ=255)

                                                                                                                                                   aaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTc           <  1:8198905/149‑1 (MQ=255)

                                                                                                                                                      cTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGct        <  1:7495274/149‑1 (MQ=255)

                                                                                                                                                       tGAGCTTGTC‑TCCTAAATATA‑AGAAACTAAAAAAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAATGAGACAAAATGAAGAAATTATTTTTCGctc       <  1:1548475/149‑1 (MQ=255)

                                                                                                                                                         aGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTg     <  2:449578/149‑1 (MQ=255)

                                                                                                                                                           cTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGc   >  1:4302711/1‑149 (MQ=255)

GGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCA  >  AM406671/1921615‑1921921

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