Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,120	C→G	7.6%	intergenic (+40/+61)	PP_1661 → / ← PP_1662	dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,859,120	0	C	G	7.6%	194.8 / 3.7	62	intergenic (+40/+61)	PP_1661/PP_1662	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACC  >  NC_002947/1858976‑1859265                                                                                                                                                  |                                                                                                                                                  gAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCTgcggca                                                                                                                                               >  2:176004/1‑149 (MQ=255)      gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTg                                                                                                                                          >  2:354470/1‑149 (MQ=255)        aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCTGCGGCAttgccgg                                                                                                                                        <  1:176004/149‑6 (MQ=255)        aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTggg                                                                                                                                        <  7:171167/149‑1 (MQ=255)         gCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgggg                                                                                                                                       <  1:354470/149‑1 (MQ=255)             tgCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGtt                                                                                                                                     <  7:190658/147‑1 (MQ=255)               cATCAATGACGAACTGGACGGGCGTTACGGAGCCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGTTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGc                                                                                                                                   >  1:33314/1‑148 (MQ=255)                aTCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                                >  2:134633/1‑149 (MQ=255)                  cAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACt                                                                                                                              >  8:46320/1‑149 (MQ=255)                              ggACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                                    >  4:261148/1‑147 (MQ=255)                               gACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCTGTAGAGGTGCGTTTCTGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGa                                                                                                                  <  2:33314/148‑1 (MQ=255)                                   ggCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                             <  3:261148/149‑1 (MQ=255)                                     cGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                           <  4:326568/149‑1 (MQ=255)                                     cGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                           <  8:185904/149‑1 (MQ=255)                                     cGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                           <  3:177453/149‑1 (MQ=255)                                      gTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTg                                                                                                                                          <  8:129747/117‑1 (MQ=255)                                      gTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTg                                                                                                                                          >  7:129747/1‑117 (MQ=255)                                      gTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgggg                                                                                                          >  2:318290/1‑149 (MQ=255)                                          cGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTc                                                                                                       <  7:304376/148‑1 (MQ=255)                                           gg‑ggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgc                                                                                                     <  1:332229/149‑1 (MQ=255)                                              ggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgccg                                                                                                   >  5:265128/1‑149 (MQ=255)                                                   gCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCaa                                                                                                <  1:134633/147‑1 (MQ=255)                                                   gCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATg                                                                                              >  1:62492/1‑149 (MQ=255)                                                   gCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATg                                                                                              >  4:32393/1‑149 (MQ=255)                                                        aGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAg                                                                                         <  8:55098/149‑1 (MQ=255)                                                         gCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                        <  4:348346/149‑1 (MQ=255)                                                           gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATg                                                                                              <  5:268725/141‑1 (MQ=255)                                                           gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCAGCAATg                                                                                              >  6:268725/1‑141 (MQ=255)                                                            aTAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGTTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGt                                                                                      >  3:349488/1‑148 (MQ=255)                                                               aCCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAtt                                                                                  <  3:32393/149‑1 (MQ=255)                                                               aCCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAt                                                                                   <  7:46320/148‑1 (MQ=255)                                                                 cTGGGCAGCCTGTCGGTAGAGGTGCGGTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgt                                                                                <  4:155562/149‑1 (MQ=255)                                                                   gggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgt                                                                              >  7:39981/1‑149 (MQ=255)                                                                       aGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGgcg                                                                          >  7:204752/1‑149 (MQ=255)                                                                                gTAGAGGTGCGTTTCGGGGCGGTGTAGGTCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                                 >  8:286765/1‑149 (MQ=255)                                                                                      gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAgg                                                           >  3:259362/1‑149 (MQ=255)                                                                                       tGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAg                                                            <  3:59540/147‑1 (MQ=255)                                                                                       tGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAg                                                            >  4:59540/1‑147 (MQ=255)                                                                                          gTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCt                                                        >  2:277363/1‑148 (MQ=255)                                                                                               ggggCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGc                                                  >  4:289974/1‑149 (MQ=255)                                                                                                gggCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCt                                                 <  6:308243/149‑1 (MQ=255)                                                                                                          ggCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTa                                       <  8:39981/149‑1 (MQ=255)                                                                                                            cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTAtt                                     >  1:159683/1‑149 (MQ=255)                                                                                                             ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                    >  3:171470/1‑149 (MQ=255)                                                                                                             ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGGCTTCAAGCTCAAGGCGGTATTg                                    >  3:55650/1‑149 (MQ=255)                                                                                                                  cATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGtgtg                               >  2:103446/1‑149 (MQ=255)                                                                                                                         aCCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtg                        >  8:232521/1‑149 (MQ=255)                                                                                                                          ccccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgt                                  <  4:247303/138‑1 (MQ=255)                                                                                                                          ccccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgt                                  >  3:247303/1‑138 (MQ=255)                                                                                                                           cccaaaaTGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtgt                       >  2:174147/7‑148 (MQ=38)                                                                                                                           cccaaaaTGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtgt                       <  1:174147/142‑1 (MQ=38)                                                                                                                           cccaaaaTGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtg                        <  6:57902/141‑1 (MQ=38)                                                                                                                           cccaaaaTGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtg                        >  5:57902/7‑147 (MQ=38)                                                                                                                            ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                     >  2:31265/1‑149 (MQ=255)                                                                                                                             cACGATTCTCCGGCAATGCCGGAGCATTTTGGGGTTGCTCTCTGAGCCGAGGCGAGTCAGGGGTTCGCCGGAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg                    <  8:204752/149‑1 (MQ=255)                                                                                                                             cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg                    >  6:177296/1‑149 (MQ=255)                                                                                                                                 aTGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAt                <  2:62492/149‑1 (MQ=255)                                                                                                                                  tGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAAGGCGGTATTGGTGTGCTTCGTGTTCGTTATg               <  1:31265/149‑1 (MQ=255)                                                                                                                                  tGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAt                <  4:259362/148‑1 (MQ=255)                                                                                                                                  tGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATg               <  3:289974/149‑1 (MQ=255)                                                                                                                                     tCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGAGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGAGCGCTGGTTTTCCAGGTATTCAAGCTCAAAGCGGTATTGGTGTGCTTCGTGTTCGTTATGGaa            >  7:62193/1‑149 (MQ=255)                                                                                                                                              gCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCccac   >  1:136270/1‑149 (MQ=255)                                                                                                                                               ccGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCccacc  <  2:47733/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGG‑GGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACC  >  NC_002947/1858976‑1859265

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