Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,766,716	G→C	11.2%	R50P (CGT→CCT) ‡	PP_3329 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,766,716	0	G	C	11.2%	173.3 / 10.2	61	R50P (CGT→CCT) ‡	PP_3329	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (3/4);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATGATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATGG  >  NC_002947/3766571‑3766862                                                                                                                                                   |                                                                                                                                                   ttCATGATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTAt                                                                                                                                                 >  5:16905/1‑149 (MQ=255)     tGATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGt                                                                                                                                             <  7:205277/149‑1 (MQ=255)     tGATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGt                                                                                                                                             <  3:360586/149‑1 (MQ=255)      gATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGCTCAGCTTCCCTGCACCTATCGGTg                                                                                                                                            >  7:387410/1‑149 (MQ=255)       aTAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGc                                                                                                                                           <  1:10707/149‑1 (MQ=255)        tAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTg                                                                                                                                            >  6:166250/1‑147 (MQ=255)           cTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGCCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGAt                                                                                                                                       >  3:235830/1‑149 (MQ=255)              gCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAggca                                                                                                                                      <  3:327523/147‑4 (MQ=255)              gCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAggca                                                                                                                                      <  5:175104/147‑4 (MQ=255)              gCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAg                                                                                                                                         >  4:116354/1‑144 (MQ=255)              gCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAg                                                                                                                                         <  3:116354/144‑1 (MQ=255)                  gCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGtt                                                                                                                                 <  5:72905/148‑1 (MQ=255)                            gCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAggc                                                                                                                                       >  3:89045/1‑130 (MQ=255)                            gCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACCGATACGTGCAggc                                                                                                                                       <  4:89045/132‑3 (MQ=255)                              tACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTgg                                                                                                                    >  6:164215/1‑149 (MQ=255)                                cgcgAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTTACTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTgg                                                                                                                    <  1:135053/148‑1 (MQ=255)                                cgcgAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCGCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTgg                                                                                                                    <  8:313630/147‑1 (MQ=255)                                          cAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTTTCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCt                                                                                                        >  8:297169/1‑149 (MQ=255)                                          cAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGc                                                                                                                                     >  8:364114/1‑120 (MQ=255)                                          cAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGc                                                                                                                                     <  7:364114/120‑1 (MQ=255)                                                 ctctGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAACTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCAt                                                                                                 >  6:161475/1‑149 (MQ=255)                                                        tCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGcc                                                                                          >  4:79733/1‑149 (MQ=255)                                                           gagCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGc                                                                                               <  7:129315/141‑1 (MQ=255)                                                           gagCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGc                                                                                               >  8:129315/1‑141 (MQ=255)                                                                         ctcCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCg                                                                           <  1:124407/148‑1 (MQ=255)                                                                            ccTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCa                                                                       <  5:166250/149‑1 (MQ=255)                                                                             cTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAg                                                                      <  3:79733/149‑1 (MQ=255)                                                                                         gtTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGtctc                                                          >  3:123863/1‑149 (MQ=255)                                                                                            aCCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCgg                                                        <  4:123863/148‑1 (MQ=255)                                                                                              cgcgCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGtct                                                           >  7:23846/1‑143 (MQ=255)                                                                                              cgcgCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGcc                                                      >  8:369630/1‑148 (MQ=255)                                                                                              cgcgCTACTGGCCAGGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGtct                                                           <  8:23846/143‑1 (MQ=255)                                                                                               gagcTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCt                                                     >  7:164528/3‑148 (MQ=255)                                                                                                cgcTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGtt                                                                                                                                 >  6:325780/1‑71 (MQ=255)                                                                                                cgcTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGtt                                                                                                                                 <  5:325780/71‑1 (MQ=255)                                                                                                 gcTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACa                                                  <  2:51071/149‑1 (MQ=255)                                                                                                   tACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGt                                                >  5:192515/1‑149 (MQ=255)                                                                                                     cTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTgc                                              >  3:198064/1‑149 (MQ=255)                                                                                                      tGGCCATTGTGTTGACGGATGCATTGCCGGTCAGCTTGCCTGAACGTATCGGTGCTGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTgcc                                             <  8:387410/149‑1 (MQ=255)                                                                                                       ggCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTgccg                                            >  5:379693/1‑149 (MQ=255)                                                                                                        gCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTgccgc                                           <  8:164528/149‑1 (MQ=255)                                                                                                        gCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTgccgc                                           >  5:267665/1‑149 (MQ=255)                                                                                                               tgtTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATAGGGCATCATCGCCTCCAGAATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGCCCTACAGTGCCGCCAGCAGt                                    >  7:338906/1‑149 (MQ=255)                                                                                                                 ttCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTgg                                  >  6:343593/1‑149 (MQ=255)                                                                                                                  tCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGc                                 >  8:386633/1‑149 (MQ=255)                                                                                                                   cAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGc                                 <  7:369630/148‑1 (MQ=255)                                                                                                                   cAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGc                                 <  6:16905/148‑1 (MQ=255)                                                                                                                   cAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGc                                 <  5:343593/148‑1 (MQ=255)                                                                                                                      ggATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTaa                             >  5:225038/1‑149 (MQ=255)                                                                                                                       caTGCAGTGCAGGGGAGCTTTCCTGCACGTTTCGGTGCAGATGCCGTTGACAAGTTATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAAc                            <  4:235830/148‑1 (MQ=255)                                                                                                                              tGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCcgcgc                     >  8:296843/1‑149 (MQ=255)                                                                                                                              tGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCcgcgc                     >  1:406830/1‑149 (MQ=255)                                                                                                                                    tCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGTCACt               >  1:285239/1‑149 (MQ=255)                                                                                                                                        cTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTcat           <  5:164215/149‑3 (MQ=255)                                                                                                                                         ttGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCCACAGCAGTGGCGTAACCCCGCGAGCCACTTCGtt          >  8:261421/1‑149 (MQ=255)                                                                                                                                           gCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGTCACTTCGtt          <  2:285239/147‑1 (MQ=255)                                                                                                                                           gCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTg         <  6:225038/148‑1 (MQ=255)                                                                                                                                                cACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTGCAGCATGCTTGCCGTCACGACGTTTTCGCTCTGTACAATGGTCTCGGCCTACAGTGCCGGCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATg   <  6:192515/149‑1 (MQ=255)                                                                                                                                                cACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATg   >  1:249982/1‑149 (MQ=255)                                                                                                                                                 aCGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATgg  <  5:161475/149‑1 (MQ=255)                                                                                                                                                 aCGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATgg  <  2:406830/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TTCATGATAGCTCGCTGGCGTTGGTTGGCTACGCGAATTGCCAAGCGTCTCTGGTTCAGAGCAACGCTTTT‑CTCCCTGTTGGGGGTGGTTACCGCGCTACTGGCCATGGTGTTCAAGGATGCAGTGCCGGTCAGCTTGCCTGCACGTATCGGTGCAGATGCCGTTGACAAGATATTGGGCATCATCGCCTCCAGCATGCTTGCCGTCACGACGTTTTCGCTCAGTACAATGGTCTCGGCCTACAGTGCCGCCAGCAGTGGCGTAACCCCGCGCGCCACTTCGTTGGTCATGG  >  NC_002947/3766571‑3766862

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