BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,921,755	Δ1 bp	intergenic (‑245/‑140)	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,753	0	A	.	99.9%	4399.4 / ‑1.5	926	intergenic (‑243/‑142)	llmg_1941/glnP	conserved hypothetical protein/glutamine ABC transporter permease and substrate binding protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (426/497); total (427/499)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 979 total aligned reads displayed.

TCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAA  >  AM406671/1921593‑1921899

tCGTTCAATTTTACAACAAAAACGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAAGGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGa                                                                                                                                                    <  1:332726/149‑1 (MQ=255)

 cGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaa                                                                                                                                                   >  2:442196/1‑149 (MQ=255)

           tACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAAc                                                                                                                                        <  1:3982042/149‑1 (MQ=255)

             caacaaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTg                                                                                                                                      >  2:1927988/1‑149 (MQ=255)

             caacaaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTg                                                                                                                                      <  1:2728271/149‑1 (MQ=255)

                caaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGc                                                                                                                                   <  1:5969676/149‑1 (MQ=255)

                caaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGc                                                                                                                                   >  1:4569151/1‑149 (MQ=255)

                   aaaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGCTCTAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTg                                                                                                                                >  2:5331465/1‑149 (MQ=255)

                   aaaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTg                                                                                                                                >  2:4133672/1‑149 (MQ=255)

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

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

                         cTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtcct                                                                                                                          <  1:5474949/149‑4 (MQ=255)

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

                            tttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTaaa                                                                                                                      >  1:695662/1‑149 (MQ=255)

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

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

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

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

                                     ataAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAggaaa                                                                                                             >  1:2833357/1‑146 (MQ=255)

                                        aaactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACt                                                                                                          >  1:500297/5‑149 (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:4261257/1‑149 (MQ=255)

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

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

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

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

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

                                                         ttatcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTc                                                                                           <  1:5589931/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:7369216/145‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                  aaCTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTa                                                                                  >  1:2184989/1‑149 (MQ=255)

                                                                     ctctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCAT‑CTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAtt                                                                                >  2:6005430/4‑149 (MQ=255)

                                                                     ctctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAtt                                                                                <  2:1604007/146‑1 (MQ=255)

                                                                        tCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCag                                                                             >  1:8299161/1‑149 (MQ=255)

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

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

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

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

                                                                           aatttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaag                                                                          >  2:6364381/1‑149 (MQ=255)

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

                                                                            atttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCagaaga                                                                         <  1:311690/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                                                                      >  2:1123445/3‑149 (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                                                                 <  1:7668802/147‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                              aTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTc                                                         <  2:3975626/149‑1 (MQ=255)

                                                                                                 gtcagAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAt                                                       <  1:4822376/145‑1 (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                                                    <  2:4788032/148‑1 (MQ=255)

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

                                                                                                     gAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAttttt                                                   <  2:7679586/149‑1 (MQ=255)

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

                                                                                                      aaaaCCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTGAATGTA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCAtttttt                                                  <  2:5769577/149‑1 (MQ=255)

                                                                                                            cctaacatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTtata                                              <  2:3778340/142‑1 (MQ=255)

                                                                                                              taacatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATAtt                                            <  2:8783636/144‑1 (MQ=255)

                                                                                                                acatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTaa                                          <  2:2425319/146‑1 (MQ=255)

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

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

                                                                                                                   tCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAggg                                       <  1:5830358/149‑1 (MQ=255)

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

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

                                                                                                                       cTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATaa                                   >  2:2041237/1‑149 (MQ=255)

                                                                                                                           gattCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCt                                <  1:2367935/147‑1 (MQ=255)

                                                                                                                           gattCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCt                                <  1:949962/147‑1 (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                      <  2:7083525/149‑1 (MQ=255)

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

                                                                                                                                        tttttGTGATATAA‑TTTTCGCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAGTATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCatat                    <  1:4473734/149‑1 (MQ=255)

                                                                                                                                        tttttGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCatat                    >  1:464590/1‑149 (MQ=255)

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

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

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

                                                                                                                                            tgtgATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCaa                >  2:8337753/1‑149 (MQ=255)

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

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

                                                                                                                                                atatAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGgag            <  2:4246114/149‑1 (MQ=255)

                                                                                                                                                 tataA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGgaga           <  1:2584663/149‑1 (MQ=255)

                                                                                                                                                    taattttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACa         >  2:5217594/4‑149 (MQ=255)

                                                                                                                                                     aattttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACaa        <  1:7138712/147‑1 (MQ=255)

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

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

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

                                                                                                                                                          tCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATga   <  2:2795310/149‑1 (MQ=255)

TCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAA  >  AM406671/1921593‑1921899

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