Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,203,423	T→A	8.7%	S58R (AGT→AGA)	PP_2810 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,203,423	0	T	A	8.7%	179.2 / 5.7	65	S58R (AGT→AGA)	PP_2810	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (2/4);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGG  >  NC_002947/3203276‑3203556                                                                                                                                                    |                                                                                                                                      aaGTGTCGTATCCGCCCGCCGTATTCTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGtt                                                                                                                                      <  1:216774/149‑1 (MQ=255)         tATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgcggc                                                                                                                              >  6:363515/1‑149 (MQ=255)         tATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgcggc                                                                                                                              <  3:174334/149‑1 (MQ=255)           tCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGa                                                                                                                            >  1:329743/1‑149 (MQ=255)             cgcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCTAGAACTGgcggcgctc                                                                                                                          <  3:334078/149‑4 (MQ=255)             cgcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCTAGAACTGgcggcgct                                                                                                                           <  7:126271/148‑3 (MQ=255)             cgcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGata                                                                                                                          <  2:329743/149‑1 (MQ=255)                ccgccgTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTc                                                                                                                       >  6:414677/1‑149 (MQ=255)                   ccgTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAg                                                                                                                    <  8:75067/149‑1 (MQ=255)                    cgTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGc                                                                                                                   <  6:275608/149‑1 (MQ=255)                         tttCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTggg                                                                                                              >  6:250010/1‑149 (MQ=255)                             aGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGca                                                                                                            >  3:310036/1‑147 (MQ=255)                               gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                         >  8:395375/1‑132 (MQ=255)                               gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                         <  7:395375/132‑1 (MQ=255)                                    tgctgGTTACCCTTGCGTCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAggg                                                                                                   <  3:260549/149‑1 (MQ=255)                                    tgctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTgg                                                                                                               >  8:363218/1‑137 (MQ=255)                                    tgctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTgg                                                                                                               <  7:363218/137‑1 (MQ=255)                                    tgctgGTTACCCTTGCAGCCTCGACCGCGATCTGTGCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTTGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTGTAGCGGCGATATTCCAGCTTGGGCACAGCCAggg                                                                                                   <  7:404688/149‑1 (MQ=255)                                      ctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAggg                                                                                                   >  4:87013/1‑147 (MQ=255)                                         gTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGgcag                                                                                                >  4:321447/1‑147 (MQ=255)                                           tACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACTCAGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACTGCCAGGGCAGCAgg                                                                                            >  7:280363/1‑149 (MQ=255)                                            aCCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAgg                                                                                            <  4:187881/148‑1 (MQ=255)                                             cccTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAgggg                                                                                          <  3:276106/149‑1 (MQ=255)                                              ccTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCCGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCAGCCAGTTCTAGCGGCGATATTCCCGCCTGGGCCCCGCCAGGGCAGCAggggg                                                                                         >  6:439054/1‑149 (MQ=255)                                                  gCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGAGGAGGCGCTGGGTTCGAAGTTGACGCCGCGGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTACAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAg                                                                                      >  6:99458/1‑148 (MQ=255)                                                   cGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAgg                                                                                     <  5:414677/148‑1 (MQ=255)                                                   cGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGa                                                                                    >  7:185231/1‑149 (MQ=255)                                                       cTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATg                                                                                  >  7:49633/1‑147 (MQ=255)                                                       cTCGACCGCGATCAGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGa                                                                                <  5:101299/149‑1 (MQ=255)                                                        tCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGACGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAg                                                                               >  6:6413/1‑149 (MQ=255)                                                           aCCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTa                                                                            >  8:264079/1‑149 (MQ=255)                                                                gATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGc                                                                        >  8:214097/1‑148 (MQ=255)                                                                     gTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCTGCAGGGGGCAGGATGGAGTTACGGCCAGGTa                                                                  >  6:194714/1‑149 (MQ=255)                                                                      ttCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAg                                                                     <  5:250010/145‑1 (MQ=255)                                                                      ttCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAg                                                                     <  4:18999/145‑1 (MQ=255)                                                                               gcctccgcCCCCTGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAg                                                         >  1:329596/6‑148 (MQ=255)                                                                                cctccgcCCCCTGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCa                                                       >  5:195600/5‑149 (MQ=255)                                                                                  tGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACt                                                     >  8:368989/1‑149 (MQ=255)                                                                                  tGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACt                                                     >  8:105758/1‑149 (MQ=255)                                                                                  tGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGATGGGGGGCGAGAGCGACGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGCGGGCAGGATGGAGTTACGGCCAGGTAGGAGGGGAGCACt                                                     >  8:394302/1‑149 (MQ=255)                                                                                    cgcCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTgg                                                   <  5:363515/149‑1 (MQ=255)                                                                                    cgcCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTg                                                    <  7:105758/148‑1 (MQ=255)                                                                                      ccccAGGGGGCGGAGGCGCTGGGTTGGAAGTTGCCGCCGCGGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCGGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGaa                                                 <  5:439054/149‑1 (MQ=255)                                                                                         cAGGGGGCGGAGGCGCTGGGTTCTAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGtt                                              <  8:295079/149‑1 (MQ=255)                                                                                             gggCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAg                                          <  3:321447/149‑1 (MQ=255)                                                                                             cccctggggCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAg                                          <  3:282431/142‑1 (MQ=255)                                                                                                 ggAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGt                                       >  1:405781/1‑148 (MQ=255)                                                                                                       gcTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGc                                <  7:214097/149‑1 (MQ=255)                                                                                                            gTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCctc                           >  5:243376/1‑149 (MQ=255)                                                                                                             ttCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCctc                           >  1:371755/1‑148 (MQ=255)                                                                                                               cGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGAGGGCGTTACGGCCAGGTACGAGGGGAGCACGGGACGTTCATGGGGGATAGGCCCcgctc                         >  5:147374/1‑147 (MQ=255)                                                                                                                gAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCTGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGATCACTGGAAGTTCAAGGGTGATAAGCCCCTCTac                        <  5:194714/148‑1 (MQ=255)                                                                                                                      gACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATc                  <  3:87013/148‑1 (MQ=255)                                                                                                                        cgccgcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCTATATTCCAGCCTGGGCACAGCGAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCACGGGTGATAAGCCCCTCTACAACATCGAc               <  5:99458/149‑1 (MQ=255)                                                                                                                        cgccgcGTGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGAc               <  8:280363/149‑1 (MQ=255)                                                                                                                         gccgcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACt              >  6:68005/1‑149 (MQ=255)                                                                                                                            gcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGa           <  2:371755/149‑1 (MQ=255)                                                                                                                             cTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGaa          <  8:185231/149‑1 (MQ=255)                                                                                                                             cTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGaa          <  6:243376/149‑1 (MQ=255)                                                                                                                             cTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGaa          >  6:199166/1‑149 (MQ=255)                                                                                                                             cTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGaa          <  8:49633/149‑1 (MQ=255)                                                                                                                                      gagaGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg  >  1:233156/1‑148 (MQ=255)                                                                                                                                          tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg  <  2:57854/143‑1 (MQ=255)                                                                                                                                          tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg  <  2:51490/143‑1 (MQ=255)                                                                                                                                          tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg  >  1:57854/2‑144 (MQ=255)                                                                                                                                          tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg  >  1:51490/2‑144 (MQ=255)                                                                                                                                                    |                                                                                                                                      AAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGG  >  NC_002947/3203276‑3203556

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