Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,392,039	G→A	5.8%	intergenic (+70/‑45)	aspS → / → PP_1214	aspartate‑‑tRNA ligase/transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,392,039	0	G	A	5.8%	255.9 / 2.7	86	intergenic (+70/‑45)	aspS/PP_1214	aspartate‑‑tRNA ligase/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/44);  new base A (2/3);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTC  >  NC_002947/1391891‑1392185                                                                                                                                                       |                                                                                                                                                      ccAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCg                                                                                                                                                       >  4:16876/1‑149 (MQ=255) ccAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCg                                                                                                                                                       >  5:283503/1‑149 (MQ=255) ccAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCg                                                                                                                                                       >  1:331711/1‑149 (MQ=255)      ccccTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATg                                                                                                                                                  <  2:229458/149‑1 (MQ=255)       cccTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATgg                                                                                                                                                 >  3:328133/1‑149 (MQ=255)            gCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATggggcgg                                                                                                                                            >  7:339944/1‑149 (MQ=255)              aTGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCATGACggggcggga                                                                                                                                          <  1:135790/149‑2 (MQ=255)                 gTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCg                                                                                                                                       >  5:63204/1‑149 (MQ=255)                 gTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCATGACggggcgggacgg                                                                                                                                       <  7:3478/149‑5 (MQ=255)                     aCGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑Ttg                                                                                                                                <  3:16876/149‑1 (MQ=255)                        ccAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGa                                                                                                                             <  5:29646/149‑1 (MQ=255)                         cAAGGCCCTGCGCGAACTGCACATCCGACGTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGAAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGa                                                                                                                             >  1:209726/1‑149 (MQ=255)                         cAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGAt                                                                                                                            >  4:45718/1‑149 (MQ=255)                           aGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTc                                                                                                                          <  2:331711/149‑1 (MQ=255)                           aGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCAACCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTc                                                                                                                          <  1:219924/149‑1 (MQ=255)                             gCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCaa                                                                                                                        <  6:63204/149‑1 (MQ=255)                              cccTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACGAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCTCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCaaa                                                                                                                       <  5:195047/149‑1 (MQ=255)                              cccTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCaaa                                                                                                                       <  2:209726/149‑1 (MQ=255)                              cccTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCaaa                                                                                                                       <  3:347524/149‑1 (MQ=255)                                 tGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGa                                                                                                                    >  5:176284/1‑149 (MQ=255)                                      gAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCTTCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCg                                                                                                                >  2:87087/1‑148 (MQ=255)                                      gAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCgg                                                                                                               >  4:267633/1‑149 (MQ=255)                                              caTCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTa                                                                                                       <  4:328133/149‑1 (MQ=255)                                              caTCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTa                                                                                                       <  8:75295/149‑1 (MQ=255)                                                tCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATg                                                                                                     >  8:266390/1‑149 (MQ=255)                                                   gAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCt                                                                                                  <  3:91919/149‑1 (MQ=255)                                                   gAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCt                                                                                                  <  3:96825/149‑1 (MQ=255)                                                         cgcgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTc                                                                                                           <  3:192601/135‑1 (MQ=255)                                                         cgcgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGt                                                                                               <  6:283503/147‑1 (MQ=255)                                                           cgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTACGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCTGGGCCGG‑‑‑TTCTGATTCAAAGGATTCGGCGACGTAATGGCTGGTCAt                                                                                            >  4:272849/1‑148 (MQ=255)                                                            gAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTc                                                                                                                 >  1:255200/1‑127 (MQ=255)                                                            gAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTc                                                                                                                 <  2:255200/126‑1 (MQ=255)                                                            gAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                                                                          >  6:159834/1‑149 (MQ=255)                                                            gAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                                                                          <  3:267633/149‑1 (MQ=255)                                                             aaCAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTcc                                                                                         <  8:339944/149‑1 (MQ=255)                                                                aGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAg                                                                                      >  2:89094/1‑149 (MQ=255)                                                                    cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTggg                                                                                  >  3:42673/1‑149 (MQ=255)                                                                     aaGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                                                                 <  2:91416/149‑1 (MQ=255)                                                                       ggTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGa                                                                               <  7:266390/149‑1 (MQ=255)                                                                        gTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGaa                                                                              >  8:193055/1‑149 (MQ=255)                                                                        gTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGaa                                                                              <  5:159834/149‑1 (MQ=255)                                                                         tCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAAc                                                                             <  5:342222/149‑1 (MQ=255)                                                                          cGAGTAAGCAGTCCCCTGGCGCATCCAGAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGTGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACa                                                                            <  7:284012/149‑1 (MQ=255)                                                                          cGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACa                                                                            >  8:70388/1‑149 (MQ=255)                                                                           gAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGTCTGGTCATTCCAAGTGGGCGAACa                                                                            >  8:325406/1‑148 (MQ=255)                                                                                aGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAg                                                                       <  5:305236/148‑1 (MQ=255)                                                                                aGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAg                                                                       >  3:46084/1‑148 (MQ=255)                                                                                aGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTGTGCCGTTCGACGCAGGTTATTTCCCGTCCCCCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGCTTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAg                                                                       <  6:272482/148‑1 (MQ=255)                                                                                 gCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCa                                                                     >  7:172159/1‑149 (MQ=255)                                                                                 gCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCGTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCa                                                                     <  2:157165/149‑1 (MQ=255)                                                                                        ccGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCaaa                                                              <  3:45718/149‑1 (MQ=255)                                                                                              gcATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAgcgc                                                        <  8:172159/149‑1 (MQ=255)                                                                                                ttCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGTGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCa                                                      <  7:114/148‑1 (MQ=255)                                                                                                  ccACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGGTGGTGATTGCAAGTGGGCGAACATGAAGCCAGGCAAAGAGCGGCAgg                                                    >  6:223763/1‑149 (MQ=255)                                                                                                  ccACAT‑GGATGCGCATTTGCGTTTCGACGCAGGTTTTTTCCCGTCACGCCCCGTCATGCGGCGGGGCCGT‑‑‑TTGGGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAAAATCACGCACCGCAAAGAGCGACAg                                                     >  8:114/1‑148 (MQ=255)                                                                                                    acaT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTGGGTTGTGATTAAAAGGATTAGGAGTAGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAg                                                     >  1:112111/1‑149 (MQ=255)                                                                                                    acaT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGCGCGCCAGGAt                                                  >  7:224572/1‑149 (MQ=255)                                                                                                       t‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGc                                                >  4:173217/1‑148 (MQ=255)                                                                                                          gATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGcc                                               >  6:241726/1‑147 (MQ=255)                                                                                                            tGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAaga                                           <  6:176284/149‑1 (MQ=255)                                                                                                              cgcCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagag                                          <  1:87087/148‑1 (MQ=255)                                                                                                               gcCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagagag                                        >  5:112726/1‑149 (MQ=255)                                                                                                                     gCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAg                                   <  3:173217/148‑1 (MQ=255)                                                                                                                     gCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAg                                   <  4:46084/148‑1 (MQ=255)                                                                                                                      cGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGt                                 <  7:193055/149‑1 (MQ=255)                                                                                                                       gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGGGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                                         >  3:79081/1‑140 (MQ=255)                                                                                                                       gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGGGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                                         <  4:79081/140‑1 (MQ=255)                                                                                                                       gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCTAGGATGCCAag                                            >  1:243341/1‑137 (MQ=255)                                                                                                                       gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCTAGGATGCCAag                                            <  2:243341/137‑1 (MQ=255)                                                                                                                       gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTc                                >  6:123586/1‑149 (MQ=255)                                                                                                                            gACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCa                            <  7:70388/148‑1 (MQ=255)                                                                                                                               gCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTAt                                                                                                      <  3:186026/71‑1 (MQ=255)                                                                                                                               gCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTAt                                                                                                      >  4:186026/1‑71 (MQ=255)                                                                                                                                  ggTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTg                     <  6:112726/149‑1 (MQ=255)                                                                                                                                   gTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGc                                      >  6:144174/1‑131 (MQ=255)                                                                                                                                   gTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGc                                      <  5:144174/131‑1 (MQ=255)                                                                                                                                    ttttttCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGa                   >  7:274996/1‑149 (MQ=255)                                                                                                                                    ttttttCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGa                   >  7:171803/1‑149 (MQ=255)                                                                                                                                         tcccgtcccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGt              >  6:216894/1‑149 (MQ=255)                                                                                                                                          cccgtcccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGt              >  8:247413/1‑148 (MQ=255)                                                                                                                                           ccgtcccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGa            <  5:216894/149‑1 (MQ=255)                                                                                                                                            ggccccgCCCCATGACGGGGCGGGACCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAaga                                                            <  2:307112/97‑1 (MQ=38)                                                                                                                                            ggccccgCCCCATGACGGGGCGGGACCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAaga                                                            >  1:307112/4‑100 (MQ=38)                                                                                                                                            ggccccgCCCCATGACGGGGCGGGACCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGaa           >  5:272975/4‑149 (MQ=255)                                                                                                                                            cgtcccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGaa           <  5:223763/149‑1 (MQ=255)                                                                                                                                                ccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATc                 <  7:259497/139‑1 (MQ=255)                                                                                                                                                ccgCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATc                 >  8:259497/1‑139 (MQ=255)                                                                                                                                                     ccGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTc  >  7:145671/1‑149 (MQ=255)                                                                                                                                                     ccGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTc  <  1:89094/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGAC‑TGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACAT‑GGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGT‑‑‑TTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTC  >  NC_002947/1391891‑1392185

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