BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,921,790	Δ1 bp	intergenic (‑280/‑105)	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,788	0	G	.	99.2%	4430.1 / 24.9	941	intergenic (‑278/‑107)	llmg_1941/glnP	conserved hypothetical protein/glutamine ABC transporter permease and substrate binding protein
Reads supporting (aligned to +/- strand): ref base G (4/4); new base . (410/521); total (415/526)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 989 total aligned reads displayed.

CAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCATGTTTGC  >  AM406671/1921626‑1921932

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

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

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

  atatAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAgga                                                                                                                                                >  2:6769160/1‑148 (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:6582404/149‑3 (MQ=255)

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

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

        aactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTa                                                                                                                                          >  2:3955896/4‑149 (MQ=255)

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

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

          ctttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTaaa                                                                                                                                        <  2:1813658/148‑1 (MQ=255)

          ctttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTaaa                                                                                                                                        >  1:8449049/2‑149 (MQ=255)

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

           tttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTaaaa                                                                                                                                       >  2:4671247/1‑149 (MQ=255)

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

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

                       tttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTAtt                                                                                                                             <  1:1255682/142‑1 (MQ=255)

                       tttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTAtt                                                                                                                             <  1:3916647/142‑1 (MQ=255)

                        ttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTNAATATA‑GGAAACTAAAAGAATGTTATTc                                                                                                                            <  1:5084323/143‑1 (MQ=255)

                         tatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTcc                                                                                                                           >  1:6381958/6‑149 (MQ=255)

                          atcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCg                                                                                                                          >  1:2861500/5‑149 (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                                                                                                                          >  1:6604889/5‑149 (MQ=255)

                           tcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGa                                                                                                                         >  1:6196947/4‑149 (MQ=255)

                           tcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGa                                                                                                                         <  1:4345857/146‑1 (MQ=255)

                            cagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAt                                                                                                                        >  1:1123669/3‑149 (MQ=255)

                             agagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAtt                                                                                                                       <  1:6778839/148‑1 (MQ=255)

                              gagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAttt                                                                                                                      >  2:4812654/1‑149 (MQ=255)

                              gagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAttt                                                                                                                      >  2:59267/1‑149 (MQ=255)

                               agaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAtttt                                                                                                                     >  1:7237523/1‑149 (MQ=255)

                               agaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGAtttt                                                                                                                     >  1:7616680/1‑149 (MQ=255)

                                 aaCTC‑TCAAATTTC‑ACAAA‑TTCAAAATGCCA‑GAGAACCT‑AACA‑TCAACTTGA‑TGCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGTATGTTATTCCGATTTTTa                                                                                                                   >  2:5342049/1‑149 (MQ=255)

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

                                   actctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAt                                                                                                                  >  2:8280058/5‑149 (MQ=255)

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

                                     tctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTc                                                                                                                >  1:4078766/3‑149 (MQ=255)

                                     tctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTc                                                                                                                <  2:7576257/147‑1 (MQ=255)

                                     tctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTc                                                                                                                <  1:1507151/147‑1 (MQ=255)

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

                                          aatttc‑acaaa‑ttcagaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaag                                                                                                           <  2:7166679/149‑1 (MQ=255)

                                           atttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaaga                                                                                                          <  1:98528/149‑1 (MQ=255)

                                              ttcacaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAAc                                                                                                        <  1:7187010/146‑1 (MQ=255)

                                              ttaacaaa‑ttcagaATGTCA‑GAAAACCT‑AAC‑GTCATCTTGA‑TTCTAGAA‑GACTTTTGTGGTGTAA‑TTTTCTCAG‑AATAAAATAAACTGATC‑TGGGTGCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTGATTCAGAAGAAc                                                                                                        <  2:3527755/147‑1 (MQ=255)

                                                    aaattcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAtttt                                                                                                   <  1:5503819/146‑1 (MQ=255)

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

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

                                                          aaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAAt                                                                                             >  2:3675012/1‑149 (MQ=255)

                                                          aaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAAt                                                                                             <  2:5081766/149‑1 (MQ=255)

                                                            aaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAAttt                                                                                           <  2:8887851/149‑1 (MQ=255)

                                                              tGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCa                                                                                         <  2:4513835/149‑1 (MQ=255)

                                                                 tcagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGCTTTTTATTCAGAAGAACATTTTTTTAATTTCAtt                                                                                       >  1:7017614/4‑149 (MQ=255)

                                                                  cagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATCTA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTGATTTCAttt                                                                                      <  1:4756079/147‑1 (MQ=255)

                                                                  cagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttt                                                                                      >  2:1448019/3‑149 (MQ=255)

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

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

                                                                     aaaaCCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAtttttt                                                                                   >  1:5264611/1‑149 (MQ=255)

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

                                                                      aaaCCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttt                                                                                      <  2:1706998/145‑1 (MQ=255)

                                                                      aaaCCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttt                                                                                      >  1:1706998/1‑145 (MQ=255)

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

                                                                                 atCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAgg                                                                         <  1:7744235/148‑1 (MQ=255)

                                                                                 atCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAgg                                                                         >  1:6092404/2‑149 (MQ=255)

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

                                                                                   cAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGa                                                                       <  1:8787317/149‑1 (MQ=255)

                                                                                        ttgattCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAg                                                                   <  1:4628056/145‑1 (MQ=255)

                                                                                            ttCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCa                                                               >  1:808503/1‑149 (MQ=255)

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

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

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

                                                                                                    cAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCa                                                        <  1:8622795/149‑1 (MQ=255)

                                                                                                     aaTTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCat                                                       >  1:3008373/1‑149 (MQ=255)

                                                                                                     aaTTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCat                                                       <  2:3421384/149‑1 (MQ=255)

                                                                                                          ttGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAATGGATAAGCTCAAATCTCATATTCa                                                  >  2:1999236/1‑149 (MQ=255)

                                                                                                               atatAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATATGCTCAAATCTCATATTCAAGgag                                             >  2:2550066/1‑149 (MQ=255)

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

                                                                                                                              agaataaaataaACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATgaagaa                                <  2:7355421/147‑1 (MQ=255)

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

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

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

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

                                                                                                                                   aaaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAattat                           <  2:2656005/149‑1 (MQ=255)

                                                                                                                                    aaataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAattatt                          <  1:8931223/149‑1 (MQ=255)

                                                                                                                                     aataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAACTCTCATATTCAAGGAGACAAAATGAAGAAATTAttt                         <  2:6165875/149‑1 (MQ=255)

                                                                                                                                     aataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAttt                         <  1:3767013/149‑1 (MQ=255)

                                                                                                                                     aataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAttt                         <  2:5469626/149‑1 (MQ=255)

                                                                                                                                     aataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAttt                         <  1:4945430/149‑1 (MQ=255)

                                                                                                                                      ataaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTAtttt                        <  1:5941516/149‑1 (MQ=255)

                                                                                                                                        aaaCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTc                      >  1:5623434/1‑149 (MQ=255)

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

                                                                                                                                          aCTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGc                    <  2:4127757/149‑1 (MQ=255)

                                                                                                                                             taGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGctct                 <  2:6907147/148‑1 (MQ=255)

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

                                                                                                                                                 ttGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCa             <  2:6436463/149‑1 (MQ=255)

                                                                                                                                                     tctCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTc          <  1:315636/147‑1 (MQ=255)

                                                                                                                                                      ctCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCa         >  1:2650966/2‑149 (MQ=255)

                                                                                                                                                       tCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCAt        <  2:3170250/149‑1 (MQ=255)

                                                                                                                                                         cTAAATATA‑GGAAACTAAAAGAATGTTATGCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCATGt      >  2:5657186/1‑149 (MQ=255)

                                                                                                                                                          tAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCATGtt     <  2:2129494/149‑1 (MQ=255)

                                                                                                                                                            aaTATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCATGTTTg   <  1:585744/149‑1 (MQ=255)

CAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAATTATTTTTCGCTCTGGCACTCATGTTTGC  >  AM406671/1921626‑1921932

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