Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,641,817	T→A	7.9%	Q463L (CAG→CTG)	PP_3209 ←	isoxanthopterin deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,641,817	0	T	A	7.9%	227.1 / 7.0	76	Q463L (CAG→CTG)	PP_3209	isoxanthopterin deaminase
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base A (4/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAA  >  NC_002947/3641674‑3641962                                                                                                                                                |                                                                                                                                                   cAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCCCCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGcac                                                                                                                                               >  7:381977/1‑147 (MQ=255) cAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑Gt                                                                                                                                              >  7:13239/1‑149 (MQ=255) cAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGcac                                                                                                                                               >  7:376588/1‑147 (MQ=255)    tCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGcc                                                                                                                                           >  3:387682/1‑149 (MQ=255)      gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCt                                                                                                                                          <  5:426680/148‑1 (MQ=255)      gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTg                                                                                                                                         <  2:95533/149‑1 (MQ=255)       aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTgg                                                                                                                                        <  2:260402/149‑1 (MQ=255)       aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTgg                                                                                                                                        >  2:319443/1‑149 (MQ=255)          ggCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTg                                                                                                                                     <  5:114679/149‑1 (MQ=255)           gCAGGGGAGCAGGGCCAGGCGTGGCCAGTTCCACACCCAGCACGTCGTCTGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGC‑ACGCCAGGCTgt                                                                                                                                    <  8:381977/149‑2 (MQ=255)           gCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGc                                                                                                                                    >  5:275432/1‑149 (MQ=255)           gCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGC‑ACGCCAGGCTgt                                                                                                                                    <  8:376588/149‑2 (MQ=255)            cAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGcc                                                                                                                                   <  7:298888/149‑1 (MQ=255)            cAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGcc                                                                                                                                   <  6:79569/149‑1 (MQ=255)             aGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGccc                                                                                                                                  <  4:387682/149‑1 (MQ=255)              ggggAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGACTGCCCa                                                                                                                                 <  1:293300/149‑1 (MQ=255)                 gAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGGCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGtt                                                                                                                              <  5:346220/149‑1 (MQ=255)                   gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCg                                                                                                                            <  3:286320/149‑1 (MQ=255)                   gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCg                                                                                                                            <  5:140397/149‑1 (MQ=255)                    cAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCgg                                                                                                                           <  8:13239/149‑1 (MQ=255)                       gCCCCGGCGGGGCCAGTTCCACACCCAGCTCGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGt                                                                                                                               <  2:389941/142‑1 (MQ=255)                       gCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGt                                                                                                                               >  1:389941/1‑142 (MQ=255)                        ccccGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAg                                                                                                                       <  1:319443/149‑1 (MQ=255)                                   ccAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                            <  3:125505/149‑1 (MQ=255)                                         ccACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATc                                                                                                       <  5:333299/148‑1 (MQ=255)                                           acacCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCa                                                                                                             >  6:215844/1‑140 (MQ=255)                                           acacCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCa                                                                                                             <  5:215844/140‑1 (MQ=255)                                              cccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCg                                                                                                 >  3:390555/1‑149 (MQ=255)                                               ccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGt                                                                                                >  7:70048/1‑149 (MQ=255)                                                      gtcgtcGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCAc                                                                                         >  1:92766/1‑149 (MQ=255)                                                       tcgtcgGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACg                                                                                        >  8:27173/1‑149 (MQ=255)                                                          tcgGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCt                                                                                      >  5:369816/1‑148 (MQ=255)                                                              cTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggc                                                                                 >  3:240360/1‑149 (MQ=255)                                                               ttCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcg                                                                                >  1:374373/1‑149 (MQ=255)                                                                 cAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcgg                                                                               <  3:292020/148‑1 (MQ=255)                                                                   aGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgccg                                                                            >  5:71699/1‑149 (MQ=255)                                                                    ggCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgccgc                                                                           >  3:3274/1‑149 (MQ=255)                                                                      cgcAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCt                                                                         <  3:254496/149‑1 (MQ=255)                                                                       gcAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTa                                                                        >  1:322844/1‑149 (MQ=255)                                                                        cAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTAc                                                                       <  4:240360/149‑1 (MQ=255)                                                                         aGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTAcc                                                                      <  2:92766/149‑1 (MQ=255)                                                                          gggCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCa                                                                     <  6:275432/149‑1 (MQ=255)                                                                                   ccTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGccc                                                            <  2:322844/149‑1 (MQ=255)                                                                                    cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                           >  6:112605/1‑149 (MQ=255)                                                                                    cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                           >  6:128882/1‑149 (MQ=255)                                                                                       cAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGc                                                          <  4:3274/147‑1 (MQ=255)                                                                                        agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                            <  4:19978/96‑1 (MQ=255)                                                                                        agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                            >  3:19978/1‑96 (MQ=255)                                                                                        agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGc                                                       <  6:369816/149‑1 (MQ=255)                                                                                         gagGCGTCCAGCTCATGCCCGTACCCCCGCGCCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTCCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATACTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCg                                                      >  6:128258/1‑149 (MQ=255)                                                                                                ccAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAgg                                                                                                                      <  5:6577/78‑1 (MQ=255)                                                                                                ccAGCTCATGCGCGTACCACCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAgg                                                                                                                      >  6:6577/1‑78 (MQ=255)                                                                                                 cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                              >  6:332058/1‑149 (MQ=255)                                                                                                 cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                              >  8:20678/1‑149 (MQ=255)                                                                                                  aGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGc                                             >  5:183155/1‑149 (MQ=255)                                                                                                     tCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCAc                                          >  6:22677/1‑149 (MQ=255)                                                                                                     tCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCAc                                          >  7:247914/1‑149 (MQ=255)                                                                                                      cATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCAGTCCCTGGCTGCCCAGTTCGGCGAGGTCCAGGCCAGGAATCGCATCGTTTTCAACGCTCAGGCGGCCGCCTACCAGCAATCCCCCCAGCGTCGGACGGCCAc                                          >  1:89395/1‑149 (MQ=255)                                                                                                       aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                            >  6:315343/1‑81 (MQ=255)                                                                                                       aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtttt                                                                                             >  5:90216/1‑96 (MQ=255)                                                                                                       aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtttt                                                                                             <  6:90216/96‑1 (MQ=255)                                                                                                            cgTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAAc                                   <  6:183155/149‑1 (MQ=255)                                                                                                                 ccccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTg                              >  8:19408/1‑149 (MQ=255)                                                                                                                 ccccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGACGGACGGCCACCGCCAACCACTg                              >  4:231024/1‑149 (MQ=255)                                                                                                                 ccccGCGGCCTGATGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCAGCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTg                              <  5:128258/149‑1 (MQ=255)                                                                                                                   ccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTggg                            <  6:71699/149‑1 (MQ=255)                                                                                                                           tGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGcc                    >  3:404604/1‑149 (MQ=255)                                                                                                                             ttGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGccc                   <  4:390555/148‑1 (MQ=255)                                                                                                                             ttGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGcc                    <  7:20678/147‑1 (MQ=255)                                                                                                                                  ggCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCa             <  2:89395/149‑1 (MQ=255)                                                                                                                                    cGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATg           <  8:247914/149‑1 (MQ=255)                                                                                                                                       tgggcAGCCAGGC‑ACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAg        >  7:355771/5‑149 (MQ=255)                                                                                                                                       tgggcAGCCAGGC‑ACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAg        >  5:425426/5‑149 (MQ=255)                                                                                                                                        gTACAGCCTGGC‑GTGCCTTGCTGCCCAGTTCGGCGATGGCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGa       <  8:70048/149‑1 (MQ=255)                                                                                                                                        gTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGa       <  7:27173/149‑1 (MQ=255)                                                                                                                                             gCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGaa  <  5:22677/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGC‑GTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAA  >  NC_002947/3641674‑3641962

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