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%	5622.0 / ‑3.1	1172	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/0); major base . (529/642); minor base G (0/1); total (529/643)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 1245 total aligned reads displayed.

TCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAA  >  AM406671/1921593‑1921903

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

  gTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACACTTTTTGTGATATAA‑TTTTCCCAGaat                                                                                                                                                      >  1:814760/1‑148 (MQ=255)

  gTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaat                                                                                                                                                      <  1:2302540/149‑2 (MQ=255)

  gTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaat                                                                                                                                                      <  1:9541222/149‑2 (MQ=255)

    tCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaataa                                                                                                                                                    >  1:8625251/1‑146 (MQ=255)

    tCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaataa                                                                                                                                                    >  2:5160323/1‑146 (MQ=255)

     cAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAGaataaa                                                                                                                                                   <  1:1251349/149‑5 (MQ=255)

       aTTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAg‑aataaaat                                                                                                                                                <  1:2644801/149‑1 (MQ=255)

        ttttACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAg‑aataaaata                                                                                                                                               >  2:6197545/1‑149 (MQ=255)

          ttACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑Aataaaataaa                                                                                                                                             <  1:4813046/149‑1 (MQ=255)

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

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

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

              aacaaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGa                                                                                                                                         <  2:7298822/149‑1 (MQ=255)

                 aaaaaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCt                                                                                                                                      <  1:8799330/149‑1 (MQ=255)

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

                 aaaaaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCt                                                                                                                                      <  2:9816619/149‑1 (MQ=255)

                  aaaaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCtt                                                                                                                                     <  2:3135707/149‑1 (MQ=255)

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

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

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

                      ggACTTTTTGACAATATAAAC‑TTTCTTCTATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAA‑CTTAACA‑TCACCTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCt                                                                                                                                 >  2:6902030/1‑149 (MQ=255)

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

                       gACTTTTTGACAATATAAAC‑TTTCTTATNTTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtc                                                                                                                                >  1:5802032/1‑148 (MQ=255)

                        aCTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtcc                                                                                                                               >  1:5140405/1‑147 (MQ=255)

                          tttttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtccta                                                                                                                             <  2:5682921/149‑5 (MQ=255)

                          tttttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTCtccta                                                                                                                             >  1:9284501/1‑145 (MQ=255)

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

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

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

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

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

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

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

                                      atacactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAGGaaa                                                                                                                 <  2:10294598/143‑1 (MQ=255)

                                        aaactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCCACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACt                                                                                                              >  1:6563107/5‑149 (MQ=255)

                                        aaactttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACt                                                                                                              >  2:2721571/5‑149 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                     ctctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAtt                                                                                    >  1:4019047/4‑149 (MQ=255)

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

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

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

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

                                                                                 cacaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAt                                                                         >  1:1588019/2‑149 (MQ=255)

                                                                                  acaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAtt                                                                        <  2:7431954/149‑1 (MQ=255)

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

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

                                                                                         tcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTTA‑TTATAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTa                                                                  <  2:10050792/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                     gAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTTTTTTAT‑TTTTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTATTTTCAttttt                                                       >  2:7098177/1‑149 (MQ=255)

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

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

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

                                                                                                              taacatCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCGGAAGAACATTTTTTTAATTTCATTTTTTTATAtt                                                <  1:6925522/144‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                ctagaCAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCaaa                                <  2:2905268/145‑1 (MQ=255)

                                                                                                                                 tagaCAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAAt                               >  1:7765052/4‑149 (MQ=255)

                                                                                                                                    aCAATTTTCGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAAtctc                            <  1:6036153/149‑1 (MQ=255)

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

                                                                                                                                        tttttGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCatat                        <  2:8293602/149‑1 (MQ=255)

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

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

                                                                                                                                              tgATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTACGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAGGCTCAGATCTCATATTCAAgg                  >  2:4472195/1‑149 (MQ=255)

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

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

                                                                                                                                                   ataattttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGAc              <  1:2855537/145‑1 (MQ=255)

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

                                                                                                                                                     aattttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACaa            >  1:7959490/3‑149 (MQ=255)

                                                                                                                                                        tttCGCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCATAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAATGAGACAAAAt         <  2:2526951/149‑1 (MQ=255)

                                                                                                                                                         ttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATg        <  2:8708711/149‑1 (MQ=255)

                                                                                                                                                         ttCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATg        <  2:8992152/149‑1 (MQ=255)

                                                                                                                                                           cccAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATgaa      >  1:2152537/1‑149 (MQ=255)

TCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTTATATTAAGGGATAAGCTCAAATCTCATATTCAAGGAGACAAAATGAAGAAA  >  AM406671/1921593‑1921903

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