Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,004,819	A→G	7.8%	H467R (CAC→CGC)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,004,819	0	A	G	7.8%	179.5 / 4.3	64	H467R (CAC→CGC)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

AGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGT  >  NC_002947/3004672‑3004957                                                                                                                                                    |                                                                                                                                            aGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCAc                                                                                                                                            <  6:49980/149‑1 (MQ=255)        cacGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGTGGCGCTTGCtgc                                                                                                                                     <  6:241669/149‑2 (MQ=255)         acgacgAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCtgtg                                                                                                                                    <  2:65110/149‑1 (MQ=255)         acgacgAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCtgtg                                                                                                                                    >  4:161115/1‑149 (MQ=255)           gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtggg                                                                                                                                  >  6:205985/1‑149 (MQ=255)           gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtggg                                                                                                                                  >  7:5746/1‑149 (MQ=255)            acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 <  5:133842/149‑1 (MQ=255)            acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 >  5:115414/1‑149 (MQ=255)            acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 <  4:109358/149‑1 (MQ=255)             cgaAGCTGCCCGTGCGTGTAGGTAAGGCACAGGGGAAGGCGTTTGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGTGCTGGCGGGTGACAATGTTCTGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                <  2:35439/149‑1 (MQ=255)             cgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                >  4:203798/1‑149 (MQ=255)              gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                               >  1:69553/1‑149 (MQ=255)              gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                               <  3:161115/149‑1 (MQ=255)                     ccccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑g                                                                                                                        <  3:203798/148‑1 (MQ=255)                      cccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑Ggt                                                                                                                      <  5:205985/149‑1 (MQ=255)                      cccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑Ggt                                                                                                                      <  2:143049/149‑1 (MQ=255)                        cTGCTTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGTAAGCGCCACTTGCTGTGGGTGGGCTCT‑Ggtgt                                                                                                                    <  5:232314/149‑1 (MQ=255)                               aaGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATATTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATc                                                                                                             >  8:43716/1‑149 (MQ=255)                                     aGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGc                                                                                                        <  2:69553/148‑1 (MQ=255)                                        cacaGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGtt                                                                                                     >  7:60984/1‑148 (MQ=255)                                             gggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTAt                                                                                               >  5:321245/1‑149 (MQ=255)                                              ggcaGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATc                                                                                              <  4:250337/146‑1 (MQ=255)                                              ggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATc                                                                                              <  6:115414/149‑1 (MQ=255)                                                 aGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGc                                                                                            <  8:60984/148‑1 (MQ=255)                                                    cGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGCAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGtt                                                                                        >  1:234460/1‑149 (MQ=255)                                                         gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGc                                                                                      <  7:56419/146‑1 (MQ=255)                                                         gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGc                                                                                      >  8:56419/1‑146 (MQ=255)                                                         gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAg                                                                                   >  2:94366/1‑149 (MQ=255)                                                            tATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTg                                                                                >  1:78492/1‑149 (MQ=255)                                                              tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  4:234502/1‑148 (MQ=255)                                                              tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCt                                                                                 <  6:321245/146‑1 (MQ=255)                                                              tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  4:137979/1‑148 (MQ=255)                                                                  aTCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCg                                                                          >  4:130061/1‑149 (MQ=255)                                                                   tCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGt                                                                         >  5:336384/1‑149 (MQ=255)                                                                         aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  1:322921/1‑137 (MQ=255)                                                                         aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               <  2:322921/137‑1 (MQ=255)                                                                         aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGcc                                                                   >  1:328750/1‑149 (MQ=255)                                                                         aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGcc                                                                   <  2:78492/149‑1 (MQ=255)                                                                              ttCAGGCGGGTTATCTATCAGGAAGCAGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACa                                                               >  1:178395/1‑148 (MQ=255)                                                                                 aGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCaa                                                           >  7:96686/1‑149 (MQ=255)                                                                                  ggCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                          <  3:130061/149‑1 (MQ=255)                                                                                  ggCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                          >  7:97287/1‑149 (MQ=255)                                                                                    cGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATgg                                                        <  6:336384/149‑1 (MQ=255)                                                                                     gggTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                       <  1:94366/149‑1 (MQ=255)                                                                                      ggTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgcg                                                      <  8:5746/149‑1 (MQ=255)                                                                                              atcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCAAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCt                                              >  8:295842/1‑149 (MQ=255)                                                                                               tcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTGGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAg                                                 <  2:234460/145‑1 (MQ=255)                                                                                                cAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            >  6:156056/1‑149 (MQ=255)                                                                                                 aGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            <  8:97287/148‑1 (MQ=255)                                                                                                 aGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            >  2:269968/1‑148 (MQ=255)                                                                                                   gAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAg                                         >  8:222561/1‑149 (MQ=255)                                                                                                                 gggTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCaaaa                           <  2:262479/149‑1 (MQ=255)                                                                                                                 gggTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCaaaa                           <  1:269968/149‑1 (MQ=255)                                                                                                                      aCATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                      <  3:234502/149‑1 (MQ=255)                                                                                                                        aaTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                      >  5:64650/1‑147 (MQ=255)                                                                                                                           gTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTa                 >  4:71456/1‑149 (MQ=255)                                                                                                                                  ttGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGt          >  8:88470/1‑149 (MQ=255)                                                                                                                                      gcgcAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc       >  8:90622/1‑148 (MQ=255)                                                                                                                                      ccacagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc       >  1:227172/4‑149 (MQ=255)                                                                                                                                        gcagcaAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGaa    <  2:328750/149‑1 (MQ=255)                                                                                                                                        acagcaAGTGGCGCTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACg      >  6:222252/2‑147 (MQ=255)                                                                                                                                          agcaAGTGGCGCTTGCTGCGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  2:53207/149‑1 (MQ=255)                                                                                                                                          agcaAGTGGCGCTTGCTGCGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  2:48347/149‑1 (MQ=255)                                                                                                                                          agcaAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  6:64650/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            AGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCT‑GGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGT  >  NC_002947/3004672‑3004957

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