Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,389,000	C→G	8.2%	G581R (GGC→CGC)	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,389,000	0	C	G	8.2%	193.0 / 6.6	61	G581R (GGC→CGC)	PP_3865	tail protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base G (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

CACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGG  >  NC_002947/4388857‑4389148                                                                                                                                                |                                                                                                                                                     cACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACGTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGccc                                                                                                                                                  <  5:88009/148‑1 (MQ=255)       cTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCACCGGGCGGCCTGCTGCgg                                                                                                                                           <  6:352303/149‑1 (MQ=255)         cGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCggg                                                                                                                                          >  2:53287/1‑148 (MQ=255)          gACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCggg                                                                                                                                          >  3:113145/1‑147 (MQ=255)               gCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                                   <  6:238718/149‑1 (MQ=255)                  ggccaggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTc                                                                                                                                <  2:47235/149‑1 (MQ=255)                     caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCTGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTcgcc                                                                                                                             <  5:307902/149‑1 (MQ=255)                     caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                                   >  6:114822/1‑143 (MQ=255)                     caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                                   <  5:114822/143‑1 (MQ=255)                          ggcgAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCcgcacg                                                                                                                        >  5:90901/1‑149 (MQ=255)                                ttATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGgcc                                                                                                                                                      <  2:349292/113‑1 (MQ=255)                                 tATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAgc                                                                                                                 >  2:289243/1‑149 (MQ=255)                                   tCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCgg                                                                                                                                           >  8:315737/1‑121 (MQ=255)                                   tCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCgg                                                                                                                                           <  7:315737/121‑1 (MQ=255)                                    cGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGcc                                                                                                              >  7:194372/1‑149 (MQ=255)                                    cGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGcc                                                                                                              >  3:305087/1‑149 (MQ=255)                                      aaGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTAGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCt                                                                                                             >  5:63125/1‑148 (MQ=255)                                         gCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAgg                                                                                                         >  5:57659/1‑149 (MQ=255)                                              ccGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                                    >  3:169837/1‑149 (MQ=255)                                                gCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTAcc                                                                                                  <  4:113145/149‑1 (MQ=255)                                                    cTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCa                                                                                               >  5:99804/1‑148 (MQ=255)                                                         cTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCt                                                                                         <  6:90901/149‑1 (MQ=255)                                                          ttCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCtt                                                                                        >  4:251455/1‑149 (MQ=255)                                                               gctgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcg                                                                                   >  8:282126/1‑149 (MQ=255)                                                                 tgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                                    >  3:35232/1‑130 (MQ=255)                                                                 tgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGCGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                                    <  4:35232/130‑1 (MQ=255)                                                                   ctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCTCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcgcg                                                                                 <  1:289243/147‑1 (MQ=255)                                                                      ctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgcc                                                                             <  7:282126/148‑1 (MQ=255)                                                                           cacccaccTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgccgc                                                                           >  7:269312/1‑145 (MQ=255)                                                                           cacccaccTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgccgc                                                                           <  8:269312/145‑1 (MQ=255)                                                                                  cTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGc                                                                                          >  1:36263/1‑123 (MQ=255)                                                                                  cTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGc                                                                                          <  2:36263/123‑1 (MQ=255)                                                                                     cgcCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTgcg                                                             <  8:194372/149‑1 (MQ=255)                                                                                        cGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGa                                                          >  3:5210/1‑149 (MQ=255)                                                                                         gATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                         <  6:57659/149‑1 (MQ=255)                                                                                             gACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTg                                                      >  4:140684/1‑148 (MQ=255)                                                                                              aCGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGc                                                     <  3:251455/148‑1 (MQ=255)                                                                                              aCGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCa                                                    >  5:163045/1‑149 (MQ=255)                                                                                               cGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCAc                                                   >  7:196062/1‑149 (MQ=255)                                                                                                gCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACg                                                  <  8:196062/149‑1 (MQ=255)                                                                                                 cAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACgg                                                 <  8:205325/149‑1 (MQ=255)                                                                                                      aaGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCAtt                                             <  4:169837/148‑1 (MQ=255)                                                                                                       aGCCCCCGGCTCGGGCAGTAGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTgg                                           <  6:63125/149‑1 (MQ=255)                                                                                                         cccccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGcc                                         >  8:173285/1‑149 (MQ=255)                                                                                                          ccccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCt                                        <  6:163045/149‑1 (MQ=255)                                                                                                            ccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTg                                       >  6:269167/1‑148 (MQ=255)                                                                                                            ccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGc                                      <  4:5210/149‑1 (MQ=255)                                                                                                                cTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCa                                  <  4:305087/149‑1 (MQ=255)                                                                                                                  cGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCAcc                                   <  4:119704/146‑1 (MQ=255)                                                                                                                  cGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCAcc                                   >  3:119704/1‑146 (MQ=255)                                                                                                                       aGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTTGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGcgcc                           >  1:209565/1‑149 (MQ=255)                                                                                                                        gTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGcc                                                                  <  6:53708/109‑1 (MQ=255)                                                                                                                        gTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGcc                                                                  >  5:53708/1‑109 (MQ=255)                                                                                                                            gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcgcgc                                                                                >  6:310072/1‑91 (MQ=255)                                                                                                                            gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTgc                      >  2:290263/1‑149 (MQ=255)                                                                                                                               cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                         <  2:253438/107‑1 (MQ=38)                                                                                                                               cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                         >  1:253438/5‑111 (MQ=38)                                                                                                                               cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                         <  4:47038/107‑1 (MQ=38)                                                                                                                               cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                         >  3:47038/5‑111 (MQ=38)                                                                                                                                gACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTgccgcc                  >  4:293463/1‑149 (MQ=255)                                                                                                                                       gcagcaGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACt           <  6:99804/149‑1 (MQ=255)                                                                                                                                         agcagGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGa         <  3:293463/149‑1 (MQ=255)                                                                                                                                              gccgccCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTAc    >  6:274694/1‑149 (MQ=255)                                                                                                                                                cgccCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACgg  <  7:173285/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGG  >  NC_002947/4388857‑4389148

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