Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,787,501	Δ1 bp	40.0%	coding (1134/1392 nt)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,787,501	0	G	.	40.0%	143.7 / 106.8	70	coding (1134/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base . (15/13);  total (37/33)
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.58e-01

GGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTC  >  NC_002947/3787355‑3787647                                                                                                                                                   |                                                                                                                                                    ggAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATgc                                                                                                                                                    <  5:101426/148‑1 (MQ=255)         gcgCCAGCACGCGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTgcgc                                                                                                                                         >  6:52562/3‑149 (MQ=255)           gCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTgcgc                                                                                                                                         >  3:303046/1‑149 (MQ=255)             cAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAgg                                                                                                                                     >  5:110774/1‑149 (MQ=255)              aGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAggg                                                                                                                                    >  6:332050/1‑149 (MQ=255)                 gcgTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGc                                                                                                                                   >  7:260936/1‑149 (MQ=255)                   gTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCgag                                                                                                                                >  2:33484/1‑148 (MQ=255)                      cTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGc                                                                                                                            <  5:332050/149‑1 (MQ=255)                       tGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAg                                                                                                                             >  6:244202/1‑149 (MQ=255)                         gcgcGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGa                                                                                                                         <  7:82618/149‑1 (MQ=255)                          cgcgGTACATGGTACCGACATGCCGACCCGGGCCGCGGGTGTTGGCCAGCATCCTCATGCCGGGCATGAACACGCCGTCTTCCACACGCTCTTTCCACACTCCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGaa                                                                                                                        <  4:49089/149‑1 (MQ=255)                              gTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTc                                                                                                                    >  1:259575/1‑149 (MQ=255)                              gTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTc                                                                                                                    <  6:110774/149‑1 (MQ=255)                               tAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCg                                                                                                                   >  6:57899/1‑149 (MQ=255)                                   tGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGc                                                                                                                 <  5:263559/149‑1 (MQ=255)                                      tACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCACAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGccc                                                                                                               >  6:161612/1‑148 (MQ=255)                                        ccGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCccacca                                                                                                          >  2:107968/1‑149 (MQ=255)                                         cGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCccacc                                                                                                           <  7:220303/149‑1 (MQ=255)                                         cGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCcaccac                                                                                                         >  2:95823/1‑149 (MQ=255)                                          gACAAGGAGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCCCGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCACGGCGAGAAGCCGAAGGTCGGCCccacca                                                                                                          >  8:284596/1‑149 (MQ=255)                                           aCATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCcaccac                                                                                                         >  5:105587/1‑149 (MQ=255)                                                 cGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGt                                                                                                 >  7:110900/1‑149 (MQ=255)                                                  gACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGatgatg                                                                                                  >  4:190459/1‑149 (MQ=255)                                                  caCCCGCGGCACGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGACGGCCCCACCACGatgatg                                                                                                  >  5:325020/2‑149 (MQ=255)                                                    cccGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGt                                                                                                 <  6:321459/148‑1 (MQ=255)                                                          ccGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCa                                                                                           <  1:390096/148‑1 (MQ=255)                                                          ccGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCAC‑TCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAgg                                                                                         <  1:107968/149‑1 (MQ=255)                                                             cgGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGggt                                                                                       >  2:246960/1‑149 (MQ=255)                                                              gggTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGggtggt                                                                                    <  5:57899/149‑1 (MQ=255)                                                              gggTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGCCCCCACCACGATGATGTTGCTCAGggtggt                                                                                    <  1:33484/149‑1 (MQ=255)                                                                gtggtgGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGggtggt                                                                                    <  8:185327/149‑1 (MQ=255)                                                                  gtTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGggtggt                                                                                    >  7:142462/1‑147 (MQ=255)                                                                  gtTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGa                                                                                  <  2:333448/149‑1 (MQ=255)                                                                   ttGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGa                                                                                  <  7:167672/148‑1 (MQ=255)                                                                   ttGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAg                                                                                 <  8:142462/149‑1 (MQ=255)                                                                       ccAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAggctg                                                                             >  6:115260/1‑149 (MQ=255)                                                                          gCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTc                                                                        <  1:95823/149‑1 (MQ=255)                                                                           cATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCtt                                                                         >  4:16524/1‑149 (MQ=255)                                                                           cATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCtt                                                                         >  1:208261/1‑149 (MQ=255)                                                                             tGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCg                                                                       <  3:16524/149‑1 (MQ=255)                                                                               cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                      <  8:260936/148‑1 (MQ=255)                                                                               cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGc                                                                                                                 <  1:185595/103‑1 (MQ=255)                                                                               cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGc                                                                                                                 >  2:185595/1‑103 (MQ=255)                                                                                tCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑a                                                                   >  3:148846/1‑149 (MQ=255)                                                                                tCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑a                                                                   >  2:322764/1‑147 (MQ=255)                                                                                  aTGCCGGGCATGTACACGCCGTCTTCCTCAGGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATa                                                                 <  7:183827/147‑1 (MQ=255)                                                                                               acacGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCGATGGGCACCATTGCGTg                                                   >  1:136645/1‑149 (MQ=255)                                                                                                cacGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  <  5:52562/148‑1 (MQ=255)                                                                                                 acGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  <  4:148846/149‑1 (MQ=255)                                                                                                   gCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  <  8:215920/147‑1 (MQ=255)                                                                                                   gCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  >  7:215920/1‑147 (MQ=255)                                                                                                      gTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGc                                           <  2:259575/149‑1 (MQ=255)                                                                                                         ttCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCggg                                        >  6:185075/1‑149 (MQ=255)                                                                                                           ccTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGc                                           <  2:208261/146‑1 (MQ=255)                                                                                                           ccTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCggg                                        >  6:369893/1‑149 (MQ=255)                                                                                                             tCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTc                                      >  4:194315/1‑149 (MQ=255)                                                                                                               aCGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  <  8:59868/135‑1 (MQ=255)                                                                                                               aCGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTgg                                                  >  7:59868/1‑135 (MQ=255)                                                                                                               aCGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAg                                  <  8:110900/149‑1 (MQ=255)                                                                                                                 gCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAg                                  <  6:325020/149‑1 (MQ=255)                                                                                                                 gCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCACAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAg                                  <  5:161612/149‑1 (MQ=255)                                                                                                                  cTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑a                                                                   <  4:24767/115‑1 (MQ=255)                                                                                                                  cTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑a                                                                   >  3:24767/1‑115 (MQ=255)                                                                                                                         cacaCTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGAc                          >  5:220339/1‑149 (MQ=255)                                                                                                                         cacaCTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGaa                                   <  6:67580/138‑1 (MQ=255)                                                                                                                         cacaCTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGaa                                   >  5:67580/1‑138 (MQ=255)                                                                                                                          acacTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACttt                       >  4:99513/1‑149 (MQ=255)                                                                                                                            acTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACttt                       <  3:190459/149‑1 (MQ=255)                                                                                                                             cTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGTCACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtttt                      <  5:115260/149‑1 (MQ=255)                                                                                                                                                gATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGt   >  5:300866/1‑149 (MQ=255)                                                                                                                                                 aTGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTc  >  8:62074/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCC‑ATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTC  >  NC_002947/3787355‑3787647

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