Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,506	G→C	7.4%	E308Q (GAG→CAG)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,506	0	G	C	7.4%	199.6 / 3.3	64	E308Q (GAG→CAG)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base C (2/3);  total (30/34)
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

AAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAA  >  NC_002947/4605360‑4605646                                                                                                                                                   |                                                                                                                                             aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg                                                                                                                                            <  6:275445/149‑1 (MQ=255)       gCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                      >  5:145485/1‑149 (MQ=255)            caGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGt                                                                                                                                 <  4:106025/149‑1 (MQ=255)               gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                      <  2:86596/141‑1 (MQ=255)               gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                      >  1:86596/1‑141 (MQ=255)                tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                             <  3:78674/149‑2 (MQ=255)                tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                             <  3:138934/149‑2 (MQ=255)                 gggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGc                                                                                                                            <  5:16966/149‑1 (MQ=255)                     aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                          >  3:153956/1‑147 (MQ=255)                      tCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAg                                                                                                                       >  1:174720/1‑149 (MQ=255)                              cATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                                <  3:156620/148‑1 (MQ=255)                                     cgcgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                        >  8:195788/1‑149 (MQ=255)                                         tCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGa                                                                                                                              >  5:143136/1‑123 (MQ=39)                                         tCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGa                                                                                                                              <  6:143136/123‑1 (MQ=39)                                                ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                               <  4:135321/147‑1 (MQ=255)                                                      acCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCTGAAGCGAGGCGCAGTCCGCTTACGAAAt                                                                                       <  1:20447/149‑1 (MQ=255)                                                         cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                    <  8:192912/149‑1 (MQ=255)                                                         cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                    >  7:192912/1‑149 (MQ=255)                                                         cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                     >  1:273995/1‑148 (MQ=255)                                                            cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTc                                                                                   <  7:13429/147‑1 (MQ=255)                                                            cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                  <  6:104726/148‑1 (MQ=255)                                                            cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                  <  7:195788/148‑1 (MQ=255)                                                             gcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAgt                                                                                <  6:145485/149‑1 (MQ=255)                                                                    cacaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                          >  5:296368/1‑148 (MQ=255)                                                                     acaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATg                                                                         >  2:24346/1‑148 (MQ=255)                                                                      caGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCGGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGaa                                                                       >  4:152305/1‑149 (MQ=255)                                                                        gaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCTACCTGGCAAGCCAGGCTGCCAAGGTCGAGCTGAAGCCCAGAAGCGAGGCTCAGACCGATTACGAAAAAGTCAGTGTcatagaag                                                                      >  8:288901/1‑141 (MQ=255)                                                                         aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGt                                                                     >  2:290728/1‑148 (MQ=255)                                                                         aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGt                                                                     >  5:83271/1‑148 (MQ=255)                                                                                 aGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAc                                                            >  8:236290/1‑149 (MQ=255)                                                                                   ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                          <  2:273995/149‑1 (MQ=255)                                                                                            cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGTACAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTGAGTGTCGATGAAGTCATCAGCACCCTaagaacga                                                 <  2:174720/149‑3 (MQ=255)                                                                                            cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAAcc                                                  <  8:45415/148‑1 (MQ=255)                                                                                             aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCAc                                                >  1:155948/1‑149 (MQ=255)                                                                                              ccccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                               <  6:296368/149‑1 (MQ=255)                                                                                                ccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGa                                             <  5:237392/149‑1 (MQ=255)                                                                                                 cGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGAc                                            <  4:153956/149‑1 (MQ=255)                                                                                                    cTGTCAACCTGTTCAAGCACAGTGCCTAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACg                                           <  6:83271/147‑1 (MQ=255)                                                                                                     tctCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGaa                                        <  1:290728/147‑1 (MQ=255)                                                                                                       tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                  >  8:249407/1‑105 (MQ=255)                                                                                                       tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                  <  7:249407/105‑1 (MQ=255)                                                                                                         aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAg                                                                      <  1:46936/115‑1 (MQ=255)                                                                                                         aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAg                                                                      >  2:46936/1‑115 (MQ=255)                                                                                                         aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                    >  4:277600/1‑149 (MQ=255)                                                                                                          aCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGaa                                        >  1:165010/1‑144 (MQ=255)                                                                                                          aCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGaa                                        <  2:165010/144‑1 (MQ=255)                                                                                                          aCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTa                                     >  2:261856/1‑147 (MQ=255)                                                                                                          aCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTa                                     <  1:261856/147‑1 (MQ=255)                                                                                                          aCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAcc                                   >  1:11259/1‑149 (MQ=255)                                                                                                           ccTGTTCAAGCTCAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGa                                                                        >  5:60245/1‑111 (MQ=255)                                                                                                           ccTGTTCAAGCTCAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGa                                                                        <  6:60245/111‑1 (MQ=255)                                                                                                           ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                                  <  7:236290/149‑1 (MQ=255)                                                                                                             tGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACCGAAGTACCGGt                                >  2:138144/1‑149 (MQ=255)                                                                                                               ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGa                              <  5:125440/149‑1 (MQ=255)                                                                                                                  aaGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                               <  5:26840/81‑1 (MQ=255)                                                                                                                  aaGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                               >  6:26840/1‑81 (MQ=255)                                                                                                                  aaGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCa                           >  4:238001/1‑149 (MQ=255)                                                                                                                          tGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGa                                         <  7:238066/127‑1 (MQ=255)                                                                                                                          tGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGACCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGa                                         >  8:238066/1‑127 (MQ=255)                                                                                                                                    aTCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGccc         >  4:194210/1‑149 (MQ=255)                                                                                                                                    aTCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGccc         >  5:297030/1‑149 (MQ=255)                                                                                                                                           tGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGc                      <  7:127988/129‑1 (MQ=255)                                                                                                                                           tGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGc                      >  8:127988/1‑129 (MQ=255)                                                                                                                                             gCCAGCGAGGCTGCCCGGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGaa  <  3:152305/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAA  >  NC_002947/4605360‑4605646

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