Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,525,770	C→A	9.4%	intergenic (+60/‑63)	PP_4858 → / → PP_4859	hypothetical protein/TetR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,525,770	0	C	A	9.4%	192.0 / 9.3	64	intergenic (+60/‑63)	PP_4858/PP_4859	hypothetical protein/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base A (2/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

AGGCGGCTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGTTT  >  NC_002947/5525623‑5525918                                                                                                                                                     |                                                                                                                                                     aGGCGGCTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCa                                                                                                                                                     <  8:72732/149‑1 (MQ=255)        tCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCg                                                                                                                                              <  7:195744/149‑1 (MQ=255)           gATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                           <  4:49638/149‑1 (MQ=255)               gTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCAATGGACCGGGcctg                                                                                                                                        <  2:196571/148‑4 (MQ=255)                tCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGt                                                                                                                                      >  6:19618/1‑149 (MQ=255)                    aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGtt                                                                                                                                     <  1:37439/147‑1 (MQ=255)                    aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                  >  5:45603/1‑149 (MQ=255)                       gAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTa                                                                                                                               >  1:312638/1‑149 (MQ=255)                       gAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTa                                                                                                                               >  7:60227/1‑149 (MQ=255)                         aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATc                                                                                                                             >  5:58763/1‑149 (MQ=255)                         aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                           <  4:295934/135‑1 (MQ=255)                         aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                           >  3:295934/1‑135 (MQ=255)                             gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTaga                                                                                                                         <  6:58763/149‑1 (MQ=255)                             gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAAGCCCGGTCAATGGACCGGGcct                                                                                                                                         <  6:86110/133‑3 (MQ=255)                             gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAAGCCCGGTCAATGGACCGGGcct                                                                                                                                         >  5:86110/1‑131 (MQ=255)                                   cGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGccc                                                                                                                   <  6:45603/149‑1 (MQ=255)                                       tCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGc                                                                                                                <  2:312638/148‑1 (MQ=255)                                           tGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGAc                                                                                                           >  4:121367/1‑149 (MQ=255)                                            gAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCC‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGAcc                                                                                                          >  2:171507/1‑149 (MQ=255)                                             aaGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCg                                                                                                         >  3:207393/1‑149 (MQ=255)                                              aGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGt                                                                                                        <  4:207393/149‑1 (MQ=255)                                                        gATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                  <  5:215819/113‑1 (MQ=255)                                                        gATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                  >  6:215819/1‑113 (MQ=255)                                                            gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGca                                                                                          >  2:1267/1‑149 (MQ=255)                                                            gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                           <  8:60227/148‑1 (MQ=255)                                                            gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                           <  5:19618/148‑1 (MQ=255)                                                             ccTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                           >  6:48601/1‑115 (MQ=255)                                                             ccTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                           <  5:48601/115‑1 (MQ=255)                                                                      aCTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCC‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGcc                                                                                <  1:171507/149‑1 (MQ=255)                                                                        tCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgc                                                                              >  7:76835/1‑149 (MQ=255)                                                                        tCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCAAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgc                                                                              >  5:287783/1‑149 (MQ=255)                                                                            ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGcatcat                                                                                      >  1:18897/1‑137 (MQ=255)                                                                            ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGcatcat                                                                                      <  2:18897/137‑1 (MQ=255)                                                                            ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCgcgcgc                                                                          <  4:235663/149‑1 (MQ=255)                                                                                      tGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAg                                                                >  3:43377/1‑149 (MQ=255)                                                                                                 gtccaaaaGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGccc                                                                                             <  1:289722/105‑1 (MQ=255)                                                                                                 gCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGccc                                                                                             >  2:289722/1‑109 (MQ=255)                                                                                                  ccGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCg                                                                        <  4:227354/129‑1 (MQ=255)                                                                                                  ccGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCg                                                                        >  3:227354/1‑129 (MQ=255)                                                                                                          aacaacTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCt                                            >  6:327873/1‑149 (MQ=255)                                                                                                            caacTTCTGTCACAACGCTTTCT‑CCCCTTCATGGCCGGTTCATTGACCGGGGTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGc                                               <  8:124248/144‑1 (MQ=255)                                                                                                            caacTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGc                                          >  6:271768/1‑149 (MQ=255)                                                                                                            caacTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAAAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGATGCCGCGCGCCGACCAGAAGCTACAGACCCGCCAGcccag                                            >  7:124248/1‑145 (MQ=255)                                                                                                                 ccgagctcaACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTGTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGa                                      <  3:253670/141‑1 (MQ=255)                                                                                                                   tGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgc                                   >  7:8627/1‑149 (MQ=255)                                                                                                                    gTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTa                                                             >  2:13400/1‑122 (MQ=255)                                                                                                                     tCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGATTTCCCAGCCTGACCGTTACAGAACCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCAGCACGCTgccg                                 >  8:117942/1‑149 (MQ=255)                                                                                                                     tCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccg                                 >  5:184144/1‑149 (MQ=255)                                                                                                                     tCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgc                                   <  4:43377/147‑1 (MQ=255)                                                                                                                       acaACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                <  8:8627/148‑1 (MQ=255)                                                                                                                       acaACGCTTGCT‑CCCCTCCAGGCCAGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgcc                               >  4:319043/1‑149 (MQ=255)                                                                                                                         aaCGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCgcctg                             >  2:7835/1‑149 (MQ=255)                                                                                                                         aaCGCTTGCT‑CCCCTCCAGGCCCGGACCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGAATCATCCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGcccggagcga                                      >  4:253670/1‑133 (MQ=255)                                                                                                                            gCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCgcctgcct                          >  8:300118/1‑149 (MQ=255)                                                                                                                                     cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                >  7:6180/1‑135 (MQ=255)                                                                                                                                     cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                <  8:6180/135‑1 (MQ=255)                                                                                                                                     cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATg                   <  5:271768/148‑1 (MQ=255)                                                                                                                                      ccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                 <  3:121367/149‑1 (MQ=255)                                                                                                                                      ccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                 <  8:76835/149‑1 (MQ=255)                                                                                                                                      ccTCCAGGCCAGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                 <  3:319043/149‑1 (MQ=255)                                                                                                                                         ccAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGagag              >  8:182168/1‑149 (MQ=255)                                                                                                                                          cAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGgaga               <  7:117942/147‑1 (MQ=255)                                                                                                                                            ggCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTgg           <  5:327873/149‑1 (MQ=255)                                                                                                                                                   tCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACACGCCAGGCCCTGCTCGATGCCgcctgccag                         >  3:253620/1‑126 (MQ=255)                                                                                                                                                     cATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGttt  <  1:7835/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGCGGCTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCT‑CCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGTTT  >  NC_002947/5525623‑5525918

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