Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,446,494	G→C	6.2%	L357L (CTC→CTG)	nicC ←	6‑hydroxynicotinate 3‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,446,494	0	G	C	6.2%	213.4 / 3.9	65	L357L (CTC→CTG)	nicC	6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base C (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCA  >  NC_002947/4446352‑4446639                                                                                                                                                |                                                                                                                                                  ccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                             <  8:184859/149‑1 (MQ=255)    tAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAggaa                                                                                                                                          <  3:45692/149‑3 (MQ=255)     aGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGgtgt                                                                                                                                         <  2:318160/149‑1 (MQ=255)       cgcAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTg                                                                                                                                       <  4:216566/149‑1 (MQ=255)          aGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                             <  5:128613/140‑1 (MQ=255)          aGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                             >  6:128613/1‑140 (MQ=255)            gTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGtt                                                                                                                                  <  3:300286/149‑1 (MQ=255)            gTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGtt                                                                                                                                  <  2:80821/149‑1 (MQ=255)               aTGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGa                                                                                                                               <  3:297514/149‑1 (MQ=255)                     gCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTg                                                                                                                                 >  4:100255/1‑141 (MQ=255)                     gCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTg                                                                                                                                 <  3:100255/141‑1 (MQ=255)                           gCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATa                                                                                                                             >  1:234292/1‑139 (MQ=255)                           gCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATa                                                                                                                             <  2:234292/139‑1 (MQ=255)                             cATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCCTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCt                                                                                                                 >  6:109966/1‑149 (MQ=255)                              aTTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTg                                                                                                                >  6:80381/1‑149 (MQ=255)                                  gcgcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGAt                                                                                                            >  2:272468/1‑149 (MQ=255)                                     cTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAgga                                                                                                                                           <  7:289046/115‑2 (MQ=255)                                     cTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAgga                                                                                                                                           >  8:289046/1‑114 (MQ=255)                                      tGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGccc                                                                                                        >  1:300331/1‑149 (MQ=255)                                        aCGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCATGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCg                                                                                                       >  1:241758/1‑149 (MQ=255)                                                 ccaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCg                                                                                              >  3:160186/1‑148 (MQ=255)                                                 ccaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCg                                                                                              <  4:160186/148‑1 (MQ=255)                                                    gCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTg                                                                                          >  7:166301/1‑149 (MQ=255)                                                    gCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGACCCAGGTCTTGGCGTTCGATACCGATTGCACCTGCGATGCCCGCTCTATGCGATTg                                                                                          >  4:206093/1‑149 (MQ=255)                                                       acgaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCt                                                                                       >  2:209971/1‑149 (MQ=255)                                                       acgaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCCTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCt                                                                                       <  5:109966/149‑1 (MQ=255)                                                         gaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                          <  5:288585/128‑1 (MQ=255)                                                         gaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                          >  6:288585/1‑128 (MQ=255)                                                           aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCAt                                                                                   <  5:80381/149‑1 (MQ=255)                                                           aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCAt                                                                                   <  1:209971/149‑1 (MQ=255)                                                            ggTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATa                                                                                  >  5:280126/1‑149 (MQ=255)                                                               gCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGc                                                                                        <  5:177775/140‑1 (MQ=255)                                                               gCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGc                                                                                        >  6:177775/1‑140 (MQ=255)                                                               gCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACaa                                                                               >  4:280873/1‑149 (MQ=255)                                                                   aTGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAAcgc                                                                            <  2:241758/148‑1 (MQ=255)                                                                    tGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGtt                                                                                           >  4:190304/1‑132 (MQ=255)                                                                    tGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGtt                                                                                           <  3:190304/132‑1 (MQ=255)                                                                      ccgccTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAAACCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                        >  2:238758/1‑149 (MQ=255)                                                                        gccTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTCCGCTTGGCTTCATACAACGCGAACGc                                                                      >  7:79675/1‑149 (MQ=255)                                                                             cccGCTTTAAAGCTCC‑TGCCCATAGAGGACGTAGACGTAAACACAGGCTGGATCTTCCTTGCTGTAGAGCCAGGTGATGTCGTTGGATACCGAATGCACCTGCGATGCCATCTCTTTGCGCTTGGCTTCATACAACGCGAACGCCGTGc                                                                 >  2:314454/1‑149 (MQ=255)                                                                              ccGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCg                                                                >  5:91625/1‑149 (MQ=255)                                                                              ccGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGGGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCg                                                                >  7:67153/1‑149 (MQ=255)                                                                                 cTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggtg                                                             <  8:166301/149‑1 (MQ=255)                                                                                    tCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggtggt                                                           <  6:286098/148‑1 (MQ=255)                                                                                       aGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACAGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                        <  7:25318/132‑1 (MQ=255)                                                                                       aGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACAGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                        >  8:25318/1‑132 (MQ=255)                                                                                         ctgc‑tgCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACa                                                     <  4:280966/149‑1 (MQ=255)                                                                                                cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGc                                                                                                               >  5:229364/1‑85 (MQ=255)                                                                                                cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGc                                                                                                               <  6:229364/85‑1 (MQ=255)                                                                                                cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCgg                                               >  6:199293/1‑149 (MQ=255)                                                                                                cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                                <  6:280126/148‑1 (MQ=255)                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTc                                                                                     <  3:45579/109‑1 (MQ=255)                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTc                                                                                     >  4:45579/1‑109 (MQ=255)                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCTGACAGCCCGGtc                                             <  1:314454/149‑1 (MQ=255)                                                                                                      gagGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCt                                         >  1:82136/1‑149 (MQ=255)                                                                                                                 cGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCgg                              <  2:82136/149‑1 (MQ=255)                                                                                                                         cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGGGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCAt                      <  8:67153/149‑1 (MQ=255)                                                                                                                           ggCTGGATCTTCCTGGCTGTAGCGCCAGGTGTTGGCGTTGGATCCCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACATCCCGGTCTCCTGCAGGCAGCGGGTCAGCATgg                    <  3:206093/149‑1 (MQ=255)                                                                                                                                 caacaccTGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGt              >  5:136460/6‑149 (MQ=255)                                                                                                                                  tCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTACGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACCCGAACGACGTGCGGTGGTCGGACAGACCGGTCTCCTGCAGGCACAGTGTCACCATCGCcgagac             >  4:81178/1‑145 (MQ=255)                                                                                                                                       cTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGa        <  2:300331/149‑1 (MQ=255)                                                                                                                                          gCTGTTGAGCCAGGTGTTGGCGTTGGATTGCGATTTCACCTGCGATGCGCGCTGTTTGCTGTTGGCTTCATACAACGCGTACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCTGCGGGTCAGCATGGCCGCGTCCTCGATg      <  3:81178/148‑1 (MQ=255)                                                                                                                                          gCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATgg     >  4:11977/1‑149 (MQ=255)                                                                                                                                             ttAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCGCGCGCTTTGCGGTTGGCTTCATACAACGCGAACGCCGGGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCa  <  5:319864/148‑1 (MQ=255)                                                                                                                                              tAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGcc   <  1:272468/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  CCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCA  >  NC_002947/4446352‑4446639

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