Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,986,052	C→G	6.7%	G783R (GGC→CGC)	PP_2612 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,986,052	0	C	G	6.7%	184.1 / 3.6	60	G783R (GGC→CGC)	PP_2612	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CATCGTTTTCCCTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGCCTTAAA  >  NC_002947/2985908‑2986199                                                                                                                                                 |                                                                                                                                                     cATCGTTTTCCCTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGc                                                                                                                                                   >  2:24103/1‑148 (MQ=255)        ttCCCTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGa                                                                                                                                            <  7:92201/148‑1 (MQ=255)          cccTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAAt                                                                                                                                          <  1:24103/148‑1 (MQ=255)            cTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAg                                                                                                                                               >  3:328566/1‑141 (MQ=255)            cTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAg                                                                                                                                               <  4:328566/141‑1 (MQ=255)            cTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAg                                                                                                                                       <  4:288840/149‑1 (MQ=255)            cTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAg                                                                                                                                       <  2:236950/149‑1 (MQ=255)             ttGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAg                                                                                                                                       >  3:308641/1‑148 (MQ=255)                  ccGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGt                                                                                                                                 >  8:102967/1‑149 (MQ=255)                   cGAATATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAg                                                                                                                                               <  3:257397/134‑1 (MQ=255)                   cGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTc                                                                                                                                >  6:333288/1‑149 (MQ=255)                   cGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATAAAACTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCATCCTGCAAAg                                                                                                                                               >  4:257397/1‑134 (MQ=255)                          gCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑g                                                                                                                        >  6:102522/1‑149 (MQ=255)                                 ctgGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGAt                                                                                                                    <  6:139096/146‑1 (MQ=255)                                       cATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGct                                                                                                           >  2:248548/1‑149 (MQ=255)                                       cATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGct                                                                                                           >  7:112886/1‑149 (MQ=255)                                       cATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTACAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGct                                                                                                           >  7:294643/1‑149 (MQ=255)                                              cTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGcttccttc                                                                                                     <  5:52271/148‑1 (MQ=255)                                              cTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGcttccttc                                                                                                     >  1:168001/1‑148 (MQ=255)                                               tAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGcttccttc                                                                                                     >  6:304640/1‑147 (MQ=255)                                               tAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGcttccttc                                                                                                     >  5:208719/1‑147 (MQ=255)                                               tAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGcttccttc                                                                                                     >  5:125414/1‑147 (MQ=255)                                               tAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTa                                                                                                   <  5:102522/149‑1 (MQ=255)                                                aaaGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTa                                                                                                   <  2:175013/148‑1 (MQ=255)                                                 aaGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAg                                                                                                  <  2:328963/148‑1 (MQ=255)                                                   gggCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGaa                                                                                                >  6:269264/1‑148 (MQ=255)                                                             gCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGaa                                                                                     >  8:268643/1‑149 (MQ=255)                                                              ccGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAAc                                                                                    <  3:262785/149‑1 (MQ=255)                                                               cGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAAcc                                                                                   <  6:208719/149‑1 (MQ=255)                                                                     tGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCGTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAg                                                                              >  5:287189/1‑148 (MQ=255)                                                                     tGAGCTTTTGAGTAGAGGAATAGACATCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCGTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAg                                                                              >  5:111389/1‑148 (MQ=255)                                                                      gAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCt                                                                                                    <  2:200706/125‑1 (MQ=255)                                                                          ttttGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGACCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCGAGAACAAGAATCGAACCTAAAGCagtag                                                                        >  8:98263/1‑149 (MQ=255)                                                                           tttGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCTCGCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCa                                                                            <  1:69983/145‑1 (MQ=255)                                                                           tttGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCa                                                                            >  2:69983/1‑144 (MQ=255)                                                                                      gAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAg                                                            >  6:185276/1‑149 (MQ=255)                                                                                           gACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCt                                                       <  4:308641/149‑1 (MQ=255)                                                                                             ccTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTg                                                      <  7:102967/148‑1 (MQ=255)                                                                                              cTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGta                                                    >  4:87899/1‑149 (MQ=255)                                                                                              cTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCAGta                                                    >  4:83300/1‑149 (MQ=255)                                                                                               tCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGta                                                    >  1:152810/1‑148 (MQ=255)                                                                                                 cACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCGTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGtat                                                   <  6:287189/147‑1 (MQ=255)                                                                                                 cACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCGTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGtat                                                   <  6:111389/147‑1 (MQ=255)                                                                                                 cACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGtata                                                  >  2:188899/1‑148 (MQ=255)                                                                                                   ccGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATaa                                               <  5:304640/149‑1 (MQ=255)                                                                                                   ccGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATaa                                               <  5:185276/149‑1 (MQ=255)                                                                                                     gATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGc                                             >  4:323705/1‑149 (MQ=255)                                                                                                      aTCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCCGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGcc                                            >  6:90287/1‑149 (MQ=255)                                                                                                      aTCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGcc                                            >  4:126213/1‑149 (MQ=255)                                                                                                       tCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCt                                           <  1:248548/149‑1 (MQ=255)                                                                                                         aaaGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCTGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGa                                         <  6:125414/149‑1 (MQ=255)                                                                                                                 tCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCGGCAAAGGGAATGAGGACTGTCCTGCCT‑GCGATGGCCAGGCTGCCTTCGAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCCGCCGCGGTATATAAGCCTGAAATAAGAt                                 >  1:81462/1‑149 (MQ=255)                                                                                                                      tCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCt                                                       <  4:29430/122‑1 (MQ=255)                                                                                                                      tCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGcc                             <  2:152810/148‑1 (MQ=255)                                                                                                                      tCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCGAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCt                                                       >  3:29430/1‑122 (MQ=255)                                                                                                                                      tgctgTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTa            <  5:269264/149‑1 (MQ=255)                                                                                                                                          gTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGc        >  2:28344/1‑149 (MQ=255)                                                                                                                                           tCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCACTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGc        <  8:112886/148‑1 (MQ=255)                                                                                                                                            ctaGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGcc       <  3:186800/146‑1 (MQ=255)                                                                                                                                               gCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAATCCTGAAATAAGATTGCCGTTCTGCCATAGACGTATTgtc       >  4:186800/1‑143 (MQ=255)                                                                                                                                                ccTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGCCTTaaa  <  3:83300/149‑1 (MQ=255)                                                                                                                                                ccTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGCCTTaaa  <  3:87899/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CATCGTTTTCCCTTGTCCCGAAGATGCTTCTGCTGGCCCATAATTCTAAAGGGCTGTGAAGCCGGGAGTGAGCTTTTGAGTAGAGGAATAGACCTCCACCGATCAAAGTTGCTCCGATCAAGGCGATTCCAGCTGCTGTCCAGCCTGCAAAGGGAATGAGGACTGTCGTGCCT‑GCGATGGCAAGGCTTCCTTCTAGAACAAGAATCGAACCTAAAGCAGTAGTCGCGCTGGCAGCCGCTGTATATAAGCCTGAAATAAGATTGCCGTTCTGCCATTGACGTATTGCCTTAAA  >  NC_002947/2985908‑2986199

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