Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,294,391	G→A	5.1%	G1211E (GGG→GAG)	PP_2893 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,294,391	0	G	A	5.1%	240.7 / 2.0	81	G1211E (GGG→GAG)	PP_2893	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base A (2/2);  total (44/37)
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.51e-01

TGGCAACAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGT  >  NC_002947/3294248‑3294540                                                                                                                                                |                                                                                                                                                      tGGCAACAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTATGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTgg                                                                                                                                                       <  1:1311962/144‑1 (MQ=255) tGGCAACAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTATGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTgg                                                                                                                                                       >  2:1311962/1‑144 (MQ=255)   gCAACAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCaa                                                                                                                                               >  2:1624175/1‑150 (MQ=255)     aaCAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAAGCACAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGc                                                                                                                                             >  1:1556383/1‑150 (MQ=255)     aaCAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCGTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGc                                                                                                                                             >  1:1806863/1‑150 (MQ=255)      aCAGCACCAGCGTGGTAGCCAGCGGTAACATCACCGAGAAGTTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTGACAAGCTGAGCAAGAAGGTCTACCTCGAGTCCGCCAGTGGGCTGGCAAGcg                                                                                                                                            <  2:1724904/150‑2 (MQ=255)           aCCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTg                                                                                                                                       <  2:1495166/150‑1 (MQ=255)           aCCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTg                                                                                                                                       >  1:1495166/1‑150 (MQ=255)             cAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGa                                                                                                                                     >  2:759789/1‑150 (MQ=255)              aGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGAGCTTGCCAGCCCACTggcg                                                                                                                                    <  1:1008913/150‑3 (MQ=255)                     tCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACAtcttct                                                                                                                             <  1:1396591/150‑1 (MQ=255)                     tCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACAtcttct                                                                                                                             <  1:1396952/150‑1 (MQ=255)                      cGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTa                                                                                                                            <  1:759789/150‑1 (MQ=255)                        ccAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCCTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACa                                                                                                                          >  2:501906/1‑150 (MQ=255)                         cAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCCACAAGGTCTACCTCGAGGCCGCCAGTGGGGTGGCAAGCTCACTGGACATCTTCTACaa                                                                                                                         >  1:20085/1‑150 (MQ=255)                         cAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACaa                                                                                                                         <  1:1624175/150‑1 (MQ=255)                           gCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAgc                                                                                                                       <  2:1056897/150‑1 (MQ=255)                                aaCATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGAc                                                                                                                  <  2:950772/150‑1 (MQ=255)                                   aTCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCGCACTGGACATCTTCTACAAGCGCGACtgc                                                                                                               >  2:986834/1‑148 (MQ=255)                                      aCCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCaa                                                                                                                                               <  2:1155722/115‑1 (MQ=255)                                      aCCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCaa                                                                                                                                               >  1:1155722/1‑115 (MQ=255)                                      aCCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAg                                                                                                            >  2:815116/1‑150 (MQ=255)                                        cGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGc                                                                                                          <  1:815116/150‑1 (MQ=255)                                         gagaAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGcc                                                                                                         >  2:522337/1‑150 (MQ=255)                                               cTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGc                                                                                                   <  1:1353039/150‑1 (MQ=255)                                                tCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCg                                                                                                  >  2:470154/1‑150 (MQ=255)                                                  aGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTa                                                                                                                            >  2:1404300/1‑122 (MQ=255)                                                  aGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTa                                                                                                                            <  1:1404300/122‑1 (MQ=255)                                                        cGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGcc                                                                                          <  1:587089/150‑1 (MQ=255)                                                              ggggTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAac                                                                                    <  1:501906/150‑1 (MQ=255)                                                               gggTGGGCGTGTTCGAACCGGCCAACACGATCGGCTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAacg                                                                                   >  1:460520/1‑150 (MQ=255)                                                                ggTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAacgg                                                                                  >  1:1209843/1‑150 (MQ=255)                                                                ggTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAAGCCTACCTCGAGACCGCTGAAGTGCTGTCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAacgg                                                                                  >  1:724379/1‑150 (MQ=255)                                                                 gTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCa                                                                                         >  1:352306/1‑142 (MQ=255)                                                                 gTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCa                                                                                         <  2:352306/142‑1 (MQ=255)                                                                   gggCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACCGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAacggacg                                                                               <  2:22798/150‑1 (MQ=255)                                                                    ggCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGGAGACAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCACAGTGGCGGCCATCAacggacgg                                                                              >  1:612289/1‑150 (MQ=255)                                                                     gCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGt                                                                             >  1:1529383/1‑150 (MQ=255)                                                                          ttCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGGGgccg                                                                        >  1:585310/1‑150 (MQ=255)                                                                                 cGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGAt                                                                 <  1:1718132/150‑1 (MQ=255)                                                                                   gCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGa                                                               >  2:1552924/1‑150 (MQ=255)                                                                                    ccAACACGCTCGCGTTGCGGTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGaa                                                              <  2:722899/150‑1 (MQ=255)                                                                                     cAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAAc                                                             >  2:305564/1‑150 (MQ=255)                                                                                      aaCACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACt                                                            >  1:331832/1‑150 (MQ=255)                                                                                         acGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCTGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGt                                                         <  2:947527/150‑1 (MQ=255)                                                                                         acGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCAGCCTCGATGAACTCGt                                                         <  2:1529383/150‑1 (MQ=255)                                                                                                  ttgcgttACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCgg                                                <  2:1209843/150‑1 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAggg                                         <  1:1835404/148‑1 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAggg                                         >  2:1835404/1‑148 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCg                                       <  2:724379/150‑1 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCg                                       >  2:791529/1‑150 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCg                                       <  2:1872798/150‑1 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCg                                       >  2:810875/1‑150 (MQ=255)                                                                                                           aaGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCg                                       >  1:1461154/1‑150 (MQ=255)                                                                                                             gCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCggt                                     <  1:554722/150‑1 (MQ=255)                                                                                                             gCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCggt                                     >  1:489210/1‑150 (MQ=255)                                                                                                             gCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCggt                                     >  2:142613/1‑150 (MQ=255)                                                                                                              cTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCggtg                                    >  2:668629/1‑150 (MQ=255)                                                                                                                    aagaagGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCggtggggtgg                              >  1:1106499/1‑150 (MQ=255)                                                                                                                    aagaagGTCTACCTCGAGGCCGCCAGTGAGCTTGCCAGCCCACTgg                                                                                                                                      >  1:1772319/1‑46 (MQ=9)                                                                                                                    aagaagGTCTACCTCGAGGCCGCCAGTGAGCTTGCCAGCCCACTgg                                                                                                                                      <  2:1772319/46‑1 (MQ=9)                                                                                                                     agaagGTCTACCTCGAGGCCGCCCGTGCGCTGCCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGggtggt                             <  2:1556383/150‑1 (MQ=255)                                                                                                                          gTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGTACGGTTGCCGCCTCGATGAACTCGATGGTGTCGTGAGTGGGCGGTGGGGaggtggtac                        >  2:398680/1‑147 (MQ=255)                                                                                                                             tACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTa                     >  2:309771/1‑150 (MQ=255)                                                                                                                               ccTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTAcc                   <  1:791529/150‑1 (MQ=255)                                                                                                                                cTCGAGGCCGCCAGTGGGCTGGCACGGTCCAGGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCg                  >  1:1506949/1‑150 (MQ=255)                                                                                                                                cTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACAtctt                                                                                                                               >  1:1435113/1‑41 (MQ=255)                                                                                                                                cTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACAtctt                                                                                                                               <  2:1435113/41‑1 (MQ=255)                                                                                                                                cTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCg                  >  2:642785/1‑150 (MQ=255)                                                                                                                                cTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCAAGCAGCGTGGCGGCGATCAACGGACGCTGGACGCCTCGATGAACTGGTGTGGGTCTGGTGAGGGCGGTGGGGTGGTGGCGAGtccgc                  >  2:1288388/1‑146 (MQ=255)                                                                                                                                 tCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCgg                 >  1:825018/1‑150 (MQ=255)                                                                                                                                   gAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGa               <  1:1598874/150‑1 (MQ=255)                                                                                                                                   gAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGa               <  2:585310/150‑1 (MQ=255)                                                                                                                                    aGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGAc              <  1:1552924/150‑1 (MQ=255)                                                                                                                                       ccgccAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGTCCGCCTCGATGAACTCGTTGGTGTCGGGTGAg                                           <  2:1658532/118‑1 (MQ=255)                                                                                                                                       ccgccAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGTCCGCCTCGATGAACTCGTTGGTGTCGGGTGAg                                           >  1:1658532/1‑118 (MQ=255)                                                                                                                                         gccAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTg                                             >  2:1649003/1‑114 (MQ=255)                                                                                                                                         gccAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTg                                             <  1:1649003/114‑1 (MQ=255)                                                                                                                                          ccAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACCCGTTGGGGTCGGGTGCGGGAGGTGGGGTGGCGGCGCGTACCGGGACGATCAg        >  2:1896344/1‑150 (MQ=255)                                                                                                                                          ccAGTGAGCTTGCCAGCCCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAg        >  1:1352200/1‑150 (MQ=255)                                                                                                                                           cAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAgg       >  1:1528849/1‑150 (MQ=255)                                                                                                                                            aGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGt      <  1:142613/150‑1 (MQ=255)                                                                                                                                             gTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGtt     <  1:470154/150‑1 (MQ=255)                                                                                                                                               gggCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTgg   <  1:305564/150‑1 (MQ=255)                                                                                                                                                ggCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTCGTGGCGCGTACCGGGACGATCAGGTTGgt  <  2:1528849/150‑1 (MQ=255)                                                                                                                                                ggCTGGCAAGCTCACTGGACATATTCTACAAGCGCGACTTCTAGCCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGGCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGgt  <  2:222483/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      TGGCAACAGCACCAGCGTGGTCGCCAGCGGTAACATCACCGAGAAGCTCAGCCTTCGTTACGGGGTGGGCGTGTTCGAACCGGCCAACACGATCGCGTTGCGTTACAAGCTGAGCAAGAAGGTCTACCTCGAGGCCGCCAGTGGGCTGGCAAGCTCACTGGACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGT  >  NC_002947/3294248‑3294540

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