Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,319	A→T	7.3%	W6R (TGG→AGG) ‡	PP_0621 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	727,319	0	A	T	7.3%	205.2 / 5.2	66	W6R (TGG→AGG) ‡	PP_0621	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCA  >  NC_002947/727181‑727466                                                                                                                                            |                                                                                                                                                    ggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCtt                                                                                                                                          <  2:884148/150‑1 (MQ=255)   ctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttc                                                                                                                                        <  2:248523/150‑1 (MQ=255)    tgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                       <  1:739621/150‑1 (MQ=255)        tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAGAGGTTGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                   >  2:1075101/1‑150 (MQ=255)         ttGTATTTGCGCCACAGCCAGAACGCGGAGGCCTTAGGGGGGAGCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCttttttattg                                                                                                                                  <  2:96862/150‑8 (MQ=255)                tGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGaa                                                                                                                           <  1:102578/150‑1 (MQ=255)                tGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGaa                                                                                                                           <  2:462470/150‑1 (MQ=255)                    ccACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AGAAAAGGGAGGCCCAAGGCCTCCCtt                                                                                                                                          <  1:618786/131‑1 (MQ=255)                    ccACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AGAAAAGGGAGGCCCAAGGCCTCCCtt                                                                                                                                          >  2:618786/1‑131 (MQ=255)                    ccACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTc                                                                                                                       >  2:471572/1‑150 (MQ=255)                     cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATAGTGTCCAGCTTGCCGATTGAAGAAGCCAACAGTCTAGGACGCGACCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                      >  2:229867/1‑150 (MQ=255)                        aGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTAGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                   <  2:512501/150‑1 (MQ=255)                                cgcgGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGAc                                                                                                           <  1:471572/150‑1 (MQ=255)                                    gcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCaaa                                                                                                       <  1:1017129/150‑1 (MQ=255)                                               ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGc                                                                                            >  1:302506/1‑150 (MQ=255)                                                  gATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCACTCAAAAAAGGGAGTCCGAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTg                                                                                          >  1:1051110/1‑150 (MQ=255)                                                       aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAg                                                                                    >  1:484534/1‑150 (MQ=255)                                                       aGAAAAGTAGGCGAACCAGGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAg                                                                                    >  1:347202/1‑150 (MQ=255)                                                          aaaGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGaa                                                                                 <  2:302506/150‑1 (MQ=255)                                                           aaGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAAc                                                                                <  1:208036/150‑1 (MQ=255)                                                              tAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGt                                                                             >  2:422646/1‑150 (MQ=255)                                                                ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCtt                                                                                                                                          <  1:104928/87‑1 (MQ=39)                                                                ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCtt                                                                                                                                          >  2:104928/1‑87 (MQ=39)                                                                   gAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtc                                                                        <  1:422646/150‑1 (MQ=255)                                                                   gAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtc                                                                        <  1:1075101/150‑1 (MQ=255)                                                                    aaCCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtct                                                                       >  2:663406/1‑150 (MQ=255)                                                                     aCCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtctc                                                                      <  2:348785/150‑1 (MQ=255)                                                                       cATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                             >  2:1079643/1‑125 (MQ=255)                                                                       cATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                             <  1:1079643/125‑1 (MQ=255)                                                                              tCCAGCTTGTCGATTGAAGAAGCCAACAGTCTATGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCttt                                                             >  2:298374/1‑150 (MQ=255)                                                                              tCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGGGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCAGCTTCttt                                                             >  1:272501/1‑150 (MQ=255)                                                                               ccAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTc                                                            <  1:663406/150‑1 (MQ=255)                                                                                cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCg                                                           >  2:422334/1‑150 (MQ=255)                                                                                 aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCTTGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  1:1151113/150‑1 (MQ=255)                                                                                 aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          >  1:885272/1‑150 (MQ=255)                                                                                 aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  2:1115018/150‑1 (MQ=255)                                                                                         gATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCg                                                  >  1:522410/1‑150 (MQ=255)                                                                                         gATTGAAGAAGCCAACAAAATCGGACGAGGCGAAT‑AGAAAAGTGGGTCCCAATGCCTCCCTTTCATTCACCGCTGAAGGTCAGTAGAACAGCCTAAACGTCATGTGGCTGAACCAGGAACGGACTGTCTCGGCTTCTTTCGGGCCTTCCg                                                  >  2:564448/1‑150 (MQ=255)                                                                                           ttGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                     <  1:866860/81‑1 (MQ=255)                                                                                           ttGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                     >  2:866860/1‑81 (MQ=255)                                                                                            tGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGc                                                <  2:1051110/150‑1 (MQ=255)                                                                                              aagaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCa                                                                                                  <  1:250065/97‑1 (MQ=255)                                                                                              aagaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCa                                                                                                  >  2:250065/1‑97 (MQ=255)                                                                                                gaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                   <  1:1134308/62‑1 (MQ=255)                                                                                                gaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                   >  2:1134308/1‑62 (MQ=255)                                                                                                 aagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                          <  1:431481/150‑1 (MQ=255)                                                                                                 aagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                          <  2:484534/150‑1 (MQ=255)                                                                                                    ccAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTggg                                       <  2:241777/150‑1 (MQ=255)                                                                                                     cAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTgg                                        <  1:61244/148‑1 (MQ=255)                                                                                                     cAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTgg                                        >  2:61244/1‑148 (MQ=255)                                                                                                               tgACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtctc                            <  1:298374/149‑1 (MQ=255)                                                                                                               ggACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCACTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtctc                            <  1:564448/150‑1 (MQ=255)                                                                                                                   gcgGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCt                        >  2:1176657/1‑150 (MQ=255)                                                                                                                     ggCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCttt                      >  1:961963/1‑150 (MQ=255)                                                                                                                     ggCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCttt                      >  2:572065/1‑150 (MQ=255)                                                                                                                         aaT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                          <  1:952368/126‑1 (MQ=255)                                                                                                                         aaT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                          >  2:952368/1‑126 (MQ=255)                                                                                                                          tgaatgaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgc                  <  2:784171/144‑1 (MQ=255)                                                                                                                             aaaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcagc               <  1:572065/150‑1 (MQ=255)                                                                                                                                   gggAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGATTCTTTCGGGCCTTCCGGCACGCTTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGctct         >  2:654702/1‑150 (MQ=255)                                                                                                                                    ggAGTCCCAACGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGAAAAACGCTCATGCATCTGAACCAGGATCGCTCTGTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCtctc        <  1:84833/150‑1 (MQ=255)                                                                                                                                      aGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCggg                                                         >  2:823635/1‑99 (MQ=255)                                                                                                                                      aGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCggg                                                         <  1:823635/99‑1 (MQ=255)                                                                                                                                        gCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg    >  1:532883/1‑150 (MQ=255)                                                                                                                                        gCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg    >  1:1093629/1‑150 (MQ=255)                                                                                                                                          ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCa  >  2:1195198/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                                    GGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCA  >  NC_002947/727181‑727466

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