Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,714,509	T→G	10.9%	L104R (CTG→CGG)	PP_1509 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,714,509	0	T	G	10.9%	159.7 / 12.9	64	L104R (CTG→CGG)	PP_1509	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base G (2/5);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CCGCTCCCGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAGGGTG  >  NC_002947/1714362‑1714652                                                                                                                                                    |                                                                                                                                                ccGCTCCCGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTgg                                                                                                                                               <  2:211032/150‑1 (MQ=255)   gCTCCCGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGc                                                                                                                                             <  1:54912/150‑1 (MQ=255)      cccGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccat                                                                                                                                          <  1:1080980/150‑3 (MQ=255)       ccGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGgc                                                                                                                                         >  1:758229/1‑150 (MQ=255)        cGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccatca                                                                                                                                        <  2:29363/150‑5 (MQ=255)        cGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccatca                                                                                                                                        <  2:1095748/150‑5 (MQ=255)        cGGCACCATTGCCCGCTGGGCGCACCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccatca                                                                                                                                        <  2:684681/150‑5 (MQ=255)                    ccGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTa                                                                                                                            >  1:624254/1‑150 (MQ=255)                    ccGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCGGCCTTCGACTTCACCCGCCTGCTGGCCCCTCTCGGCCGGCTCGCATCACCCTa                                                                                                                            >  2:113508/1‑150 (MQ=255)                     cGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTAc                                                                                                                           >  1:1018352/1‑150 (MQ=255)                        tGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTgg                                                                                                                        <  1:963649/150‑1 (MQ=255)                                         cGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccatca                                                                                                                                        <  2:972365/117‑5 (MQ=255)                                         cGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCAGCGGGGccatca                                                                                                                                        >  1:972365/1‑113 (MQ=255)                                                tttCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGc                                                                                                >  1:1421739/1‑150 (MQ=255)                                                     cTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGgc                                                                                           >  1:864715/1‑150 (MQ=255)                                                          cTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGtt                                                                                      >  2:1062561/1‑150 (MQ=255)                                                            gCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTg                                                                                    >  1:841477/1‑150 (MQ=255)                                                             ccAAGTTTCTGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTgg                                                                                   <  2:1187929/150‑1 (MQ=255)                                                             ccAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTgg                                                                                   <  2:1143907/150‑1 (MQ=255)                                                              cAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCGGTGAGAGCAGCTTGAGCCTGCCTTCGACTTCACCCGGATGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTggg                                                                                  >  1:1403648/1‑150 (MQ=255)                                                               aaGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGt                                                                                 <  2:1168996/150‑1 (MQ=255)                                                                 gTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGa                                                                               <  2:454160/150‑1 (MQ=255)                                                                      cgcgGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCttt                                                                          <  2:1414818/150‑1 (MQ=255)                                                                           gATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCTCTTCGACCCGGCGCGTTTGGGTGAGCTTTGTCAc                                                                     >  2:1181868/1‑150 (MQ=255)                                                                             tGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGtt                                                                                                                <  1:963173/105‑1 (MQ=255)                                                                             tGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGtt                                                                                                                >  2:963173/1‑105 (MQ=255)                                                                                   ccATGAATGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCCCCCTACTGGCTGCAGTTGCCGGCCCAGTTCGTCCCCGCGCGTTTGGGTGAGCTTTGTCACTTCCCCGa                                                             >  2:487659/1‑150 (MQ=255)                                                                                    cATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGAt                                                            <  2:758229/150‑1 (MQ=255)                                                                                    cATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGAt                                                            <  1:487659/150‑1 (MQ=255)                                                                                    cATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGACTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGAt                                                            <  1:150071/150‑1 (MQ=255)                                                                                      tGAAGGCGACCTGCGAGACAAGATTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCa                                                          <  2:683714/150‑1 (MQ=255)                                                                                        aaGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGa                                                        <  1:1015172/150‑1 (MQ=255)                                                                                          ggCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCAccc                                                                                                                              <  2:204462/78‑1 (MQ=255)                                                                                          ggCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCAccc                                                                                                                              >  1:204462/1‑78 (MQ=255)                                                                                               ccTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGt                                                 <  2:864715/150‑1 (MQ=255)                                                                                               ccTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGt                                                 <  2:1421739/150‑1 (MQ=255)                                                                                                cTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTg                                                >  2:1033603/1‑150 (MQ=255)                                                                                                 tGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGt                                                                        >  2:975527/1‑125 (MQ=255)                                                                                                 tGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGt                                                                        <  1:975527/125‑1 (MQ=255)                                                                                                 tGCAAGACCAGCTTGAGCCTGCCTTCAACTTCTCCCGGTGGATGGCCCCGCTGGTCCGTCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGc                                               <  2:1403648/150‑1 (MQ=255)                                                                                                   cGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGt                                                                        >  1:520805/1‑123 (MQ=255)                                                                                                   cGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGt                                                                        <  2:520805/123‑1 (MQ=255)                                                                                                   cGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGccc                                             <  1:700584/150‑1 (MQ=255)                                                                                                            cTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCccgtggccgtg                                    <  2:329400/150‑1 (MQ=255)                                                                                                               gAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGaa                                 <  1:1060476/150‑1 (MQ=255)                                                                                                                   ctgcctTCGACTTCACCCGGCTGAGGGCCCCGCTGGGCCGGCGCGTATCACCCTCGTGAGCGCATGGGCCAGGCCATTTCGGCCCGGCTCGTTTGGGTGAGATTTGTCACTTCCTCGATCATATGGCCGTGCCCGTGGCCGTGGAAGTgc                             <  1:1033603/150‑1 (MQ=255)                                                                                                                            aCTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAgg                    >  2:1203620/1‑150 (MQ=255)                                                                                                                            aCTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAAATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAgg                    <  1:113508/150‑1 (MQ=255)                                                                                                                             cTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGc                   >  2:754505/1‑150 (MQ=255)                                                                                                                                  cccGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCtttt              <  2:624254/150‑1 (MQ=255)                                                                                                                                   ccGGCTGATGGCCCCGCTGGGCGGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGGTTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAggacgcttt             >  2:570643/1‑143 (MQ=255)                                                                                                                                   ccGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGAt                                                       <  2:8658/108‑1 (MQ=255)                                                                                                                                   ccGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGAt                                                       >  1:8658/1‑108 (MQ=255)                                                                                                                                    cGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCtttttt            >  2:70750/1‑150 (MQ=255)                                                                                                                                     ggCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTg           <  1:1062561/150‑1 (MQ=255)                                                                                                                                     ggCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTg           <  2:841477/150‑1 (MQ=255)                                                                                                                                          aTGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGCACTGCGCCACCAGGCCTTTTTTGCCAAg      >  2:240289/1‑150 (MQ=255)                                                                                                                                            ggCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAggg    >  2:417816/1‑150 (MQ=255)                                                                                                                                            ggCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAggg    <  1:570643/150‑1 (MQ=255)                                                                                                                                            ggCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGGTGAAGTGCGCAACCAGGCCTTTTTTGCCAAggg    <  1:1181868/150‑1 (MQ=255)                                                                                                                                             gCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCATTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGt                                                                                 <  1:1403826/72‑1 (MQ=255)                                                                                                                                             gCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGt                                                                                 >  2:1403826/1‑72 (MQ=255)                                                                                                                                             gCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAGGGt   >  1:1392612/1‑150 (MQ=255)                                                                                                                                              cccaGCGGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAGGGTg  >  1:359675/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                CCGCTCCCGGCACCATTGCCCGCTGGGCGCAGCTCATGCCCGAGCACTTTCGCTTCACTGCCAAGTTTCCGCGGGATGTCAGCCATGAAGGCGACCTGCGAGACCAGCTTGAGCCTGCCTTCGACTTCACCCGGCTGATGGCCCCGCTGGGCCGGCGCGTATCACCCTACTGGCTGCAGTTGCCGGCCCAGTTCGGCCCGGCGCGTTTGGGTGAGCTTTGTCACTTCCTCGATCAGATTGCCGTGCCCGTGGCCGTGGAAGTGCGCAACCAGGCCTTTTTTGCCAAGGGTG  >  NC_002947/1714362‑1714652

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