Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,174,174	T→A	8.2%	L101Q (CTG→CAG)	PP_2786 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,174,174	0	T	A	8.2%	252.7 / 10.6	85	L101Q (CTG→CAG)	PP_2786	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base A (4/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

ATCGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAACGC  >  NC_002947/3174029‑3174323                                                                                                                                                  |                                                                                                                                                      atcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGc                                                                                                                                                   <  2:827949/150‑1 (MQ=255) atcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGc                                                                                                                                                   <  1:1560282/150‑1 (MQ=255)  tcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATACACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGcc                                                                                                                                                  <  2:682759/150‑1 (MQ=255)       aagcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGaa                                                                                                                                             <  2:2023673/150‑1 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           >  2:2429807/1‑150 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           >  2:143541/1‑150 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           >  2:2429624/1‑150 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           <  1:2279098/150‑1 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           >  2:1803675/1‑150 (MQ=255)         gcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                           <  2:1924771/150‑1 (MQ=255)           aagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACa                                                                                                                                         >  1:266098/1‑150 (MQ=255)            agGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAg                                                                                                                                        >  1:505309/1‑150 (MQ=255)             ggTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGc                                                                                                                                       <  1:395113/150‑1 (MQ=255)             ggTAACAGCCGCATCCGCGAGTTAAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGGGACGTCGCCTTTCACACCGCAAACCCGAGCGAGGCCCTGGGCCGCGAAAACAGc                                                                                                                                       >  1:13044/1‑150 (MQ=255)              gTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGcc                                                                                                                                      >  2:156373/1‑150 (MQ=255)                    gCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGa                                                                                                                                <  1:1422644/150‑1 (MQ=255)                     ccGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGAt                                                                                                                               >  2:1875128/1‑150 (MQ=255)                      cGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACa                                                                                                                                         <  1:1796347/139‑1 (MQ=255)                      cGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACa                                                                                                                                         >  2:1796347/1‑139 (MQ=255)                            cgcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTg                                                                                                                        >  2:1336552/1‑150 (MQ=255)                            cgcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTg                                                                                                                        >  2:1339322/1‑150 (MQ=255)                             gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACCGCCAGTGGATGTATTTGc                                                                                                                       >  2:1022216/1‑150 (MQ=255)                             gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGc                                                                                                                       >  1:1894846/1‑150 (MQ=255)                             gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGc                                                                                                                       <  2:1133850/150‑1 (MQ=255)                               gAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCg                                                                                                                     <  1:866207/150‑1 (MQ=255)                               gAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGACAAGTGAAGTGTAATACGTCGCCTTCCATATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCg                                                                                                                     <  1:1485528/150‑1 (MQ=255)                                     aGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGc                                                                                                               >  1:2170609/1‑150 (MQ=255)                                     aGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGc                                                                                                               <  1:1661404/150‑1 (MQ=255)                                          ccTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACa                                                                                                          <  2:280136/150‑1 (MQ=255)                                             gCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGTATGTATTTGCCGGTCAGCCGACAGAc                                                                                                       >  1:604998/1‑150 (MQ=255)                                               aGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGACGACGGCCGCAGGATGTATTTGCCGGTCAGCCGACAGAccc                                                                                                     >  2:1145493/1‑150 (MQ=255)                                               aGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGAccc                                                                                                     >  2:576680/1‑150 (MQ=255)                                                gAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGAc                                                                                                                                          >  1:1846280/1‑112 (MQ=255)                                                gAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGAc                                                                                                                                          <  2:1846280/112‑1 (MQ=255)                                                    gATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGt                                                                                                                                   >  1:1150591/1‑115 (MQ=255)                                                    gATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGt                                                                                                                                   <  2:1150591/115‑1 (MQ=255)                                                         cccGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTc                                                                                           >  1:398822/1‑150 (MQ=255)                                                                aGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGcc                                                                                                              >  1:1893411/1‑124 (MQ=255)                                                                aGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGcc                                                                                                              <  2:1893411/124‑1 (MQ=255)                                                                  cGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGa                                                                                  <  1:1875128/150‑1 (MQ=255)                                                                  cGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGa                                                                                  <  2:604998/150‑1 (MQ=255)                                                                     cAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACaa                                                                               >  2:2030896/1‑150 (MQ=255)                                                                      aaaTACACCATGTATTTCTCCGCGCCAAGGGATGTGCGTAACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACaaa                                                                              <  2:490151/150‑1 (MQ=255)                                                                            aCCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAGCCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGa                                                                                  <  1:2185549/140‑1 (MQ=255)                                                                            aCCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGa                                                                                  >  2:2185549/1‑140 (MQ=255)                                                                            aCCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCggg                                                                        <  2:266098/150‑1 (MQ=255)                                                                               aTGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGTCg                                                                     >  1:1553521/1‑150 (MQ=255)                                                                               aTGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCg                                                                     >  1:267741/1‑150 (MQ=255)                                                                                     ttCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATg                                                               <  2:267741/150‑1 (MQ=255)                                                                                           gcgcCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATg                                                               >  1:1188320/1‑144 (MQ=255)                                                                                           gcgcCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATg                                                               <  2:1188320/144‑1 (MQ=255)                                                                                           gcgcCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGc                                                         <  1:2030896/150‑1 (MQ=255)                                                                                              ccAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGac                                                      >  1:160640/1‑150 (MQ=255)                                                                                              ccAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGac                                                      <  1:576680/150‑1 (MQ=255)                                                                                                aaGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCGTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGacta                                                    >  2:1286943/1‑150 (MQ=255)                                                                                                 aGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGactac                                                   <  1:143541/150‑1 (MQ=255)                                                                                                   tGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTc                                                 <  2:782390/150‑1 (MQ=255)                                                                                                     aTGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCAt                                                                <  2:272852/133‑1 (MQ=255)                                                                                                     aTGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCAt                                                                >  1:272852/1‑133 (MQ=255)                                                                                                       gtgCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTAt                                             <  2:532521/150‑1 (MQ=255)                                                                                                                gTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTc                                    <  2:505309/150‑1 (MQ=255)                                                                                                                 tCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGCCCTCg                                   >  2:39094/1‑150 (MQ=255)                                                                                                                      ttCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAg                              >  2:131011/1‑150 (MQ=255)                                                                                                                      ttCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAg                              >  2:2182240/1‑150 (MQ=255)                                                                                                                          acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATg                                                                                                                              <  2:204030/50‑1 (MQ=255)                                                                                                                          acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATg                                                                                                                              >  1:204030/1‑50 (MQ=255)                                                                                                                          acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTc                          <  1:131011/150‑1 (MQ=255)                                                                                                                            aTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCggg                                                                        >  1:330319/1‑102 (MQ=255)                                                                                                                            aTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCggg                                                                        <  2:330319/102‑1 (MQ=255)                                                                                                                             tCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTcgcg                       >  1:955693/1‑150 (MQ=255)                                                                                                                                    cccGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGa                <  1:2429807/150‑1 (MQ=255)                                                                                                                                    cccGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGa                <  1:2429624/150‑1 (MQ=255)                                                                                                                                    cccGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGa                <  1:1286943/150‑1 (MQ=255)                                                                                                                                     ccGCGCGAGGCCCTGGGCCTAGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGAc               >  1:1210493/1‑150 (MQ=255)                                                                                                                                      tcttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACaa                               >  1:1287339/5‑133 (MQ=255)                                                                                                                                      tcttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACaa                               <  2:1287339/129‑1 (MQ=255)                                                                                                                                        ttcgAGGCCCAGGGCCTCGCAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGa                                  <  2:2449763/126‑1 (MQ=255)                                                                                                                                        ttcgAGGCCCAGGGCCTCGCAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGa                                  >  1:2449763/3‑128 (MQ=255)                                                                                                                                        ttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCCCCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCGGGCCTCGACAGGATTCCCGCCGAGGACtct            >  2:931711/3‑148 (MQ=255)                                                                                                                                        ttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGacta                                                    >  1:182444/3‑110 (MQ=255)                                                                                                                                        ttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGacta                                                    <  2:182444/108‑1 (MQ=255)                                                                                                                                         gcgAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATc           <  2:1894846/150‑1 (MQ=255)                                                                                                                                           gAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGc         <  1:1803675/150‑1 (MQ=255)                                                                                                                                             ggCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACCACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCa       >  1:312273/1‑150 (MQ=255)                                                                                                                                                 cTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAACg   >  2:1276103/1‑150 (MQ=255)                                                                                                                                                  tGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAACGc  >  2:1938313/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                      ATCGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCCGCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAACGC  >  NC_002947/3174029‑3174323

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