Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,001,196	A→T	11.5%	Q57L (CAG→CTG)	PP_0863 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,001,196	0	A	T	11.5%	244.8 / 21.2	87	Q57L (CAG→CTG)	PP_0863	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base T (5/5);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

TGAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001047‑1001345                                                                                                                                                      |                                                                                                                                                       tGAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCca                                                                                                                                                        >  1:2147635/1‑150 (MQ=255) tGAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCca                                                                                                                                                        <  1:2117209/150‑1 (MQ=255)  gAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcag                                                                                                                                                       <  1:2446988/150‑1 (MQ=255)   aGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagc                                                                                                                                                      <  1:164307/150‑1 (MQ=255)    ggAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagca                                                                                                                                                     <  2:1355922/150‑1 (MQ=255)     gAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagcag                                                                                                                                                    >  2:740985/1‑150 (MQ=255)     gAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagcag                                                                                                                                                    >  1:77135/1‑150 (MQ=255)        gtggtgGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGAt                                                                                                                                                 <  1:472664/150‑1 (MQ=255)           gtggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGcc                                                                                                                                              <  2:46158/150‑1 (MQ=255)           gtggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGcc                                                                                                                                              >  1:708116/1‑150 (MQ=255)           gtggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGcc                                                                                                                                              >  1:2248499/1‑150 (MQ=255)            tggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCa                                                                                                                                             <  2:2244852/150‑1 (MQ=255)            tggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCa                                                                                                                                             <  1:2444873/150‑1 (MQ=255)                      cGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCg                                                                                                                                   <  2:534721/150‑1 (MQ=255)                      cGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCGGGATGGGCGTTGCATCCAGCAGGATGCCACGGTTGCCCg                                                                                                                                   >  2:316595/1‑150 (MQ=255)                        gCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGc                                                                                                                                 >  1:232659/1‑150 (MQ=255)                            gCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTg                                                                                                                             <  2:688636/150‑1 (MQ=255)                              aGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGt                                                                                                                           >  2:40185/1‑150 (MQ=255)                               GGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGtt                                                                                                                          >  2:77135‑M1/1‑150 (MQ=255)                                  CATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCgg                                                                                                                       <  1:1689141‑M1/150‑1 (MQ=255)                                  CATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCgg                                                                                                                       <  1:1272901‑M1/150‑1 (MQ=255)                                     gCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTAt                                                                                                                    >  1:195447/1‑150 (MQ=255)                                     gCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTAt                                                                                                                    >  1:655889/1‑150 (MQ=255)                                         gAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                                                >  1:2210514/1‑150 (MQ=255)                                          aGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTg                                                                                                               >  1:852240/1‑150 (MQ=255)                                              aaCTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGcaca                                                                                                           >  1:1329629/1‑150 (MQ=255)                                               aCTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAg                                                                                                          <  2:195447/150‑1 (MQ=255)                                                cTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAgg                                                                                                         >  2:157136/1‑150 (MQ=255)                                                cTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAgg                                                                                                         <  2:2147635/150‑1 (MQ=255)                                                cTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAgg                                                                                                         >  1:1477464/1‑150 (MQ=255)                                                 tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATc                                                                                                                   <  1:361761/139‑1 (MQ=255)                                                 tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATc                                                                                                                   >  2:361761/1‑139 (MQ=255)                                                 tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                                                                        <  2:499840/150‑1 (MQ=255)                                                 tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                                                                        >  1:1042513/1‑150 (MQ=255)                                                 tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                                                                        >  1:2035810/1‑150 (MQ=255)                                                  cccGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCGGGCCTTTGTTCGGTATCGCTGCACAGGGt                                                                                                       <  1:1519538/150‑1 (MQ=255)                                                   ccGNCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTTGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTg                                                                                                      >  2:1519538/1‑150 (MQ=255)                                                   ccGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTg                                                                                                      >  1:1587827/1‑150 (MQ=255)                                                    cGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTgg                                                                                                     >  2:2339606/1‑150 (MQ=255)                                                      gCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCTCGGGACAGGGTGGCa                                                                                                   <  2:1165988/150‑1 (MQ=255)                                                        aaGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGt                                                                                                 >  2:1857590/1‑150 (MQ=255)                                                          ggccGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGc                                                                                               <  1:1889007‑M2/146‑1 (MQ=255)                                                              gccCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                                                           <  1:1711156/150‑1 (MQ=255)                                                              gccCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                                                           <  2:232659/150‑1 (MQ=255)                                                                ccAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagca                                                                                                                                                     >  2:2446452/1‑90 (MQ=255)                                                                ccAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagca                                                                                                                                                     <  1:2446452/90‑1 (MQ=255)                                                                ccAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGgc                                                                                         <  2:852240/150‑1 (MQ=255)                                                                    gCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAg                                                                                                          <  1:1177549/129‑1 (MQ=255)                                                                    gCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAg                                                                                                          >  2:1177549/1‑129 (MQ=255)                                                                             gcggcAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTc                                                                            >  2:759956/1‑150 (MQ=255)                                                                              cggcAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCg                                                                           <  1:1839267/150‑1 (MQ=255)                                                                                gcAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCggg                                                                         <  2:655889/150‑1 (MQ=255)                                                                                         cAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCa                                                                                                   >  2:1114144/1‑115 (MQ=255)                                                                                         cAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCa                                                                                                   <  1:1114144/115‑1 (MQ=255)                                                                                           gggTGTGGTCGGGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑Ga                                                             <  2:882130/150‑1 (MQ=255)                                                                                            ggTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCAAGCATGAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAg                                                            <  1:1843625/150‑1 (MQ=255)                                                                                                tgGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATg                                                        <  2:1689141/150‑1 (MQ=255)                                                                                                  gTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCgg                                                                                                                                  >  2:175564/1‑75 (MQ=255)                                                                                                  gTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCgg                                                                                                                                  <  1:175564/75‑1 (MQ=255)                                                                                                   tCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTgg                                                                                                                            >  2:1246546/1‑80 (MQ=255)                                                                                                   tCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTgg                                                                                                                            <  1:1246546/80‑1 (MQ=255)                                                                                                   gcGAGGCGCAGTTGCTTCTGGGGCAGAGTCTCCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCt                                                     <  2:300965/149‑1 (MQ=255)                                                                                                      aGGCGGTGTTGCTGCTGGGGGAGAGCCTGCTGGATGGGAGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTTGTTCGGTATCGCGGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGACTG‑GAGCATGGCTgtg                                                  <  2:374773/150‑3 (MQ=255)                                                                                                        gcgcAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTgaggg                                                >  2:1125032/1‑149 (MQ=255)                                                                                                              ttGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCg                                                                                                                        >  2:1255272/1‑73 (MQ=255)                                                                                                              ttGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCg                                                                                                                        <  1:1255272/73‑1 (MQ=255)                                                                                                              ttGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAg                                          <  1:1857590/150‑1 (MQ=255)                                                                                                              ttGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAg                                          <  1:581306/150‑1 (MQ=255)                                                                                                               tgctgcTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAgg                                         <  2:1587827/150‑1 (MQ=255)                                                                                                                gctgctGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCCTGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGt                                        <  2:2307496/150‑1 (MQ=255)                                                                                                                gctgctGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGt                                        <  1:1889743/150‑1 (MQ=255)                                                                                                                       gggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGCGTGGCAGTGCCATGGCGCACCACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTg                                 <  2:1329629/150‑1 (MQ=255)                                                                                                                         gCAGATCCTGCTGGATGGGCGGGGCATCCAGCAGGAGGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGTGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGt                               <  2:572239/150‑1 (MQ=255)                                                                                                                                 tgctgGATGGGCGTGGCATCCAGCAGGATGTCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAggcgg                       >  2:1410760/1‑150 (MQ=255)                                                                                                                                 tgctgGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAGTTTGGCGCAggcgg                       >  2:872954/1‑150 (MQ=255)                                                                                                                                 tgctgGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAggcgg                       >  1:288370/1‑150 (MQ=255)                                                                                                                                 tgctgGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAggcgg                       <  1:799807/150‑1 (MQ=255)                                                                                                                                   ctgGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGcc                     >  1:1761596/1‑150 (MQ=255)                                                                                                                                      gATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGtt                  <  2:288370/150‑1 (MQ=255)                                                                                                                                      gATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGGTGGTTCGGGATCGCTGCACAGGGTGGCAGGTCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGtt                  >  1:1993294/1‑150 (MQ=255)                                                                                                                                         gggCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCAc               <  2:1761596/150‑1 (MQ=255)                                                                                                                                           ccGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTa             >  1:821158/2‑150 (MQ=255)                                                                                                                                            cGTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCAc               >  2:1968126/1‑147 (MQ=255)                                                                                                                                            cGTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCAc               <  1:1968126/147‑1 (MQ=255)                                                                                                                                            cGTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTAc            >  1:1396224/1‑150 (MQ=255)                                                                                                                                            cGTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTAc            >  2:783833/1‑150 (MQ=255)                                                                                                                                             gTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGt                                                <  2:1860743/113‑1 (MQ=255)                                                                                                                                             gTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACg           >  1:1294611/1‑150 (MQ=255)                                                                                                                                             gTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACAGGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGt                                                >  1:1860743/1‑113 (MQ=255)                                                                                                                                               ggCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTgg                                                                                                     <  1:1834646/59‑1 (MQ=255)                                                                                                                                               ggCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTgg                                                                                                     >  2:1834646/1‑59 (MQ=255)                                                                                                                                                gCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGAAAGCA‑CGCGGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGccc        >  1:145765/1‑150 (MQ=255)                                                                                                                                                 cATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCg       >  2:1119047/1‑150 (MQ=255)                                                                                                                                                 cATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCg       >  1:1037808/1‑150 (MQ=255)                                                                                                                                                   tCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTg     >  2:1567917/1‑150 (MQ=255)                                                                                                                                                   tCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGATGGGGTGGCGGAGTGTATATGGCGGAGGCGGCCGTGCACTACGCCCGTg     >  2:624260/1‑150 (MQ=255)                                                                                                                                                     cagcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCg   >  1:907837/1‑150 (MQ=255)                                                                                                                                                      agcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCgg  >  1:338832/1‑150 (MQ=255)                                                                                                                                                      agcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCgg  >  1:736357/1‑150 (MQ=255)                                                                                                                                                      agcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCgg  >  1:1793931/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       TGAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTG‑GAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001047‑1001345

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