Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,631	T→C	8.0%	D467G (GAC→GGC) ‡	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,631	0	T	C	8.0%	209.5 / 6.3	75	D467G (GAC→GGC) ‡	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (4/2);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCG  >  NC_002947/1383486‑1383771                                                                                                                                                           |                                                                                                                                              tGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGcgtcgt                                                                                                                                            >  8:362520/1‑149 (MQ=255)  gggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGt‑t                                                                                                                                          <  6:108798/149‑1 (MQ=255)   ggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCTGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑Tc                                                                                                                                         <  1:330574/149‑1 (MQ=255)        gTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCATCATCTGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCa                                                                                                                                       <  7:226189/148‑1 (MQ=38)        gTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAgat                                                                                                                                    >  2:203157/1‑149 (MQ=255)        gTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAg                                                                                                                                      <  6:411379/147‑1 (MQ=255)             cTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGc                                                                                                                               >  2:412616/1‑149 (MQ=255)                gTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCt                                                                                                                            >  4:365181/1‑149 (MQ=255)                 tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCtt                                                                                                                           <  6:255752/149‑1 (MQ=255)                 tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCtt                                                                                                                           >  6:154018/1‑149 (MQ=255)                  cgcgGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTg                                                                                                                          >  8:400999/1‑149 (MQ=255)                     ggTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTgtcg                                                                                                                       <  8:262886/149‑1 (MQ=255)                       tcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTgtcgt                                                                                                                      >  7:385266/1‑148 (MQ=255)                           ccAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAg                                                                                                                 >  1:273996/1‑149 (MQ=255)                            cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAg                                                                                                                 <  7:400999/148‑1 (MQ=255)                            cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGc                                                                                                                >  3:95345/1‑149 (MQ=255)                             agcagcACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGct                                                                                                               <  3:389909/149‑1 (MQ=255)                                  caCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCt                                                                                                          <  7:362520/149‑1 (MQ=255)                                       cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGa                                                                                                       <  4:30746/147‑1 (MQ=255)                                       cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAt                                                                                                      <  4:114450/148‑1 (MQ=255)                                       cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATg                                                                                                     >  2:159421/1‑149 (MQ=255)                                        gtAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTACCTAGCGCAGTGGGACCAGTGCA‑‑ATATTGAAAGGAGCC‑‑‑‑‑‑‑ACGCCGTCGT‑TCCAgatg                                                                                                                                   >  7:172845/2‑118 (MQ=255)                                          aGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATgg                                                                                                    >  4:155798/1‑147 (MQ=255)                                            ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGc                                                                                                >  4:415803/1‑149 (MQ=255)                                            ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGc                                                                                                >  5:12206/1‑149 (MQ=255)                                              ggcggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTg                                                                                              >  2:188419/1‑149 (MQ=255)                                               gcggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGc                                                                                             <  3:365181/149‑1 (MQ=255)                                                cggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCa                                                                                            <  5:154018/149‑1 (MQ=255)                                                 ggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCaa                                                                                           <  5:59292/149‑1 (MQ=255)                                                   cagtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAAc                                                                                          >  1:390613/1‑148 (MQ=255)                                                      tcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACcaca                                                                                      >  3:93225/1‑149 (MQ=255)                                                         gtTCGGCTTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACgg                                                                                   >  6:277752/1‑149 (MQ=255)                                                          ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACggg                                                                                  <  1:203157/149‑1 (MQ=255)                                                                aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCGGTATAAAGAAGGCCGTGTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGAAGCTGCAACCACACg                                                                                    >  1:141527/1‑149 (MQ=255)                                                                  aCAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGcc                                                                           >  3:388164/1‑148 (MQ=255)                                                                    aGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCACTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAAt                                                                        >  6:349523/1‑149 (MQ=255)                                                                         gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGcc                                                                     <  8:385266/147‑1 (MQ=255)                                                                          tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                  <  2:390613/149‑1 (MQ=255)                                                                             gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGc                                                                >  8:155664/1‑148 (MQ=255)                                                                                gCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAAc                                                             <  1:188419/148‑1 (MQ=255)                                                                                    cgcgGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCAga                                                                                                                                     <  1:148702/72‑1 (MQ=21)                                                                                    cgcgGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCAga                                                                                                                                     >  2:148702/1‑72 (MQ=21)                                                                                     gcgGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATc                                                        >  1:272684/1‑148 (MQ=255)                                                                                       ggACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGAc                                                     <  6:12206/149‑1 (MQ=255)                                                                                        gACCACGTCGGTTTCGTAGCGAAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                      >  8:158647/1‑147 (MQ=255)                                                                                              gTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGa                                              >  4:118372/1‑149 (MQ=255)                                                                                              gTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGa                                              >  4:142268/1‑149 (MQ=255)                                                                                                 ggTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGAcc                                            <  2:273996/148‑1 (MQ=255)                                                                                                  gTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCgg                                          >  2:291449/1‑149 (MQ=255)                                                                                                   tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGc                                         >  1:70575/1‑149 (MQ=255)                                                                                                      cGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCtt                                       <  3:142268/148‑1 (MQ=255)                                                                                                      cGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCtt                                       <  1:412616/148‑1 (MQ=255)                                                                                                          gcgcAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAAGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCCACAGTAAGACCGGCATGCCGt                                  >  5:310628/1‑149 (MQ=255)                                                                                                            gcAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtc                                <  5:277752/149‑1 (MQ=255)                                                                                                              aGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtctc                              <  3:155798/149‑1 (MQ=255)                                                                                                               gTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCg                             <  7:253241/149‑1 (MQ=255)                                                                                                                tGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCAGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGc                            >  6:363674/1‑149 (MQ=255)                                                                                                                 gggACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGcc                           >  5:177715/1‑149 (MQ=255)                                                                                                                    aCCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCgg                                          <  8:224132/131‑1 (MQ=255)                                                                                                                    aCCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCgg                                          >  7:224132/1‑131 (MQ=255)                                                                                                                    aCCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAg                        >  8:224708/1‑149 (MQ=255)                                                                                                                      cAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTTTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCa                      >  4:382766/1‑149 (MQ=255)                                                                                                                        gggcatcaTCTGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCAGAT‑‑TGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCAcc                    >  6:12799/8‑149 (MQ=37)                                                                                                                        gggCA‑‑ATCTTGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCAcc                    <  7:266683/149‑1 (MQ=255)                                                                                                                        gggCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGc                                     <  8:237676/132‑1 (MQ=255)                                                                                                                        gggCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGc                                     >  7:237676/1‑132 (MQ=255)                                                                                                                         ggCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCt                   <  2:70575/149‑1 (MQ=255)                                                                                                                            a‑‑aTCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCCTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTgg                 >  1:306362/1‑148 (MQ=255)                                                                                                                                  tGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACtt            >  6:111313/1‑149 (MQ=255)                                                                                                                                  tGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCAGAT‑‑TGCCCTGGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACtt            >  6:380798/1‑147 (MQ=38)                                                                                                                                    gAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCg                             >  3:203284/1‑130 (MQ=255)                                                                                                                                    gAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCg                             <  4:203284/130‑1 (MQ=255)                                                                                                                                     aaCGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAg         >  1:27437/1‑149 (MQ=255)                                                                                                                                     aaCGGAGCC‑‑‑‑‑‑‑AGGGCGTCGA‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAg         <  1:21481/149‑1 (MQ=255)                                                                                                                                      aCGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAgc        <  4:388164/149‑1 (MQ=255)                                                                                                                                            cc‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCACCCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCg  >  4:48488/1‑149 (MQ=255)                                                                                                                                                           |                                                                                                                                              TGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCG  >  NC_002947/1383486‑1383771

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