Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	809,807	C→G	6.0%	P124R (CCC→CGC) ‡	PP_0696 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	809,807	0	C	G	6.0%	205.7 / 4.6	68	P124R (CCC→CGC) ‡	PP_0696	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base G (2/2);  total (36/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 = 6.90e-01

ATCGATTTTGAGGTCAGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTGCTGCTCAAGGACATGCT  >  NC_002947/809659‑809954                                                                                                                                                     |                                                                                                                                                    aTCGATTTTGAGGTCAGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGccc                                                                                                                                                    >  1:86326/1‑150 (MQ=255)        ttGAGGTCAGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTc                                                                                                                                             <  1:1900738/150‑1 (MQ=255)              tCAGGCCGGCATCGGTAACCAGCGCACCTTCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGcc                                                                                                                                       >  1:989825/1‑150 (MQ=255)              tCAGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGcc                                                                                                                                       >  1:1680271/1‑150 (MQ=255)                aGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCg                                                                                                                                     >  2:1236948/1‑150 (MQ=255)                      gCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGCGGTTCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCg                                                                                                                               >  1:1384461/1‑150 (MQ=255)                        aTCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGt                                                                                                                             <  2:299164/150‑1 (MQ=255)                         tCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATACCGCCGCAATCAACTGGGGTAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTg                                                                                                                            >  2:543574/1‑150 (MQ=255)                          cGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgc                                                                                                                           >  1:1243499/1‑150 (MQ=255)                           ggTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCCCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgcc                                                                                                                          <  2:721142/150‑1 (MQ=255)                            gTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgccg                                                                                                                         >  2:1390208/1‑150 (MQ=255)                            gTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgccg                                                                                                                         >  2:1322006/1‑150 (MQ=255)                                  aGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACcg                                                                                                                   >  1:509008/1‑150 (MQ=255)                                  aGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACcg                                                                                                                   >  1:1627617/1‑150 (MQ=255)                                     gAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCgcgc                                                                                                                <  1:20158/150‑1 (MQ=255)                                         cgcCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCtttt                                                                                                            >  1:164878/1‑150 (MQ=255)                                          gcCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGCTCCCGCGCCCGGCCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTc                                                                                                           >  1:1707293/1‑150 (MQ=255)                                           ccTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCg                                                                                                          >  1:1096528/1‑150 (MQ=255)                                             tGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCt                                                                                                        <  2:1627617/150‑1 (MQ=255)                                                 aCCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAg                                                                                                    >  1:1248760/1‑150 (MQ=255)                                                   cccccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGtt                                                                                                  >  2:237199/1‑150 (MQ=255)                                                    ccccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgcc                                                                                                                          >  2:30613/1‑125 (MQ=255)                                                    ccccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTgcc                                                                                                                          <  1:30613/125‑1 (MQ=255)                                                     ccccGGCCCGCGCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGGGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgc                                                                                                <  2:333806/150‑1 (MQ=255)                                                     cccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCGGGCTGCTctt                                                                                                                                           <  1:1897395/107‑3 (MQ=255)                                                     cccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCGGGCTGCTctt                                                                                                                                           >  2:1897395/1‑105 (MQ=255)                                                     cccAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgc                                                                                                <  1:1390208/150‑1 (MQ=255)                                                       cAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgctg                                                                                              <  2:164878/150‑1 (MQ=255)                                                        aGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATGCCGCGCCCGGGTAGACGTTCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgctgc                                                                                             <  2:989825/150‑1 (MQ=255)                                                        aGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGATCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgctgc                                                                                             >  1:1454165/1‑150 (MQ=255)                                                           cccGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTg                                                                                          >  1:1037817/1‑150 (MQ=255)                                                               cTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCgg                                                                                      <  2:86326/150‑1 (MQ=255)                                                                          ggCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCctgc                                                                           >  1:12066/1‑150 (MQ=255)                                                                          ggCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCctgc                                                                           <  1:173582/150‑1 (MQ=255)                                                                            ccGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCctgctg                                                                         <  2:1361286/150‑1 (MQ=255)                                                                            ccGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCctgctg                                                                         >  1:1361286/1‑150 (MQ=255)                                                                                 aaTCGAGATCCCGCGCCCGGGCAGTGTGGAAAAAAGCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCtggt                                                                    <  2:1096528/150‑1 (MQ=255)                                                                                   tCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCtggtgg                                                                  >  1:1099457/1‑150 (MQ=255)                                                                                   tCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCtggtgg                                                                  <  2:1243499/150‑1 (MQ=255)                                                                                   tCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCtggtgg                                                                  <  2:1099457/150‑1 (MQ=255)                                                                                          cccGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGc                                                           <  1:237199/150‑1 (MQ=255)                                                                                               gcCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGTCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTg                                                      <  2:1037817/150‑1 (MQ=255)                                                                                               gcCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTg                                                      <  1:191174/150‑1 (MQ=255)                                                                                                    ggCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAg                                                 <  2:1372978/150‑1 (MQ=255)                                                                                                     gCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgctgc                                                                                             <  2:165234/105‑1 (MQ=255)                                                                                                     gCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgctgc                                                                                             >  1:165234/1‑105 (MQ=255)                                                                                                     gCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGc                                                >  1:1020148/1‑150 (MQ=255)                                                                                                          gtgGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGTGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTc                                           <  2:509008/150‑1 (MQ=255)                                                                                                             gCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgct                                                                                               >  2:852822/1‑95 (MQ=255)                                                                                                             gCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTtgct                                                                                               <  1:852822/95‑1 (MQ=255)                                                                                                             gCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCag                                        >  1:1405390/1‑150 (MQ=255)                                                                                                              cAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCaga                                       <  2:508050/150‑1 (MQ=255)                                                                                                                aaaCCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGc                                     <  2:1248760/150‑1 (MQ=255)                                                                                                                    ccgccgTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCcgcg                                 >  2:1675519/1‑150 (MQ=255)                                                                                                                       ccgTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGAcc                              >  1:715805/1‑150 (MQ=255)                                                                                                                       ccgTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTAGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGAcc                              <  1:209783/150‑1 (MQ=255)                                                                                                                         gTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGTAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGAcccc                            >  1:1578729/1‑150 (MQ=255)                                                                                                                            aaGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCctct                         >  1:73094/1‑150 (MQ=255)                                                                                                                                   tCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTc                                                                                                           >  1:1185069/1‑61 (MQ=255)                                                                                                                                   tCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTc                                                                                                           <  2:1185069/61‑1 (MQ=255)                                                                                                                                   tCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACctgc                  <  1:472355/150‑1 (MQ=255)                                                                                                                                   tCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTAGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACctgc                  >  2:218178/1‑150 (MQ=255)                                                                                                                                       ggCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGGCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTgctgct              <  2:1707293/150‑1 (MQ=255)                                                                                                                                          acaAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTGCTGCTCaa           >  2:578363/1‑150 (MQ=255)                                                                                                                                          acaAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTGCTGCTCaa           >  2:1741553/1‑150 (MQ=255)                                                                                                                                           caAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTGCTGCTCAAg          <  1:1236948/150‑1 (MQ=255)                                                                                                                                             ggagCAGCGGGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTgg                                                            <  1:1100776/97‑1 (MQ=255)                                                                                                                                             ggagCAGCGGGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTgg                                                            >  2:1100776/2‑98 (MQ=255)                                                                                                                                                   gcCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCCCCTCCCCGCTGCTCACGGACAtgct  >  2:775485/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATCGATTTTGAGGTCAGGCCGGCATCGGTAACCAGCGAAGCGCCTGTCACCCCCCAGGCCCGCTCCGGCGAGCGGCCGAAAATCGAGATCCCGCGCCCGGGCAGTGTGGCAAAACCCGCCGTCAAGCCTGTCGAGGCACAAGAGCAGCCCGCTGCTCCGCGCCCGGCGCCGGTGCCGCCACCGCGCTTTTCGCTGCAGTTGCTGCGTGCCGGCAGTTGCCTGCTGCTGGTGGAACTGGCCACTGGTCAGCCGTTCCAGAGCCGCGACCCCTCTTACCTGCTGCTCAAGGACATGCT  >  NC_002947/809659‑809954

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