Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,175,196	C→T	10.1%	L91F (CTC→TTC)	PP_1030 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,175,196	0	C	T	10.1%	159.8 / 7.8	59	L91F (CTC→TTC)	PP_1030	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/20);  new base T (3/3);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTGC  >  NC_002947/1175053‑1175343                                                                                                                                                |                                                                                                                                                    gCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCgag                                                                                                                                                >  2:221752/1‑148 (MQ=255) gCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCgag                                                                                                                                                >  2:225817/1‑148 (MQ=255) gCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCgaa                                                                                                                                                <  3:66238/149‑1 (MQ=255) gCACCGGTAACCGCATGGGCCATGCCCTGGCGGAAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCgag                                                                                                                                                >  2:221755/1‑148 (MQ=255)   aCCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCCCCGAAGGCTTCCTCgaag                                                                                                                                               >  2:300718/1‑148 (MQ=255)          aCCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGt                                                                                                                                       >  7:124125/1‑149 (MQ=255)           ccGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCATACGCTGCATCTATACCCCCACCGAAGGCTTCTTCGAGGAAGccttcg                                                                                                                                      <  1:225817/149‑5 (MQ=255)           ccGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCGAGGAAGccttcg                                                                                                                                      <  1:221752/149‑5 (MQ=255)           ccGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCGAGGAAGccttcg                                                                                                                                      <  1:221755/149‑5 (MQ=255)               aTGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGccc                                                                                                                                  >  2:281708/1‑149 (MQ=255)                  ggCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGa                                                                                                                               <  5:16065/149‑1 (MQ=255)                  ggCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGa                                                                                                                               <  8:124125/149‑1 (MQ=255)                   gCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAg                                                                                                                              <  6:304226/149‑1 (MQ=255)                    ccATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGt                                                                                                                             >  5:241631/1‑149 (MQ=255)                     cATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTg                                                                                                                            >  2:213589/1‑149 (MQ=255)                                 gcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCg                                                                                                                                 >  4:302849/1‑132 (MQ=255)                                 gcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCg                                                                                                                                 <  3:302849/132‑1 (MQ=255)                                      ggCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGaggtagt                                                                                                             >  3:138148/1‑146 (MQ=255)                                           gcAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGTAGTGCCGGCGCGAGGTAGGTGAGCCg                                                                                                      <  1:108455/149‑1 (MQ=255)                                               aGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCgttgt                                                                                                  <  6:241631/149‑1 (MQ=255)                                                     ccTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGa                                                                                              >  4:211185/1‑147 (MQ=255)                                                     ccTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGa                                                                                              >  3:84411/1‑147 (MQ=255)                                                         cACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCt                                                                                        <  2:177881/149‑1 (MQ=255)                                                          aCCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTg                                                                                       >  1:123232/1‑149 (MQ=255)                                                           cccGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTgg                                                                                      <  1:281708/149‑1 (MQ=255)                                                            ccGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTgg                                                                                      >  4:142010/1‑148 (MQ=255)                                                            ccGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACACCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTgg                                                                                      >  6:217454/1‑148 (MQ=255)                                                              ggAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTgg                                                                                      >  3:265373/1‑146 (MQ=255)                                                                     cTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGAt                                                                            <  5:217454/149‑1 (MQ=255)                                                                      tGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATg                                                                           >  6:164756/1‑149 (MQ=255)                                                                         ccTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAACGCCTGGAAGAGTGGATGCCg                                                                        >  4:244689/1‑149 (MQ=255)                                                                          cTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCgg                                                                       >  1:241452/1‑149 (MQ=255)                                                                           ttCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCgg                                                                       >  4:301858/1‑148 (MQ=255)                                                                                cGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCTACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGa                                                                 <  5:66007/149‑1 (MQ=255)                                                                                 gAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAg                                                                <  7:9932/149‑1 (MQ=255)                                                                                       gTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGa                                                                                     >  3:160719/1‑122 (MQ=255)                                                                                       gTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGa                                                                                     <  4:160719/122‑1 (MQ=255)                                                                                          aaCGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACtt                                                       >  5:254703/1‑149 (MQ=255)                                                                                               ccTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCt                                                  >  4:190718/1‑149 (MQ=255)                                                                                               ccTGGAAGTGATGCTCAAACGCAGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCt                                                  >  7:193532/1‑149 (MQ=255)                                                                                                cTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTg                                                 >  7:95477/1‑149 (MQ=255)                                                                                                    aaGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCa                                             <  3:211185/149‑1 (MQ=255)                                                                                                         aTGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACtgtg                                        >  4:88193/1‑149 (MQ=255)                                                                                                          tGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGc                                       <  8:95477/149‑1 (MQ=255)                                                                                                             taaaaCGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTAc                                                    <  8:246038/131‑1 (MQ=255)                                                                                                             tCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGc                                    >  4:282909/1‑149 (MQ=255)                                                                                                             tCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGc                                    >  8:166504/1‑149 (MQ=255)                                                                                                               aaaCGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCtt                                  <  1:300718/149‑1 (MQ=255)                                                                                                                aaCGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTc                                 >  7:25714/1‑149 (MQ=255)                                                                                                                    cTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAACCCTGGAAGAGTGGATGCCTGCAGTGAGCGACACCTTTAACTGACCACTGTGCGGCTTCGaa                              >  8:64768/1‑148 (MQ=255)                                                                                                                       cATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGa                            <  3:88193/147‑1 (MQ=255)                                                                                                                       cATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGAt                           <  6:254703/148‑1 (MQ=255)                                                                                                                            ataCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACAt                      >  6:293850/1‑148 (MQ=255)                                                                                                                             taCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACa                       <  3:282909/146‑1 (MQ=255)                                                                                                                                ccccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATg                 >  4:58155/1‑149 (MQ=255)                                                                                                                                ccccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATg                 >  1:4053/1‑149 (MQ=255)                                                                                                                                ccccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATg                 >  8:270715/1‑149 (MQ=255)                                                                                                                                 cccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATgg                >  4:299492/1‑149 (MQ=255)                                                                                                                                   cACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCt              >  6:151698/1‑149 (MQ=255)                                                                                                                                                cTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTGc  <  4:84411/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTGC  >  NC_002947/1175053‑1175343

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