Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,621	A→C	6.2%	A470A (GCT→GCG)	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,621	0	A	C	6.2%	254.2 / 6.7	81	A470A (GCT→GCG)	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (41/35);  new base C (3/2);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

ACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGC  >  NC_002947/1383476‑1383765                                                                                                                                                    |                                                                                                                                                         aCTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGcc                                                                                                                                                       <  1:104720/149‑1 (MQ=255)   ttGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑Ag                                                                                                                                              >  4:389909/1‑149 (MQ=255)     gctctcTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑Aggg                                                                                                                                            <  5:274799/144‑1 (MQ=255)     gATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑Agg                                                                                                                                             <  5:230680/148‑1 (MQ=255)       tGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGcg                                                                                                                                          <  3:416830/149‑1 (MQ=255)           tGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑Aggg                                                                                                                                            >  6:274799/1‑143 (MQ=255)           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)                   ttttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                          >  5:83190/1‑128 (MQ=255)                   ttttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                          <  6:83190/128‑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:154018/1‑149 (MQ=255)                           tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCtt                                                                                                                     <  6:255752/149‑1 (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                                                                                                                               >  2:148702/1‑72 (MQ=21)                                                                                              cgcgGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCC‑‑‑‑‑‑‑CTGGCTCCGT‑TCCAga                                                                                                                               <  1:148702/72‑1 (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:142268/1‑149 (MQ=255)                                                                                                        gTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGa                                        >  4:118372/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                                 <  1:412616/148‑1 (MQ=255)                                                                                                                cGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCtt                                 <  3:142268/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                               >  7:237676/1‑132 (MQ=255)                                                                                                                                  gggCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGc                               <  8:237676/132‑1 (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)                                                                                                                                                    |                                                                                                                                                         ACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCC‑‑‑‑‑‑‑AGGGCGTCGT‑TCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGC  >  NC_002947/1383476‑1383765

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