BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,340,125	0	A	G	42.4%	73.5 / 115.6	100	intergenic (+66/‑265)	ribAB‑II/PP_3814	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand): ref base A (26/31); new base G (18/24); total (45/55)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 106 total aligned reads displayed.

CCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCC  >  NC_002947/4339976‑4340273

ccgccgCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGa                                                                                                                                                      <  1:5285223/150‑1 (MQ=255)

ccgccgCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGa                                                                                                                                                      <  1:1097760/150‑1 (MQ=255)

ccgccgCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGa                                                                                                                                                      <  2:1319336/150‑1 (MQ=255)

ccgccgCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTGAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGa                                                                                                                                                      <  2:5745815/150‑1 (MQ=255)

 cgccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGac                                                                                                                                                     <  1:2838299/150‑1 (MQ=255)

   ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                   <  2:1630924/150‑1 (MQ=255)

    cgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAcacc                                                                                                                                                  <  2:710629/150‑1 (MQ=255)

      cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                <  2:836457/150‑1 (MQ=255)

      cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                <  1:2199368/150‑1 (MQ=255)

       tGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACaccacc                                                                                                                                               <  1:5526380/150‑1 (MQ=255)

              cgcgGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAt                                                                                                                                        <  1:1719457/150‑1 (MQ=255)

               gcgGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCCCGCATGCGAGGGTGTCACCACCGACGCagt                                                                                                                                       >  2:3981308/1‑148 (MQ=255)

               gcgGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                       >  2:3981661/1‑150 (MQ=255)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:332714/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  1:4897911/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  1:6023036/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  1:798558/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:1167770/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:1638148/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  1:2874857/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:4676470/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                    <  1:4998402/150‑1 (MQ=255)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                    <  2:1195635/150‑1 (MQ=255)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                    <  1:642617/150‑1 (MQ=255)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTAATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:786289/150‑1 (MQ=25)

                  gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTAATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                    <  2:786280/150‑1 (MQ=25)

                   ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                   <  2:2318163/150‑1 (MQ=25)

                   ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                                   <  2:2461925/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  2:3778302/150‑1 (MQ=25)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  2:3619162/150‑1 (MQ=25)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  1:4591818/150‑1 (MQ=25)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  2:5701952/150‑1 (MQ=25)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  2:40501/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  2:1637360/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  2:4495474/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:1168596/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:1197414/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:1888959/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:2077493/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:2077502/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:338848/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                  <  1:3531904/150‑1 (MQ=255)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCAGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  1:3627432/150‑1 (MQ=25)

                    gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCAGCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  2:5044165/150‑1 (MQ=25)

                     cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGTACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                 <  2:5434318/150‑1 (MQ=255)

                     cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCggg                                                                                                                                 <  2:3483219/150‑1 (MQ=25)

                     cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                 <  2:5433904/150‑1 (MQ=255)

                       ggCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTa                                                                                                                               <  2:1381406/150‑1 (MQ=255)

                          gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAAc                                                                                                                            <  1:2511719/150‑1 (MQ=25)

                          gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAAc                                                                                                                            <  2:3199961/150‑1 (MQ=25)

                          gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAc                                                                                                                            <  1:5076512/150‑1 (MQ=255)

                          gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAc                                                                                                                            <  2:136685/150‑1 (MQ=255)

                           ggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                           <  2:692916/150‑1 (MQ=255)

                               aTGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCt                                                                                                                       <  2:4468055/150‑1 (MQ=255)

                                    cTGGAGGTGGGGCAGAGCA‑CCCGTGTCGCCCCGGGGTGGGATGGCTGCGAACACAAACGCACTTGATTTCGGGGGGCCGTGTTGGCGCGGGTTTACCGGCCAATGCGCGGGTGAAACCACCGACGAATGCGCGGGTAAACCCGCTCCTAc                                                                                                                  <  2:1418485/150‑1 (MQ=19)

                                                     cATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  >  2:2278559/1‑117 (MQ=18)

                                                     cATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                  <  1:2278559/117‑1 (MQ=18)

                                                          cGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                   >  2:5736697/1‑111 (MQ=18)

                                                          cGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                   <  1:5736697/111‑1 (MQ=18)

                                                                                                                        ggAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCAAACCTCTTGCCAAAATCTTGGAATACCATAATATg                               >  2:5532949/1‑150 (MQ=255)

                                                                                                                         gAGCGGGTTTACCCGCGAAAGCGAGGGTGAAACCACCGACGCAATCGCGGGTAAACCCGCGCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCAAACCTCTTGCCAAAATCTTGGAATACCATAAAATGa                              >  2:347644/1‑150 (MQ=38)

                                                                                                                             gggTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                          >  1:5478382/1‑150 (MQ=25)

                                                                                                                               gTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                        >  2:1469280/1‑150 (MQ=25)

                                                                                                                               gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                              >  2:1062335/1‑128 (MQ=255)

                                                                                                                               gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                              <  1:1062335/128‑1 (MQ=255)

                                                                                                                               gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                        >  1:185095/1‑150 (MQ=255)

                                                                                                                                 ttACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGt                      >  1:4662674/1‑150 (MQ=255)

                                                                                                                                  tccccGCGATTGCGGGGTTTACACCACCGCCGATTTGGGGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa                     <  1:3981661/148‑1 (MQ=255)

                                                                                                                                    cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                   >  2:3844228/1‑150 (MQ=255)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGACTCCCATAATATGATATCCCGTAGAc                  >  2:5557308/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:45375/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:1670124/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:497122/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:4877158/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:5187633/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:1091859/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:4106653/1‑150 (MQ=25)

                                                                                                                                     ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:130222/1‑150 (MQ=255)

                                                                                                                                     ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  2:2545446/1‑150 (MQ=255)

                                                                                                                                     ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:325928/1‑150 (MQ=255)

                                                                                                                                     ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:473799/1‑150 (MQ=255)

                                                                                                                                     ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                  >  1:510846/1‑150 (MQ=255)

                                                                                                                                      cgcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                 >  1:1260374/1‑150 (MQ=255)

                                                                                                                                       gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAATCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCGCCATACCTGTTGCCAAAATCTTGGAATTGCATAATATGTTATCCCGTTGAcgt                >  2:5282493/1‑148 (MQ=21)

                                                                                                                                       gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                              <  2:5204302/120‑1 (MQ=21)

                                                                                                                                       gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                              >  1:5204302/1‑120 (MQ=21)

                                                                                                                                       gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                >  2:898460/1‑150 (MQ=25)

                                                                                                                                       gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                >  1:796132/1‑150 (MQ=25)

                                                                                                                                          aaTGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa             >  2:5305831/1‑150 (MQ=255)

                                                                                                                                            gtcGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt           >  2:222652/3‑150 (MQ=255)

                                                                                                                                             gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTc          >  2:5409460/1‑150 (MQ=255)

                                                                                                                                             gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTc          >  1:4384122/1‑150 (MQ=255)

                                                                                                                                             gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTc          >  1:1886679/1‑150 (MQ=255)

                                                                                                                                              cGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCt         >  1:2332192/1‑150 (MQ=255)

                                                                                                                                               gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTg        >  1:5811009/1‑150 (MQ=255)

                                                                                                                                               gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTg        >  1:4850243/1‑150 (MQ=255)

                                                                                                                                                aGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGc       <  2:1260374/150‑1 (MQ=255)

                                                                                                                                                 gggTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCa      >  2:4587166/1‑150 (MQ=255)

                                                                                                                                                   gTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAgg    >  1:4369481/1‑150 (MQ=255)

                                                                                                                                                     gACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGcc  >  1:2430595/1‑150 (MQ=255)

CCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGT‑TCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCC  >  NC_002947/4339976‑4340273

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