Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,603,666	A→T	6.1%	intergenic (+997/‑505)	PP_4074 → / → PP_4076	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,603,666	0	A	T	6.1%	196.1 / 2.2	65	intergenic (+997/‑505)	PP_4074/PP_4076	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/2);  total (29/36)
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.30e-01

CAGCTCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAA  >  NC_002947/4603520‑4603810                                                                                                                                                   |                                                                                                                                                  cAGCTCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGa                                                                                                                                                 <  4:55634/149‑1 (MQ=255)    cTCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATcc                                                                                                                                              <  5:232761/149‑1 (MQ=255)    cTCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATcc                                                                                                                                              <  4:323356/149‑1 (MQ=255)     tCAGCCTTCCTGTTGGCCAGCTGCATTAAAGACAGCTTTGCCCGCTACCGCTGGTGCCCAAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAAACGCTGCGCCAGCTGCAGGCCAAGATcc                                                                                                                                              >  7:356327/1‑148 (MQ=255)             ccTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAg                                                                                                                                      <  1:28672/148‑1 (MQ=255)              cTGTTTGCCATGCGCAATTACGCCAGCTTTTCCCGCTACCGCTGGGGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAgg                                                                                                                                     <  4:221666/148‑1 (MQ=255)              cTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAgg                                                                                                                                     <  7:170346/148‑1 (MQ=255)               tGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGt                                                                                                                                    >  5:395301/1‑148 (MQ=255)               tGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGt                                                                                                                                    >  4:104164/1‑148 (MQ=255)                 tttgCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTAc                                                                                                                                  <  8:356327/145‑1 (MQ=255)                 ttGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTAc                                                                                                                                  <  4:184434/148‑1 (MQ=255)                    gCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGTCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGAt                                                                                                                              <  7:269058/149‑1 (MQ=255)                    gCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGa                                                                                                                               <  3:104164/148‑1 (MQ=255)                      cAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGt                                                                                                                                    >  5:355478/1‑141 (MQ=255)                       aGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTc                                                                                                                           <  4:359378/149‑1 (MQ=255)                                 ccGACAGCTTTGCCCGCTACCGAGGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                 <  6:358322/148‑1 (MQ=255)                                   gACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                 <  3:402306/147‑1 (MQ=255)                                   gACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATAGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                 >  4:402306/1‑147 (MQ=255)                                      aGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTAcc                                                                                                                                         <  4:69184/120‑1 (MQ=255)                                      aGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTAcc                                                                                                                                         >  3:69184/1‑120 (MQ=255)                                              ccGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAgg                                                                                                                                     >  7:252425/1‑116 (MQ=255)                                              ccGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAgg                                                                                                                                     <  8:252425/116‑1 (MQ=255)                                              ccGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTgg                                                                                                    <  6:395301/149‑1 (MQ=255)                                              ccGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTgg                                                                                                    >  6:214316/1‑149 (MQ=255)                                               cGCTACCGCTGGTGCCCCAACATCTTCGGCCCGCAATCGGGAGGTGCCGGCGAAGGCCTGCGGGTGCATCTGTACGTTTCGCTGGGCCAGCTGCAGGCCTAGTTGCGTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                 <  7:301514/135‑1 (MQ=255)                                               cGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTCCCGGTGCATCTGTACGAATAGCTGGCCCAGCTGCAGGCCAAGATACCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                 >  8:301514/1‑135 (MQ=255)                                                gCTACCGGTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGCAGGCCTGCCGGAGCATCGGTACGAATCGATGGGCCCGCTGCAGGCGAAGATCCCTAGCGGGGTACTGAGTGCGGGCCGCAAGGAGTTGGAGGTgggg                                                                                                  >  8:290324/1‑147 (MQ=255)                                                gCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGcc                                                                                                  >  3:31453/1‑149 (MQ=255)                                                gCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGcc                                                                                                  >  3:117281/1‑149 (MQ=255)                                                   aCCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCg                                                                                                 >  1:336333/1‑147 (MQ=255)                                                       cTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAgg                                                                                           <  7:314151/149‑1 (MQ=255)                                                       cTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAgg                                                                                           <  8:188556/149‑1 (MQ=255)                                                          gTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAggg                                                                                          <  5:214316/147‑1 (MQ=255)                                                          gTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCt                                                                                        <  2:222600/149‑1 (MQ=255)                                                             ccccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCttt                                                                                      >  1:240624/1‑148 (MQ=255)                                                               ccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGt                                                                                                             <  2:43895/123‑1 (MQ=255)                                                               ccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGt                                                                                                             >  1:43895/1‑123 (MQ=255)                                                               ccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATc                                                                                   >  1:238195/1‑149 (MQ=255)                                                                 aaCATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGc                                                                                 >  4:299543/1‑149 (MQ=255)                                                                    atcatcGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCt                                                                             <  4:59409/149‑1 (MQ=255)                                                                       atcGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGa                                                                           <  7:311376/148‑1 (MQ=255)                                                                          ggCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCAt                                                                       <  3:299087/149‑1 (MQ=255)                                                                              cGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATgcgc                                                                   <  3:228222/149‑1 (MQ=255)                                                                                      ggAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAg                                                           <  2:336333/149‑1 (MQ=255)                                                                                         ggTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGa                                                        <  7:290324/149‑1 (MQ=255)                                                                                          gggCCGTAGAAGGCCTTCAGGTGCTTCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGAc                                                       <  7:401823/147‑1 (MQ=255)                                                                                            gCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACaa                                                     >  4:281268/1‑149 (MQ=255)                                                                                             ccGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAAc                                                    >  3:4272/1‑149 (MQ=255)                                                                                                      gCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGttc                                           >  2:314540/1‑149 (MQ=255)                                                                                                         tGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCATGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGttcttc                                        >  2:244136/1‑149 (MQ=255)                                                                                                           ccGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTtcttct                                       <  3:281268/148‑1 (MQ=255)                                                                                                           ccGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTtctc                                      >  4:386687/1‑149 (MQ=255)                                                                                                             ggTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGTGACAACGCCGCGTTCTTCTCAg                                    <  3:299543/149‑1 (MQ=255)                                                                                                                 cATCTGTACGAATCGCTGGGCCAGCTGCATGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAACGAGTTCGAGCTGGCCGAAGATGGCTTTATCGCTCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGcc                                  <  1:314540/148‑1 (MQ=17)                                                                                                                  aTCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAAc                               >  3:74987/1‑149 (MQ=255)                                                                                                                  aTCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAAc                               >  5:367130/1‑149 (MQ=255)                                                                                                                    cTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTc                             >  3:66625/1‑149 (MQ=255)                                                                                                                                 tGGGCCAGCTGCAGGCCATGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCa                                 >  5:149220/1‑132 (MQ=255)                                                                                                                                 tGGGCCAGCTGCAGGCCATGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCa                                 <  6:149220/132‑1 (MQ=255)                                                                                                                                    gCCAGCTGCAGGCCATGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTACTCAGCCAACTCGGTGCAGAAGCCGaaa             <  1:244136/149‑1 (MQ=255)                                                                                                                                    gCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGaa              >  7:238480/1‑148 (MQ=255)                                                                                                                                            cAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTccc     <  2:6283/149‑1 (MQ=255)                                                                                                                                            cAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTccc     <  4:31453/149‑1 (MQ=255)                                                                                                                                             aGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTccc     >  4:174852/1‑148 (MQ=255)                                                                                                                                               gCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCaa  >  5:283211/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CAGCTCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGC‑CCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAA  >  NC_002947/4603520‑4603810

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