Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,753,064	G→C	7.9%	A12P (GCC→CCC) ‡	PP_1555 →	putative Phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,753,064	0	G	C	7.9%	230.1 / 7.2	74	A12P (GCC→CCC) ‡	PP_1555	putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (3/3);  total (39/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 = 7.62e-01

AAGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCC  >  NC_002947/1752922‑1753210                                                                                                                                               |                                                                                                                                                   aaGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGccgccg                                                                                                                                              <  2:236748/149‑1 (MQ=255) aaGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGAAGAGATTCGCCGGCAGGccgccg                                                                                                                                              >  6:331273/1‑149 (MQ=255)  aGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCgg                                                                                                                                             <  7:276858/149‑1 (MQ=255)  aGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGCGCGTACCTGGCCCTGCCAAAAACGGAGGTTCCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCCGCAGgccgcccg                                                                                                                                             >  7:162365/1‑147 (MQ=255)        ccTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCGGCCTGCCGGCGaa                                                                                                                                         <  4:83482/147‑1 (MQ=255)        ccTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCGGCCTGCCGGCGaa                                                                                                                                         <  2:52452/147‑1 (MQ=255)          tGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGc                                                                                                                                     >  5:332102/1‑149 (MQ=255)            gCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCg                                                                                                                                    <  7:198387/148‑1 (MQ=255)            gCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGt                                                                                                                                   <  5:191414/149‑1 (MQ=255)              cGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGa                                                                                                                                 >  4:144728/1‑149 (MQ=255)              cGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGa                                                                                                                                 <  3:300229/149‑1 (MQ=255)                       tCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCgcgc                                                                                                                        >  6:342012/1‑149 (MQ=255)                         aGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGccc                                                                                                                      >  5:145039/1‑149 (MQ=255)                          ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAg                                                                                                                                        >  8:8784/1‑130 (MQ=255)                          ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAg                                                                                                                                        <  7:8784/130‑1 (MQ=255)                          ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCt                                                                                                                     >  6:205780/1‑149 (MQ=255)                           gggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTg                                                                                                                    <  4:25599/149‑1 (MQ=255)                            ggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTgg                                                                                                                   <  6:332102/149‑1 (MQ=255)                              cATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTggg                                                                                                                  >  5:368194/1‑148 (MQ=255)                                tttCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCaa                                                                                                               >  5:172908/1‑149 (MQ=255)                                       gCCATGGCGTGCGTGAGGGCGTACCTGGCCCTGCCAAAATCGGATGTTGCAGCATGACCTACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        <  8:162365/149‑1 (MQ=255)                                       gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGCCCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  2:279436/1‑149 (MQ=255)                                       gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  7:329050/1‑149 (MQ=255)                                       gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  8:317736/1‑149 (MQ=255)                                         cATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCt                                                                                                      >  8:157489/1‑149 (MQ=255)                                           tGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGt                                                                                                    <  7:157489/149‑1 (MQ=255)                                                 gCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTcc                                                                                              >  5:311291/1‑149 (MQ=255)                                                    tGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCt                                                                                           >  3:210694/1‑149 (MQ=255)                                                     gAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGCCCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCtt                                                                                          >  4:18484/1‑149 (MQ=255)                                                              ccTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTc                                                                                 >  1:283721/1‑149 (MQ=255)                                                               cTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCTACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTcc                                                                                >  1:243760/1‑149 (MQ=255)                                                                         cAAAAACGGAGGTTGCAGAATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCGTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTGCATCGCTCCAGcgcgcgcg                                                                      >  8:189057/1‑149 (MQ=255)                                                                          aaaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCgg                                                                     >  4:46552/1‑149 (MQ=255)                                                                           aaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCggg                                                                    >  6:197354/1‑149 (MQ=255)                                                                           aaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCggg                                                                    <  5:197354/149‑1 (MQ=255)                                                                            aaaCGGAGGCTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   <  8:326946/149‑1 (MQ=255)                                                                               cGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGAGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGt                                                                >  7:149911/1‑149 (MQ=255)                                                                                ggAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTg                                                               <  3:144728/149‑1 (MQ=255)                                                                                ggAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTg                                                               <  6:311291/149‑1 (MQ=255)                                                                                 gAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTgg                                                              <  5:205780/149‑1 (MQ=255)                                                                                  aGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGt                                                             <  2:243760/149‑1 (MQ=255)                                                                                   ggTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTa                                                            >  5:70801/1‑149 (MQ=255)                                                                                       gcaggaTGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCt                                                         <  7:317736/148‑1 (MQ=255)                                                                                           caTGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATcc                                                    <  8:329050/149‑1 (MQ=255)                                                                                                 cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              >  8:213993/1‑149 (MQ=255)                                                                                                 cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              <  5:342012/149‑1 (MQ=255)                                                                                                 cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              >  3:357328/1‑149 (MQ=255)                                                                                                      ccgccgCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   <  5:334206/123‑1 (MQ=255)                                                                                                      ccgccgCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   >  6:334206/1‑123 (MQ=255)                                                                                                                 aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  8:149911/148‑1 (MQ=255)                                                                                                                 aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  6:70801/148‑1 (MQ=255)                                                                                                                 aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  6:145039/148‑1 (MQ=255)                                                                                                                   aaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGt                                        <  3:273514/137‑1 (MQ=255)                                                                                                                   aaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGt                                        >  4:273514/1‑137 (MQ=255)                                                                                                                    aTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCgg                           <  3:46552/149‑1 (MQ=255)                                                                                                                        gCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTc                        <  1:166954/148‑1 (MQ=255)                                                                                                                        gCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTc                        >  2:166954/1‑148 (MQ=255)                                                                                                                          gATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATa                     <  3:18484/149‑1 (MQ=255)                                                                                                                           aTGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATTGGCTCATAc                    <  4:210694/149‑1 (MQ=255)                                                                                                                               agagATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                >  3:216515/1‑149 (MQ=255)                                                                                                                                   aTTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  4:70494/1‑149 (MQ=255)                                                                                                                                   aTTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCATGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  2:305642/1‑149 (MQ=255)                                                                                                                                    ttggCCGGCGGCCTGCCGGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  3:87748/4‑148 (MQ=255)                                                                                                                                    ttggCCGGCGGCCTGCCGGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  1:3031/4‑148 (MQ=255)                                                                                                                                     tCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  4:299352/1‑147 (MQ=255)                                                                                                                                      ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                <  2:199244/141‑1 (MQ=255)                                                                                                                                      ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                >  1:199244/2‑142 (MQ=255)                                                                                                                                       gCCGGCATGCCGCCGGCCAAGTTCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCtgagctgagc        <  5:331273/149‑1 (MQ=255)                                                                                                                                       gCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCtgagctgagc        >  8:106608/1‑149 (MQ=255)                                                                                                                                         cGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATAcc                   <  2:293038/136‑1 (MQ=255)                                                                                                                                         cGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATAcc                   >  1:293038/1‑136 (MQ=255)                                                                                                                                          ggCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagctga          >  7:319623/1‑144 (MQ=255)                                                                                                                                          ggCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagctga          <  8:319623/144‑1 (MQ=255)                                                                                                                                            cAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGc   >  1:339756/1‑149 (MQ=255)                                                                                                                                             aGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGcc  <  3:299352/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   AAGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCC  >  NC_002947/1752922‑1753210

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