Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,410,937	A→G	6.2%	I333T (ATC→ACC)	PP_1232 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,410,937	0	A	G	6.2%	191.8 / 2.7	65	I333T (ATC→ACC)	PP_1232	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CACCAGCAGGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCTTGGCGG  >  NC_002947/1410801‑1411084                                                                                                                                         |                                                                                                                                                    cACCAGCAGGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGgc                                                                                                                                          <  2:293235/148‑1 (MQ=255)  aCCAGCAGGTCGGCACGCCCCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCGGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGgcgc                                                                                                                                        >  7:363352/1‑149 (MQ=255)    cagcagGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCt                                                                                                                                      <  1:171306/149‑1 (MQ=255)    cagcagGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCt                                                                                                                                      >  6:242481/1‑149 (MQ=255)     agcagGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCtt                                                                                                                                     >  6:201547/1‑149 (MQ=255)          gTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGc                                                                                                                                  >  6:122496/1‑147 (MQ=255)           tCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCcag                                                                                                                               <  5:122496/149‑1 (MQ=255)                acgcacCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCagcagc                                                                                                                           >  4:184513/1‑148 (MQ=255)                acgcacCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCagcagc                                                                                                                           <  3:184513/148‑1 (MQ=255)                acgcacCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCg                                                                                                                          <  4:368673/149‑1 (MQ=255)                    acCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCtt                                                                                                                      >  7:272667/1‑149 (MQ=255)                    acCGGTTTCTACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGGCGGCAAGACGGTTGTTGGAGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATCTCGTTGGGCGCCTTGGCCAGCAGCGGCtt                                                                                                                      <  8:340105/149‑1 (MQ=255)                     ccTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTg                                                                                                                     <  5:201547/149‑1 (MQ=255)                      cTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGa                                                                                                                    <  2:176128/149‑1 (MQ=255)                         tttCCACGGGTAGGTGCCCGCGGACAGCCCCGGCAGCCGTTCCGCGCGCTGCTGGGCGTCGGCAAGCCGGGTGTTCGTGACGGCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCGTGAg                                                                                                                   <  8:132394/147‑1 (MQ=255)                              aCGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCTAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCa                                                                                                            >  2:356481/1‑149 (MQ=255)                                 ggTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGGGCTGCTGGGCGTCGGCAAGACGGTTGTTTGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGc                                                                                                                                  <  7:229110/124‑1 (MQ=255)                                 ggTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGc                                                                                                                                  >  8:229110/1‑124 (MQ=255)                                  gTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACggg                                                                                                        >  2:314805/1‑149 (MQ=255)                                      ttGCCCGGGTTCAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCtt                                                                                                                      <  7:167380/131‑1 (MQ=255)                                      ttGCCCGGGTTCAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCtt                                                                                                                      >  8:167380/1‑131 (MQ=255)                                      ttGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTc                                                                                                    >  5:268983/1‑149 (MQ=255)                                         cccGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGtt                                                                                                 <  8:272667/149‑1 (MQ=255)                                           cGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTc                                                                                                >  2:235677/1‑148 (MQ=255)                                            gggTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGtt                                                                                                 <  5:199462/146‑1 (MQ=255)                                            gggTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCaa                                                                                              <  5:242481/149‑1 (MQ=255)                                             ggTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCATCGTTCAAt                                                                                             >  6:172763/1‑149 (MQ=255)                                              gTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGGTGTCGTTg                                                                                                                                             <  7:339442/100‑1 (MQ=255)                                              gTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGGTGTCGTTg                                                                                                                                             >  8:339442/1‑100 (MQ=255)                                                 cAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTCTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgc                                                                                          <  5:145065/148‑1 (MQ=255)                                                  aGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGGTGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgccc                                                                                        >  2:271039/1‑149 (MQ=255)                                                  aGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgccc                                                                                        >  3:188412/1‑149 (MQ=255)                                                   gCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgcccg                                                                                       >  2:364136/1‑149 (MQ=255)                                                    ccccTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgcccgc                                                                                      >  8:287211/1‑149 (MQ=255)                                                     cccTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTgcccgcc                                                                                     <  8:40309/149‑1 (MQ=255)                                                          cATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCtt                                                                                  <  1:350511/147‑1 (MQ=255)                                                           aTCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGAc                                                                               >  1:51815/1‑149 (MQ=255)                                                               gTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGc                                                                                                                                  <  5:150327/94‑1 (MQ=255)                                                               gTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGc                                                                                                                                  >  6:150327/1‑94 (MQ=255)                                                               catCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCtgt                                                                           <  3:69599/147‑1 (MQ=255)                                                                  cGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGc                                                                         <  4:169059/148‑1 (MQ=255)                                                                              gggCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTTGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAggccaggcca                                                            >  2:117872/1‑149 (MQ=255)                                                                                       cAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTgcg                                                   >  6:215943/1‑149 (MQ=255)                                                                                          gACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTgcggcg                                                <  2:51815/149‑1 (MQ=255)                                                                                                     gTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGTGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTgg                                     <  5:215943/149‑1 (MQ=255)                                                                                                     gTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTgg                                     <  5:172763/149‑1 (MQ=255)                                                                                                      tGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTggg                                    >  1:23304/1‑149 (MQ=255)                                                                                                        aTGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGGTGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGt                                   <  1:271039/148‑1 (MQ=255)                                                                                                          gTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCtt                                >  1:191735/1‑149 (MQ=255)                                                                                                           tCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTc                               <  1:364136/149‑1 (MQ=255)                                                                                                            ccAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTc                               >  6:297400/1‑148 (MQ=255)                                                                                                            ccAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTc                               <  5:297400/148‑1 (MQ=255)                                                                                                            ccAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCg                              >  5:246913/1‑149 (MQ=255)                                                                                                             cAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGt                             >  7:47969/1‑149 (MQ=255)                                                                                                              aGGGCGATCTTGGCCATGTTGTAGGGTCTGGAGGTGTGCGCCTTGGCCAGAAGCTGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTc                            <  8:363352/149‑1 (MQ=255)                                                                                                                 tCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAg                         <  4:188412/149‑1 (MQ=255)                                                                                                                    aTCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAaggc                                                                   <  2:239408/104‑1 (MQ=255)                                                                                                                    aTCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAaggc                                                                   >  1:239408/1‑104 (MQ=255)                                                                                                                    aTCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTg                      >  2:378636/1‑149 (MQ=255)                                                                                                                         tgCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTgg                                     >  4:171485/1‑129 (MQ=255)                                                                                                                         tgCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTgg                                     <  3:171485/129‑1 (MQ=255)                                                                                                                                  gTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAgcgc         <  1:380316/148‑1 (MQ=255)                                                                                                                                  gTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCt        >  6:163421/1‑149 (MQ=255)                                                                                                                                     ggTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCTTgg     >  5:34119/1‑149 (MQ=255)                                                                                                                                        gATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCTTggcgg  >  7:3809/1‑149 (MQ=255)                                                                                                                                         aTGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCTTggcgg  <  7:287211/148‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                    CACCAGCAGGTCGGCACGCACCTGTTTCAACGGGTAGTTGCCCGGGTACAGCCCCTGCATCCGTTCGGCGCGCTGCTGGGCGTCGGCAAGACGGTTGTTGGTGATGTCCAGGTCGATCTGTGCCAAGTTGTAGGTGATGTCGTTGGGCGCCTTGGCCAGCAGCGGCTTGAGCATTTCACGGGCCTCGTTCAATTGCCCGCCCTTGACCTGTGCCAAGGCCAGGCCATAACGTGCGGCGTCCAGCTTTGGGTCTTCGTCCAGTTGGGCGCGGAAGCGCTTGGCGG  >  NC_002947/1410801‑1411084

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