Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	564,382	Δ1 bp	11.4%	intergenic (+34/‑139)	rplQ → / → katA	50S ribosomal protein L17/catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	564,382	0	T	.	11.4%	206.6 / 22.0	70	intergenic (+34/‑139)	rplQ/katA	50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base . (5/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGA  >  NC_002947/564237‑564524                                                                                                                                                   |                                                                                                                                               cggcTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGcc                                                                                                                                             <  8:157473/149‑1 (MQ=255)       ccTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGttt                                                                                                                                       <  6:13253/149‑1 (MQ=255)         tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGttttt                                                                                                                                     >  4:203554/1‑149 (MQ=255)            gcATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGt                                                                                                                                  <  2:983/149‑1 (MQ=255)             cATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                 <  4:24905/149‑2 (MQ=255)                    aaGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACAt                                                                                                                          <  3:17112/149‑1 (MQ=255)                          ggTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTg                                                                                                                    >  4:256349/1‑149 (MQ=255)                           gTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGt                                                                                                                   <  1:233429/149‑1 (MQ=255)                           gTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGt                                                                                                                   <  1:228812/149‑1 (MQ=255)                            tttCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGtt                                                                                                                  >  7:53087/1‑149 (MQ=255)                                  cgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTc                                                                                                            >  3:167228/1‑149 (MQ=255)                                     tGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                             >  6:89238/1‑145 (MQ=255)                                     tGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                             <  5:89238/145‑1 (MQ=255)                                       gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCTTCCGGTCGGCGGTTCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                         <  7:139515/147‑1 (MQ=255)                                       gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTTCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                        <  5:24303/148‑1 (MQ=255)                                       gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                        <  6:232965/148‑1 (MQ=255)                                         gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGAt                                                                                                     <  5:152186/149‑1 (MQ=255)                                               gcgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGACGGCGGTGCTGGAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAg                                                                                                >  1:221221/1‑148 (MQ=255)                                                cgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGt                                                                                               >  7:229329/1‑148 (MQ=255)                                                        gCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                       >  6:303412/1‑148 (MQ=255)                                                                 gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                             <  5:303412/149‑1 (MQ=255)                                                                 gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                             <  8:53087/149‑1 (MQ=255)                                                                 gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                 >  3:59545/1‑96 (MQ=37)                                                                 gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                 <  4:59545/97‑2 (MQ=37)                                                                   gCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATGTATTTTTTATTTATTATTATGTTATTCATTTAATGAATTTGtacaa                                                                           >  7:257747/1‑145 (MQ=255)                                                                   gCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                           <  4:167228/149‑1 (MQ=255)                                                                    gtGGTTGCTCGGCCGGTCGGCGGGGCTGTAGAATCTGCTGAGTAAGACGTTCAGTCTGCTACCGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTAATGGTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                           <  5:117872/147‑1 (MQ=255)                                                                          gATCGTCCGGTCGTCGGTGCTGTAGAAGCTGCTTAGTAAGACGTTCTGTCTGCTACAGAAAACCGTGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                     <  7:166560/148‑1 (MQ=255)                                                                           aTCGTCCGGTAGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCtt                                                                   <  5:130507/149‑1 (MQ=255)                                                                                   gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                     <  7:191370/139‑1 (MQ=255)                                                                                   gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                     >  8:191370/1‑139 (MQ=255)                                                                                      ggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATa                                                        >  3:79460/1‑149 (MQ=255)                                                                                           tGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctt                                                   <  8:221134/149‑2 (MQ=255)                                                                                              tGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc                                                 >  4:64481/1‑148 (MQ=255)                                                                                               gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctct                                                >  3:91372/1‑148 (MQ=255)                                                                                               gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                  <  6:309730/146‑1 (MQ=255)                                                                                               gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                  >  5:309730/1‑146 (MQ=255)                                                                                                tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                               >  7:172569/1‑148 (MQ=255)                                                                                                 aGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                               >  2:108047/1‑147 (MQ=255)                                                                                                  gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAg                                            >  2:235365/1‑149 (MQ=255)                                                                                                      ctgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGa                                        <  7:275937/149‑1 (MQ=255)                                                                                                       tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                       >  5:325224/1‑149 (MQ=255)                                                                                                       tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                       >  7:256262/1‑149 (MQ=255)                                                                                                        gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                       >  3:63291/1‑148 (MQ=255)                                                                                                        gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACACTCTCTACAGCCGAt                                       >  4:251117/1‑148 (MQ=255)                                                                                                              tAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                >  4:98079/1‑149 (MQ=255)                                                                                                              tAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                >  4:93552/1‑149 (MQ=255)                                                                                                                    gTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCtt                           >  8:205413/1‑148 (MQ=255)                                                                                                                     ttCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGa                         <  1:235365/149‑1 (MQ=255)                                                                                                                       ctttctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                        <  3:226996/145‑1 (MQ=255)                                                                                                                          tctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAcc                    >  4:59797/1‑149 (MQ=255)                                                                                                                              cTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGaaa                <  2:221221/149‑1 (MQ=255)                                                                                                                                aCAAAAAACCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg              >  1:39429/1‑149 (MQ=255)                                                                                                                                 cAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAAt               <  7:205413/147‑1 (MQ=255)                                                                                                                                  aaaaaaCCGGGCACTAG‑GCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                       <  3:164533/72‑1 (MQ=25)                                                                                                                                  aaaaaaCCGGGCACTAG‑GCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                       >  4:164533/3‑74 (MQ=25)                                                                                                                                  aGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg              >  7:251663/1‑147 (MQ=255)                                                                                                                                  aGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg              >  5:138882/1‑147 (MQ=255)                                                                                                                                   gAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgag           <  4:63291/149‑1 (MQ=255)                                                                                                                                   aaaaaCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg              >  6:199768/2‑146 (MQ=255)                                                                                                                                    aaaaCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTAATGTAATGAATTTGTCACTGTAATCTTGCTCGTCAATACTCTCTCTACAGAAGATTAGCCGGCCCTTGAAGACAGAAATGgaga          >  1:128446/1‑149 (MQ=255)                                                                                                                                    aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                 >  8:86617/1‑94 (MQ=255)                                                                                                                                    aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                 <  7:86617/94‑1 (MQ=255)                                                                                                                                    aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCTCTTGAAGACCGAAATGga            >  6:16193/1‑147 (MQ=255)                                                                                                                                    aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga            >  7:251278/1‑147 (MQ=255)                                                                                                                                    aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga            >  4:164586/1‑147 (MQ=255)                                                                                                                                         cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga            <  3:164988/142‑1 (MQ=255)                                                                                                                                         cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga            >  4:164988/1‑142 (MQ=255)                                                                                                                                          gggC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt    <  4:91372/149‑1 (MQ=255)                                                                                                                                            gC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGa  <  6:325224/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               CGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGA  >  NC_002947/564237‑564524

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