Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,625,587	Δ1 bp	5.4%	intergenic (‑57/+20)	PP_3195 ← / ← PP_3196	hypothetical protein/nodulation protein d1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,625,587	0	T	.	5.4%	226.4 / 6.9	74	intergenic (‑57/+20)	PP_3195/PP_3196	hypothetical protein/nodulation protein d1
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base . (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAATTCACGCATA  >  NC_002947/3625442‑3625736                                                                                                                                                   |                                                                                                                                                      cGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTgg                                                                                                                                                   <  1:1486080/150‑1 (MQ=255)  ggTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGAATGCGAGGTATAGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTggg                                                                                                                                                  <  1:2233638/150‑1 (MQ=255)   gTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGTCGACAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCAAG‑CACCGCTGCGGTGTCTgggg                                                                                                                                                 >  1:760923/1‑150 (MQ=255)    tCGTGGAGTCCGACGCCCACTTGAGCGACAGGACGACCAGGCGAGGTCCCAAGGCAAGGTACACAGGCAGGAGAGGGTCGATTTAAGCTTGGTCGGTGAGGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGc                                                                                                                                                >  2:882269/1‑150 (MQ=255)         gAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTggg                                                                                                                                                  >  1:997348/1‑143 (MQ=255)         gAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTggg                                                                                                                                                  <  2:997348/143‑1 (MQ=255)              cGGCGCCCACTTCAGCAACAGGCAGGAAGGGTGCCGTTTCACCGGATGGTACGCGGGCAGGGCAGGGTTGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgccgcc                                                                                                                                      <  2:1130150/150‑1 (MQ=255)                   cccACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATAGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCa                                                                                                                                 >  2:720666/1‑150 (MQ=255)                     cACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAgg                                                                                                                               >  2:2126044/1‑150 (MQ=255)                     cACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAgg                                                                                                                               >  1:1751227/1‑150 (MQ=255)                      aCTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGc                                                                                                                              >  2:1575259/1‑150 (MQ=255)                        ttGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCgg                                                                                                                            >  2:1017866/1‑150 (MQ=255)                           aGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTgg                                                                                                                         >  2:1486028/1‑150 (MQ=255)                           aGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTgg                                                                                                                         >  2:1486579/1‑150 (MQ=255)                            gCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgc                                                                                                                        >  1:2334324/1‑150 (MQ=255)                             cGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATg                                                                                                                                            >  1:269857/1‑129 (MQ=255)                             cGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATg                                                                                                                                            <  2:269857/129‑1 (MQ=255)                                caggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGt                                                                                                                           <  2:1147981/143‑1 (MQ=255)                                caggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGt                                                                                                                           >  1:1147981/1‑143 (MQ=255)                                 aggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                                   <  1:1701971/150‑1 (MQ=255)                                 aggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                                   <  2:582431/150‑1 (MQ=255)                                 aggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                                   <  2:1994240/150‑1 (MQ=255)                                  ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTgg                                                                                                                                                   <  1:1101148/117‑1 (MQ=255)                                  ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTgg                                                                                                                                                   >  2:1101148/1‑117 (MQ=255)                                   gcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGc                                                                                                                              >  1:178796/1‑137 (MQ=255)                                   gcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGc                                                                                                                              <  2:178796/137‑1 (MQ=255)                                   gcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCaga                                                                                                                 <  1:1958270/150‑1 (MQ=255)                                   gcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCaga                                                                                                                 <  1:836615/150‑1 (MQ=255)                                      ggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCa                                                                                                              >  1:1914397/1‑150 (MQ=255)                                                 gTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGt                                                                                                   <  2:438303/150‑1 (MQ=255)                                                       aGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGgc                                                                                             >  1:399665/1‑150 (MQ=255)                                                        ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGGTG‑CTGGGGcatc                                                                                                                                             <  2:1543559/101‑4 (MQ=37)                                                        ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGGTG‑CTGGGGcatc                                                                                                                                             >  1:1543559/1‑98 (MQ=37)                                                         gCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                                   <  1:1096397/126‑1 (MQ=255)                                                         gCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                                   >  2:1096397/1‑126 (MQ=255)                                                           aaGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCAc                                                                                         >  2:1177749/1‑150 (MQ=255)                                                            aGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                                        >  2:123518/1‑150 (MQ=255)                                                            aGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                                        <  1:1017866/150‑1 (MQ=255)                                                              gTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGa                                                                                      >  1:1743280/1‑150 (MQ=255)                                                               tACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGaa                                                                                     <  1:2126044/150‑1 (MQ=255)                                                                     ggcagggcagGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCACGCAGAGCATCTGGCTGCGt                                                                                                   >  2:430004/1‑130 (MQ=255)                                                                     ggcagggcagGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCACGCAGAGCATCTGGCTGCGt                                                                                                   <  1:430004/130‑1 (MQ=255)                                                                        agggcagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                            >  2:2271970/1‑150 (MQ=255)                                                                            cagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGAt                                                                        <  2:1743280/150‑1 (MQ=255)                                                                              gggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGCTGCGCCAT‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATgc                                                                      <  2:1271616/150‑1 (MQ=255)                                                                                  ggggCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTc                                                                  <  1:1575259/150‑1 (MQ=255)                                                                                  ggggCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTc                                                                  >  1:77346/1‑150 (MQ=255)                                                                                    ggCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                            >  2:708388/1‑138 (MQ=255)                                                                                    ggCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                            <  1:708388/138‑1 (MQ=255)                                                                                    ggCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTcgc                                                                >  1:411509/1‑150 (MQ=255)                                                                                      ccAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGGGCATCTGTCTGCGTGCTGGCGCGCTGAATGACATCGCTGATGCGCTCGCGa                                                              >  1:757660/1‑150 (MQ=255)                                                                                       cAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGaa                                                             <  2:760923/150‑1 (MQ=255)                                                                                          tCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCCG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGcc                                                          >  1:1023878/1‑150 (MQ=255)                                                                                          tCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATTGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTg                                                                                               >  2:1499699/1‑113 (MQ=255)                                                                                          tCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGcc                                                          >  2:98006/1‑150 (MQ=255)                                                                                          tCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCGGGGGCGATTGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTg                                                                                               <  1:1499699/113‑1 (MQ=255)                                                                                                tCAGGGATGCGAGGTATCGCGAACCGCGATGCCCCAG‑CACCGCTTCTGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCtga                                                    <  1:1892638/150‑2 (MQ=255)                                                                                                   gggATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCa                                                 >  2:2232689/1‑150 (MQ=255)                                                                                                    ggATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAACCATCTGTGCAt                                                >  2:108050/1‑150 (MQ=255)                                                                                                       tGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATggg                                             <  2:234505/150‑1 (MQ=255)                                                                                                       tGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACAGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATggg                                             <  2:1751227/150‑1 (MQ=255)                                                                                                          gAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTc                                          <  1:2232689/150‑1 (MQ=255)                                                                                                               aTCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGt                                     >  1:2255228/1‑150 (MQ=255)                                                                                                                                cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGcgagc                                    <  2:2196595/134‑5 (MQ=255)                                                                                                                                cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGcgagc                                    >  1:2196595/1‑130 (MQ=255)                                                                                                                                cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcac                    >  1:143783/1‑150 (MQ=255)                                                                                                                                 ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCTGATCAGGGGGCGGCCAGATCGGGATCTACCTGCACGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcaca                   >  1:1252536/1‑150 (MQ=255)                                                                                                                                 ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGAGATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcaca                   >  1:1601185/1‑150 (MQ=255)                                                                                                                                  cAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGCGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATCCGCTCGCGAAGCCATCTGTGCATGCGGTCGCCGTCCATGCTGCGGTGCCCCAg                  >  2:1187969/1‑150 (MQ=255)                                                                                                                                  cAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAg                  <  1:2271970/150‑1 (MQ=255)                                                                                                                                       aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtgctg                             <  2:2299002/135‑1 (MQ=255)                                                                                                                                       aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtgctg                             >  1:2299002/1‑135 (MQ=255)                                                                                                                                          gcagCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt               >  2:231371/4‑145 (MQ=255)                                                                                                                                          gcagCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt               <  1:231371/142‑1 (MQ=255)                                                                                                                                                   tCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAATTCACGCATa  >  2:1143580/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAATTCACGCATA  >  NC_002947/3625442‑3625736

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