BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,751,022	C→G	6.8%	D113E (GAC→GAG)	PP_1551 →	replication protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,751,022	0	C	G	6.8%	322.5 / 12.3	103	D113E (GAC→GAG)	PP_1551	replication protein O
Reads supporting (aligned to +/- strand): ref base C (45/51); new base G (3/4); total (48/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 103 total aligned reads displayed.

GCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGC  >  NC_002947/1750876‑1751168

gCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCg                                                                                                                                                 >  2:922242/1‑150 (MQ=255)

  tCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGc                                                                                                                                               >  2:2060627/1‑150 (MQ=255)

  tCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGAGCCAGGGGTGACATCGGCGTGTATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGc                                                                                                                                               <  1:901283/150‑1 (MQ=255)

    aCCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGGCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCAGACTCGGCTg                                                                                                                                             >  1:2302007/1‑150 (MQ=255)

      ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           <  2:1793988/150‑1 (MQ=255)

      ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           <  2:1185324/150‑1 (MQ=255)

      ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           <  2:1141430/150‑1 (MQ=255)

      ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGAATGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           >  1:1420929/1‑150 (MQ=255)

         gACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgttttg                                                                                                                                        <  1:922242/150‑6 (MQ=255)

         gACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGAATGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTg                                                                                                                                        <  2:1420929/150‑1 (MQ=255)

          acacCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgttttgg                                                                                                                                       <  2:590907/150‑7 (MQ=255)

              cGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCg                                                                                                                                   >  1:646467/1‑150 (MQ=255)

                  tCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATc                                                                                                                               >  2:268856/1‑150 (MQ=255)

                       ggCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTc                                                                                                                          >  2:1628877/1‑150 (MQ=255)

                       ggCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTc                                                                                                                          >  1:2282416/1‑150 (MQ=255)

                        gCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTg                                                                                                                                        >  1:9688/1‑135 (MQ=255)

                        gCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTg                                                                                                                                        <  2:9688/135‑1 (MQ=255)

                            tttCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCagaga                                                                                                                     >  1:524417/1‑149 (MQ=255)

                             ttCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTc                                                                                                                    <  2:2282416/150‑1 (MQ=255)

                              tCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           <  1:2161645/126‑1 (MQ=255)

                              tCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                           >  2:2161645/1‑126 (MQ=255)

                               ccAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTGGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGa                                                                                                                  <  2:248913/150‑1 (MQ=255)

                                cAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCt                                                                                                                                              >  2:2364527/1‑121 (MQ=255)

                                cAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCt                                                                                                                                              <  1:2364527/121‑1 (MQ=255)

                                 aGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGGATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                <  2:685551/150‑1 (MQ=255)

                                 aGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                >  1:787592/1‑150 (MQ=255)

                                 aGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                >  2:1580215/1‑150 (MQ=255)

                                 aGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                >  2:2368255/1‑150 (MQ=255)

                                  gCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCACTTCAGCCGAGTCgg                                                                                                                                                >  2:1059381/1‑117 (MQ=255)

                                  gCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCACTTCAGCCGAGTCgg                                                                                                                                                <  1:1059381/117‑1 (MQ=255)

                                        tgcGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCa                                                                                                         >  2:410159/1‑150 (MQ=255)

                                               ccGTGTCCTGTTCCGCGAAGGTGGTACCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGACTCAGAGTCGAAACAGACCAAAAccgc                                                                                                  >  2:256476/1‑150 (MQ=255)

                                                cGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGACTCGGCTCAGAGTCGAAACAGACCAAAAccgcc                                                                                                 <  1:1687833/150‑1 (MQ=255)

                                                  tgtCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGa                                                                                               <  1:602840/150‑1 (MQ=255)

                                                      cTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTcc                                                                                           >  1:1054976/1‑150 (MQ=255)

                                                      cTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTcc                                                                                           <  1:1702602/150‑1 (MQ=255)

                                                       tGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTccc                                                                                          >  2:1325348/1‑150 (MQ=255)

                                                       tGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTccc                                                                                          <  1:2060627/150‑1 (MQ=255)

                                                         ttCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCctt                                                                                        <  1:306553/150‑1 (MQ=255)

                                                           ccGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCGAAGCTCATGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGAACGAAAAGGCCGACTCCcttct                                                                                      >  1:150635/1‑150 (MQ=255)

                                                           ccGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGAGCCGACACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCcttct                                                                                      <  2:72423/150‑1 (MQ=255)

                                                            cgcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATTGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCcttctt                                                                                     >  1:559097/1‑150 (MQ=255)

                                                             gcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                >  2:2108795/1‑122 (MQ=255)

                                                             gcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                <  1:2108795/122‑1 (MQ=255)

                                                             gcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAg                                                                                                             <  2:1814383/125‑1 (MQ=255)

                                                             gcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAg                                                                                                             >  1:1814383/1‑125 (MQ=255)

                                                             gcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCttt                                                                                    <  1:1628877/150‑1 (MQ=255)

                                                              cgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTa                                                                                   <  1:2317616/150‑1 (MQ=255)

                                                              cgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTa                                                                                   <  2:54716/150‑1 (MQ=255)

                                                              cgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTa                                                                                   >  1:54716/1‑150 (MQ=255)

                                                              cgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTa                                                                                   >  2:2297237/1‑150 (MQ=255)

                                                              cgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTa                                                                                   <  1:2297237/150‑1 (MQ=255)

                                                                    tggtgCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaa                                                                             <  2:524417/150‑1 (MQ=255)

                                                                    tggtgCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaa                                                                             <  2:1358826/150‑1 (MQ=255)

                                                                    tggtgCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaa                                                                             <  2:787592/150‑1 (MQ=255)

                                                                    tggtgCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTCCTTTCTtatac                                                                             <  2:646467/150‑2 (MQ=255)

                                                                        gCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaaa                                                                         <  1:2368255/150‑1 (MQ=255)

                                                                         cccGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaaag                                                                        <  2:1054976/150‑1 (MQ=255)

                                                                         cccGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaaag                                                                        >  2:1819074/1‑150 (MQ=255)

                                                                          ccGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaaaga                                                                       >  1:866590/1‑150 (MQ=255)

                                                                            ggggTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                >  2:1078545/1‑107 (MQ=255)

                                                                            ggggTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaaa                                                                                                                <  1:1078545/107‑1 (MQ=255)

                                                                                 gACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTCCTTTATTCTaagaaa                                                                         <  2:884888/141‑1 (MQ=255)

                                                                                 gACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTCCTTTATTCTaagaaa                                                                         >  1:884888/1‑141 (MQ=255)

                                                                                    aTCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTAt                                                             >  2:616835/1‑150 (MQ=255)

                                                                                    aTCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTAt                                                             <  1:1819074/150‑1 (MQ=255)

                                                                                        gCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCCCAGAGTCCAAACAGACCAAAACCGCCGACTCCCTTCCTTATTCTAAGAAAAAAACCCCCTCTGtca                                                         >  2:1018428/1‑148 (MQ=255)

                                                                                          gtaaaTGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAAt                                                       <  1:410159/147‑1 (MQ=255)

                                                                                           tGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATc                                                      >  2:1055972/1‑150 (MQ=255)

                                                                                              aTGACCCGAAAAACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATcttc                                                   <  1:371635/150‑1 (MQ=255)

                                                                                                  cccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATcttccttc                                               <  1:1580215/150‑1 (MQ=255)

                                                                                                  cccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATcttccttc                                               <  1:1724450/150‑1 (MQ=255)

                                                                                                   ccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgtt                                                                                                                                           >  2:366858/1‑55 (MQ=255)

                                                                                                     gAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATTGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGa                                            <  2:559097/150‑1 (MQ=255)

                                                                                                     gAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATcttcaaaagga                                            >  1:705128/1‑143 (MQ=255)

                                                                                                     gAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGa                                            >  1:2248072/1‑150 (MQ=255)

                                                                                                        ggACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTCAGAATGAAACCCCCTATGTAAATCTTCCTTCGGAAGa                                         <  1:1665091/150‑1 (MQ=255)

                                                                                                         gACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGAt                                        >  2:85558/1‑150 (MQ=255)

                                                                                                         gACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGAt                                        >  2:236213/1‑150 (MQ=255)

                                                                                                            tGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTTGGCTGAAGTGGTCCGGATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTACCTTCTTTATTCTAAGAAAAAAACCCCCTATGTAAATCTTACTTCGCAAGATGtt                                     >  2:2248368/1‑150 (MQ=255)

                                                                                                            tGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGtt                                     <  1:1055972/150‑1 (MQ=255)

                                                                                                             ggATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTa                                    <  1:1646898/150‑1 (MQ=255)

                                                                                                                tttACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACAt                                 <  2:2018310/150‑1 (MQ=255)

                                                                                                                 ttACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTa                                    >  2:306718/1‑146 (MQ=255)

                                                                                                                 ttACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTa                                    <  1:306718/146‑1 (MQ=255)

                                                                                                                     gTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGa                                         >  1:2130867/1‑137 (MQ=255)

                                                                                                                     gTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGa                                         <  2:2130867/137‑1 (MQ=255)

                                                                                                                           ccGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGc                      <  1:616835/150‑1 (MQ=255)

                                                                                                                             gAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCga                    >  1:1982493/1‑150 (MQ=255)

                                                                                                                               aaCAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTa                                    >  2:2076829/1‑132 (MQ=255)

                                                                                                                               aaCAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTa                                    <  1:2076829/132‑1 (MQ=255)

                                                                                                                               aaCAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagc                  <  1:268856/150‑1 (MQ=255)

                                                                                                                                aCAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagcc                 >  2:1432628/1‑150 (MQ=255)

                                                                                                                                  aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCACCCAGCgcgcagg               >  1:2356080/1‑150 (MQ=255)

                                                                                                                                  aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGCAAGATGTTACATGCCCCCCCACCcaaccgg               >  2:423008/1‑150 (MQ=255)

                                                                                                                                     ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccggagc            <  1:236213/150‑1 (MQ=255)

                                                                                                                                     ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAAAGTTCCTTAGGAAGGTGTTACATGCCCCCCCAGCcagccggagc            >  1:380226/1‑150 (MQ=255)

                                                                                                                                      cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccggagcc           <  1:256476/150‑1 (MQ=255)

                                                                                                                                           cAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTc      <  1:423008/150‑1 (MQ=255)

                                                                                                                                               cGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGTCGGCTCCGGc  <  2:866590/150‑1 (MQ=255)

GCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGC  >  NC_002947/1750876‑1751168

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