Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,096,065	A→G	8.3%	T133A (ACC→GCC)	PP_5607 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,096,065	0	A	G	8.3%	165.3 / 4.6	60	T133A (ACC→GCC)	PP_5607	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (3/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGG  >  NC_002947/4095928‑4096210                                                                                                                                          |                                                                                                                                                   cGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAgc                                                                                                                                         <  5:425931/149‑1 (MQ=255)   ggCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAgcgc                                                                                                                                       <  1:258416/149‑1 (MQ=255)    gCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgag                                                                                                                                      >  3:7979/1‑147 (MQ=255)    gCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                       <  4:7979/148‑2 (MQ=255)     ccTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGccc                                                                                                                                     >  3:390453/1‑149 (MQ=255)         gtgtTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATg                                                                                                                                 >  7:396485/1‑149 (MQ=255)          tgtTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATgg                                                                                                                                <  1:265267/149‑1 (MQ=255)              gctgcATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCg                                                                                                                             <  4:389103/148‑1 (MQ=255)                  catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTgg                                                                                                                          <  5:188475/147‑1 (MQ=255)                  catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCt                                                                                                                        <  1:394376/149‑1 (MQ=255)                    tggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAg                                                                                                                      >  5:363527/1‑149 (MQ=255)                         tggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAAc                                                                                                                 <  4:323063/149‑1 (MQ=255)                                   gggCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGcc                                                                                                        <  8:97371/148‑1 (MQ=255)                                       aCCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTc                                                                                                   >  6:437798/1‑149 (MQ=255)                                         cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                 >  4:114670/1‑149 (MQ=255)                                         cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                 >  6:164589/1‑149 (MQ=255)                                         cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                 >  6:121616/1‑149 (MQ=255)                                           gCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccccc                                                                                               >  2:167494/1‑149 (MQ=255)                                                  ccccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                       >  6:381189/1‑101 (MQ=255)                                                  ccccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                       <  5:381189/102‑2 (MQ=255)                                                    ccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgcg                                                                                      >  3:58524/1‑149 (MQ=255)                                                         gACCAGGGGTGTGCACAGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTACTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGtt                                                                                 <  1:298918/149‑1 (MQ=255)                                                            cAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTc                                                                              <  5:437798/149‑1 (MQ=255)                                                                ggTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCTTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTg                                                                         >  1:292208/1‑149 (MQ=255)                                                                 gtgtGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTgg                                                                        <  6:378978/149‑1 (MQ=255)                                                                 gtgtGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTgg                                                                        <  7:360570/149‑1 (MQ=255)                                                                      cACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGc                                                                   <  8:41812/149‑1 (MQ=255)                                                                       aCCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCg                                                                  >  8:210966/1‑149 (MQ=255)                                                                        ccGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGc                                                                   >  5:193798/1‑147 (MQ=255)                                                                        ccGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGc                                                                   <  6:193798/147‑1 (MQ=255)                                                                           gTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCgtggt                                                              >  1:121646/1‑149 (MQ=255)                                                                              cAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGtt                                                           <  6:363527/149‑1 (MQ=255)                                                                                 ccATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCgtgg                                                               >  6:181349/1‑142 (MQ=255)                                                                                 ccATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCgtgg                                                               <  5:181349/142‑1 (MQ=255)                                                                                  cATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGc                                                        >  1:57023/1‑148 (MQ=255)                                                                                        ccTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGatcatc                                                 <  7:210966/149‑1 (MQ=255)                                                                                        ccTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGatc                                                    <  5:121616/146‑1 (MQ=255)                                                                                         cTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGAtcatca                                                >  5:275333/1‑149 (MQ=255)                                                                                                  aCTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGt                                       <  3:114670/149‑1 (MQ=255)                                                                                                   cTGCATAACCTGTTTGGTCTGCGGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGtt                                      >  7:271667/1‑149 (MQ=255)                                                                                                      cATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGAt                                   <  2:426907/149‑1 (MQ=255)                                                                                                       aTAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATg                                  <  2:292208/149‑1 (MQ=255)                                                                                                          aCCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGc                               <  2:57023/149‑1 (MQ=255)                                                                                                             tGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGa                                                      >  8:164548/1‑123 (MQ=255)                                                                                                             tGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGa                                                      <  7:164548/123‑1 (MQ=255)                                                                                                                  ggTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCACTTAGGTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcgg                        >  1:141234/1‑149 (MQ=255)                                                                                                                  ggTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcgg                        <  1:156210/148‑1 (MQ=255)                                                                                                                     ctgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                      >  5:182078/1‑147 (MQ=255)                                                                                                                      tgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTGTGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                      >  4:72788/1‑146 (MQ=255)                                                                                                                        ctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCaga                  <  5:95363/148‑1 (MQ=255)                                                                                                                           aTGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGt                <  4:58524/147‑1 (MQ=255)                                                                                                                           aTGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTg               <  6:182078/148‑1 (MQ=255)                                                                                                                             gTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGc             >  2:14172/1‑148 (MQ=255)                                                                                                                              ttCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGc             >  2:170081/1‑147 (MQ=255)                                                                                                                              ttCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCa           <  4:390453/149‑1 (MQ=255)                                                                                                                               atagggcGCTGGCCAAGGTCAGCGACCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGAAGAGTGGCCAg          >  6:367050/7‑149 (MQ=255)                                                                                                                                ccTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAg          <  6:275333/148‑1 (MQ=255)                                                                                                                                ccTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAg          <  5:164589/148‑1 (MQ=255)                                                                                                                                      gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTc    <  2:121646/148‑1 (MQ=255)                                                                                                                                       cTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCgg  <  1:167494/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                   CGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTC‑CCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGG  >  NC_002947/4095928‑4096210

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