Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,446,491	G→C	6.3%	Y358* (TAC→TAG)	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,491	0	G	C	6.3%	204.5 / 3.7	63	Y358* (TAC→TAG)	nicC	6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base C (2/2);  total (31/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 = 8.06e-01

CCTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGG  >  NC_002947/4446344‑4446636                                                                                                                                                     |                                                                                                                                                  ccTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTg                                                                                                                                                   >  4:300286/1‑148 (MQ=255)         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                                                                                                                                          >  6:128613/1‑140 (MQ=255)                  aGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                          <  5:128613/140‑1 (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                                                                                                                          <  2:234292/139‑1 (MQ=255)                                   gCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATa                                                                                                                          >  1:234292/1‑139 (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                                                                                     >  6:177775/1‑140 (MQ=255)                                                                       gCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGc                                                                                     <  5:177775/140‑1 (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                                                                     >  8:25318/1‑132 (MQ=255)                                                                                               aGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACAGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                     <  7:25318/132‑1 (MQ=255)                                                                                                 ctgc‑tgCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACa                                                  <  4:280966/149‑1 (MQ=255)                                                                                                        cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGc                                                                                                            <  6:229364/85‑1 (MQ=255)                                                                                                        cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGc                                                                                                            >  5:229364/1‑85 (MQ=255)                                                                                                        cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCgg                                            >  6:199293/1‑149 (MQ=255)                                                                                                        cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                             <  6:280126/148‑1 (MQ=255)                                                                                                          cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTc                                                                                  >  4:45579/1‑109 (MQ=255)                                                                                                          cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTc                                                                                  <  3:45579/109‑1 (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)                                                                                                                                                     |                                                                                                                                                  CCTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGC‑TGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGG  >  NC_002947/4446344‑4446636

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