BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,921,704	Δ1 bp	intergenic (‑194/‑191)	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,702	0	A	.	99.7%	4731.4 / 5.2	998	intergenic (‑192/‑193)	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 . (484/509); minor base G (2/1); total (488/510)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1058 total aligned reads displayed.

AGTTCTGTTTACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTT  >  AM406671/1921547‑1921852

   tCTGTTTACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCTaaca                                                                                                                                                     >  1:404076/1‑147 (MQ=255)

   tCTGTTTACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTCTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCTaaca                                                                                                                                                     >  2:7856603/1‑147 (MQ=255)

      gTTTACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCTaacatca                                                                                                                                                  <  1:6301004/149‑6 (MQ=255)

       tttACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCTaacatcaa                                                                                                                                                 <  2:5598604/149‑7 (MQ=255)

       tttACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATGTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAGAACCTaacatcaa                                                                                                                                                 >  1:5013432/1‑143 (MQ=255)

        ttACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCTaacatcaac                                                                                                                                                >  2:84683/1‑142 (MQ=255)

           cTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTg                                                                                                                                           >  2:2584663/1‑149 (MQ=255)

            ttAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGa                                                                                                                                          >  2:9164895/1‑149 (MQ=255)

            ttAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGa                                                                                                                                          <  1:4548667/149‑1 (MQ=255)

            ttAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGa                                                                                                                                          >  2:6335938/1‑149 (MQ=255)

            ttAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGa                                                                                                                                          >  2:3429943/1‑149 (MQ=255)

                 tAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGAttcta                                                                                                                                     >  2:6505823/1‑146 (MQ=255)

                  aGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGAttctag                                                                                                                                    >  2:3767013/1‑145 (MQ=255)

                   gAATACGTAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGa                                                                                                                                  >  2:4473734/1‑149 (MQ=255)

                    aaTACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACCATATGAAC‑TTTCTTATATCTTA‑TCA‑GAGAACTC‑TCAAATCTC‑ACAAA‑TTCACAATGTCAG‑‑AAACCCCAACA‑TCAACTTGA‑TTCTAGac                                                                                                                                 >  1:8400868/1‑148 (MQ=255)

                    aaTACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGac                                                                                                                                 >  1:5368661/1‑148 (MQ=255)

                     aTACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGaca                                                                                                                                <  1:8583478/149‑3 (MQ=255)

                     aTACATAAATATCTCTATCAAGCACTCGTTCAATTTAACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGaca                                                                                                                                >  2:6160339/1‑147 (MQ=255)

                      tACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGacaa                                                                                                                               >  2:8665899/1‑146 (MQ=255)

                       aCATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAGacaat                                                                                                                              >  1:2096957/1‑145 (MQ=255)

                               atCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGa                                                                                                                     <  2:3984192/149‑1 (MQ=255)

                                tctcTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGat                                                                                                                    >  2:2158556/1‑149 (MQ=255)

                                  tctATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGatat                                                                                                                  >  2:1127411/1‑149 (MQ=255)

                                    tATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCGACTTGA‑TTCTAG‑ACAATTTTTGTGATATaa                                                                                                                <  1:2555036/149‑1 (MQ=255)

                                     aTCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAAC‑CCTCAAATTTC‑AGAAA‑TTCAAAATGTCA‑GAGAAGCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGAtatact                                                                                                               >  2:6640814/1‑147 (MQ=255)

                                      tCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAAtt                                                                                                              <  1:8856575/149‑1 (MQ=255)

                                      tCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAAtt                                                                                                              <  2:162911/149‑1 (MQ=255)

                                          gCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAAttttcc                                                                                                          <  1:4868364/149‑3 (MQ=255)

                                          gCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATCTTTGTGATATAAttttcc                                                                                                          >  2:2677413/1‑147 (MQ=255)

                                           cACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAATTTTccc                                                                                                         >  1:734140/1‑149 (MQ=255)

                                            aCTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAAttttccca                                                                                                        <  2:3561377/149‑5 (MQ=255)

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

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

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

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

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

                                                             acaaAAAGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAg                                                                                     <  1:4441262/149‑1 (MQ=255)

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

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

                                                                  aaGGACTTTTTGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGAttcta                                                                                                                                     <  2:773217/100‑4 (MQ=255)

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

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

                                                                         ttttGACAATATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTaa                                                                        <  1:981360/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)

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

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

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

                                                                                aaTATAAAC‑TTTCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATAgg                                                                 >  1:6189449/1‑149 (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                                                                <  2:4926106/149‑2 (MQ=255)

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

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

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

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

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

                                                                                          tttCTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCACCTTGA‑TTCTAG‑ACAATTTTTGTGATAT‑ACTTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAACCTaaaa                                                       >  2:4650272/1‑149 (MQ=255)

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

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

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

                                                                                             cTTATATTTTA‑TCA‑GAGAACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGaa                                                    <  2:2050008/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)

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

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

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

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

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

                                                                                                         atcagagaACTC‑TCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAat                                                                      >  2:5457968/5‑122 (MQ=255)

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

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

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

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

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

                                                                                                                  actctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTCG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATACA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAt                                  <  2:7533943/145‑1 (MQ=255)

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

                                                                                                                    tctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTTTGAT‑TTATTTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTATATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTAttt                                >  2:8990790/3‑148 (MQ=255)

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

                                                                                                                     ctCAAATTTC‑ACAAA‑TTCAAAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCa                               <  2:2355689/148‑1 (MQ=255)

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

                                                                                                                       caaatttc‑acaaa‑ttcaaaaTGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCaga                             <  1:3215023/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:6562444/149‑1 (MQ=255)

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

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

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

                                                                                                                             ttcacaaa‑ttcaaaATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCA‑TAATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAAc                        >  2:5736252/4‑149 (MQ=255)

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

                                                                                                                                    aattcaaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACAttttt                  <  1:7668934/147‑1 (MQ=255)

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

                                                                                                                                        caaAATGTCA‑GAAAACCT‑AACA‑TCAACTTGA‑TTCTAG‑ACAATTTTTGTGATATAA‑TTTTCCCAG‑AATAAAATAAACTGAGCTTGTC‑TCCTAAATATA‑GGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTaa              <  1:1359454/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

AGTTCTGTTTACTTAAATAGAATACATAAATATCTCTATCAAGCACTCGTTCAATTTTACAACAAAAAGGACTTTTTGACAATATAAACTTTTCTTATATTTTATTCAGGAGAACTCTTCAAATTTCAACAAATTTCAAAATGTCAGGAAAACCTAAACATTCAACTTGATTTCTAGAACAATTTTTGTGATATAATTTTTCCCAGAAATAAAATAAACTGAGCTTGTCTTCCTAAATATAGGGAAACTAAAAGAATGTTATTCCGATTTTTATTCAGAAGAACATTTTTTTAATTTCATTTTTTT  >  AM406671/1921547‑1921852

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