Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,389,006	G→C	7.7%	P579A (CCG→GCG) ‡	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,389,006	0	G	C	7.7%	208.6 / 7.0	65	P579A (CCG→GCG) ‡	PP_3865	tail protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base C (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGGGCTTTC  >  NC_002947/4388863‑4389154                                                                                                                                                |                                                                                                                                                     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)                           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                                                                                 <  8:269312/145‑1 (MQ=255)                                                                     cacccaccTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgccgc                                                                                 >  7:269312/1‑145 (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                                                                        >  5:53708/1‑109 (MQ=255)                                                                                                                  gTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGcc                                                                        <  6:53708/109‑1 (MQ=255)                                                                                                                      gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcgcgc                                                                                      >  6:310072/1‑91 (MQ=255)                                                                                                                      gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTgc                            >  2:290263/1‑149 (MQ=255)                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                               >  1:253438/5‑111 (MQ=38)                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                               <  4:47038/107‑1 (MQ=38)                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                               >  3:47038/5‑111 (MQ=38)                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                               <  2:253438/107‑1 (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)                                                                                                                                            cccGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGa               >  2:190052/1‑140 (MQ=255)                                                                                                                                            cccGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGa               <  1:190052/140‑1 (MQ=255)                                                                                                                                              cGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGGGCtt    <  5:274694/149‑1 (MQ=255)                                                                                                                                                gTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGcc                                               <  2:148653/104‑1 (MQ=255)                                                                                                                                                gTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGcc                                               >  1:148653/1‑104 (MQ=255)                                                                                                                                                gTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGGGCTTTc  <  3:140684/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGAGTTACGGGCTTTC  >  NC_002947/4388863‑4389154

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