Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,631	T→C	6.5%	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	6.5%	171.9 / 2.6	61	D467G (GAC→GGC) ‡	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
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.96e-01

TGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTA  >  NC_002947/1383486‑1383776                                                                                                                                                  |                                                                                                                                                  tGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGcgtc                                                                                                                                                  >  8:301192/1‑147 (MQ=255) tGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGcgtc                                                                                                                                                  <  7:301192/147‑1 (MQ=255)    gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTcc                                                                                                                                             <  4:332631/149‑1 (MQ=255)          tttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAgatga                                                                                                                                       >  4:276898/1‑149 (MQ=255)              ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGcc                                                                                                                                   >  8:87283/1‑149 (MQ=255)                gTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGGCGGTTGCGTAGCGCAGTGGGACCAGAGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                                  <  5:239849/148‑1 (MQ=255)                 tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCtt                                                                                                                                <  3:59483/149‑1 (MQ=255)                   gcgGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgt                                                                                                                              >  4:89542/1‑149 (MQ=255)                    cgGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtc                                                                                                                             <  5:264249/149‑1 (MQ=255)                       tcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtcgtc                                                                                                                          <  3:109189/149‑1 (MQ=255)                          tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCa                                                                                                                         >  6:148857/1‑147 (MQ=255)                           ccAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                                      <  3:276898/149‑1 (MQ=255)                            cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                                      <  2:120852/148‑1 (MQ=255)                            cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGACGCCCTGGCTCCGTTCCAga                                                                                                                                          <  3:21276/128‑1 (MQ=39)                            cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGACGCCCTGGCTCCGTTCCAga                                                                                                                                          >  4:21276/1‑128 (MQ=39)                                   aCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAg                                                                                                                                           >  4:298657/1‑120 (MQ=255)                                   aCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAg                                                                                                                                           <  3:298657/120‑1 (MQ=255)                                       cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCACATGATGACCTTGTCGGCATAGCTCTGCTCGATg                                                                                                          >  6:110395/1‑149 (MQ=255)                                        aTAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAt                                                                                                           >  3:103661/1‑147 (MQ=255)                                        aTAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAt                                                                                                           <  4:103661/147‑1 (MQ=255)                                           gCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAg                                                                                                      <  7:87283/149‑1 (MQ=255)                                            ccGGCCGCAGTCAGTTCTGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGc                                                                                                     <  8:294437/149‑1 (MQ=255)                                              ggcggcAGTCAGTTCGGCATACTGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTg                                                                                                   <  6:10939/149‑1 (MQ=255)                                                  gcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCGTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAAc                                                                                               <  6:333964/149‑1 (MQ=255)                                                   cagtcagtTCGGCATACAGCTTGTCGGTTGCTTCGAGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCaa                                                                                                <  4:167046/147‑1 (MQ=255)                                                     gtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACca                                                                                             >  6:91267/1‑148 (MQ=255)                                                         gtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACgg                                                                                        >  8:276442/1‑149 (MQ=255)                                                         gtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACgg                                                                                        >  4:275327/1‑149 (MQ=255)                                                          ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAg                                                                                                      >  7:298188/1‑134 (MQ=255)                                                          ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAg                                                                                                      <  8:298188/134‑1 (MQ=255)                                                                aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGc                                                                                                        <  7:226340/126‑1 (MQ=255)                                                                aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGc                                                                                                        >  8:226340/1‑126 (MQ=255)                                                                aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                                 >  3:75717/1‑149 (MQ=255)                                                                aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                                 >  7:87583/1‑149 (MQ=255)                                                                   cAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCaa                                                                              <  5:336335/149‑1 (MQ=255)                                                                   cAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCaa                                                                              <  5:91267/149‑1 (MQ=255)                                                                      cTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGc                                                                           <  3:89542/149‑1 (MQ=255)                                                                         gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCg                                                                         <  7:276442/148‑1 (MQ=255)                                                                          tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCATGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                       <  5:148857/149‑1 (MQ=255)                                                                            ggTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAg                                                                      <  2:161018/148‑1 (MQ=255)                                                                             gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCa                                                                    >  7:293251/1‑149 (MQ=255)                                                                             gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGACGTAGc                                                                     >  7:321209/1‑148 (MQ=255)                                                                                 cTTCGCGGACCACGTCGGTTTCGTAGCCCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTGCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATCGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAAca                                                                 >  3:22144/1‑147 (MQ=255)                                                                                       ggACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                           >  3:320191/1‑148 (MQ=255)                                                                                               tCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGc                                                                                                        >  3:224979/1‑95 (MQ=255)                                                                                               tCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGc                                                                                                        <  4:224979/95‑1 (MQ=255)                                                                                                        tAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGcc                                         <  4:320191/149‑1 (MQ=255)                                                                                                               gTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCGATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCg                                  >  5:156369/1‑149 (MQ=255)                                                                                                                   gACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATGGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                       <  4:247410/108‑1 (MQ=255)                                                                                                                   gACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATGGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                       >  3:247410/1‑108 (MQ=255)                                                                                                                        gggCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                       >  7:1145/1‑103 (MQ=255)                                                                                                                        gggCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                       <  8:1145/103‑1 (MQ=255)                                                                                                                         ggCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCt                        <  8:321209/149‑1 (MQ=255)                                                                                                                          gCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTg                       >  3:308867/1‑149 (MQ=255)                                                                                                                           cAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTgg                      <  8:87583/149‑1 (MQ=255)                                                                                                                              tCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCAc                   >  5:74335/1‑149 (MQ=255)                                                                                                                                tGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACtt                 >  5:270754/1‑149 (MQ=255)                                                                                                                                tGGAACGACGCCCTGGCTCCGTTCCAGAT‑‑TGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGtt                                      >  3:231812/1‑126 (MQ=37)                                                                                                                                tGGAACGACGCCCTGGCTCCGTTCCAGAT‑‑TGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGtt                                      <  4:231812/126‑1 (MQ=37)                                                                                                                                     cGGAGCCAGGGCGTCGTTCAAGATGATGCCCTTGTCGTCATAGCACTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAgcg            >  8:276736/1‑149 (MQ=255)                                                                                                                                          ccAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGTACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCg       >  8:316016/1‑149 (MQ=255)                                                                                                                                          ccAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCg       <  8:293251/149‑1 (MQ=255)                                                                                                                                               gCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTa  <  5:110395/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTA  >  NC_002947/1383486‑1383776

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