Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,805,118	T→A	5.6%	K925M (AAG→ATG)	PP_4234 ←	bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,805,118	0	T	A	5.6%	220.8 / 2.7	71	K925M (AAG→ATG)	PP_4234	bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

TTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804972‑4805258                                                                                                                                                   |                                                                                                                                             ttCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTg                                                                                                                                            >  4:322700/1‑149 (MQ=255) ttCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTg                                                                                                                                            >  8:38018/1‑149 (MQ=255)      gCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGAt                                                                                                                                       <  7:335650/149‑1 (MQ=255)       cAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATg                                                                                                                                      <  4:351165/149‑1 (MQ=255)       cAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGccatc                                                                                                                                       <  3:230441/148‑4 (MQ=255)          gTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCa                                                                                                                                   >  7:236526/1‑149 (MQ=255)             gTCGAGCATCGCCAGGGCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATgc                                                                                                                                <  4:16631/149‑1 (MQ=255)                gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATgcgc                                                                                                                              >  7:273770/1‑148 (MQ=255)                gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAgcgc                                                                                                                              <  2:62549/148‑1 (MQ=255)                   catcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTGAGTTGCAGTGGACGGTTAGCCATTGAGGGCCGCCTGGATATGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACg                                                                                                                          >  3:247480/1‑149 (MQ=255)                    atcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACgg                                                                                                                         <  1:194392/149‑1 (MQ=255)                    atcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACg                                                                                                                          <  8:298790/148‑1 (MQ=255)                       gCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCAt                                                                                                                      >  8:67052/1‑149 (MQ=255)                            gTCCGCGGGACCCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGcc                                                                                                                  >  1:36912/1‑148 (MQ=255)                            gTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCg                                                                                                                 <  3:273362/149‑1 (MQ=255)                                  gggAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCa                                                                                                           >  2:266846/1‑149 (MQ=255)                                   ggAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCAc                                                                                                          <  3:16276/149‑1 (MQ=255)                                        cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgac                                                                                                      >  4:299688/1‑148 (MQ=255)                                        cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTAAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgac                                                                                                      >  7:112361/1‑148 (MQ=255)                                         tCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacga                                                                                                    >  1:376008/1‑149 (MQ=255)                                           gACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCgacgac                                                                                                   <  3:168768/148‑1 (MQ=255)                                             ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGTCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCgg                                                                                                >  5:365476/1‑149 (MQ=255)                                             ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCgg                                                                                                >  7:217974/1‑149 (MQ=255)                                                    gACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGatcatc                                                                                         <  2:104065/149‑1 (MQ=255)                                                     aCCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAtcatca                                                                                        <  3:353870/149‑1 (MQ=255)                                                      ccGGTTGACCGTTGCGCGTCAGTTGCAGTGGACGTTTAGCCATTGAGGTCCTCCTTGATACGGGCTGGGGTGATGGGGATCTCGCGCACGCGCTTGCGGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATcac                                                                                       <  4:214707/149‑1 (MQ=255)                                                      ccGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGTCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATcac                                                                                       <  6:365476/149‑1 (MQ=255)                                                         gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATcaccac                                                                                    >  3:159011/1‑149 (MQ=255)                                                          ttGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                   >  3:80353/1‑149 (MQ=255)                                                            gACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTc                                                                                 >  7:281756/1‑149 (MQ=255)                                                                gTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGc                                                                              >  7:200115/1‑148 (MQ=255)                                                                       gTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCt                                                                      >  5:360856/1‑149 (MQ=255)                                                                        tCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCtt                                                                     >  7:321356/1‑149 (MQ=255)                                                                        tCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCtt                                                                     >  7:68830/1‑149 (MQ=255)                                                                          aGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTgg                                                                   <  1:292392/149‑1 (MQ=255)                                                                           gTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                               >  8:362812/1‑104 (MQ=255)                                                                           gTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                               <  7:362812/104‑1 (MQ=255)                                                                           gTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCGGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCAATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGa                                                                  <  4:240836/149‑1 (MQ=255)                                                                              gCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAggg                                                               >  6:26360/1‑149 (MQ=255)                                                                              gCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAgg                                                                <  3:322700/148‑1 (MQ=255)                                                                               cAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAggg                                                               <  8:273770/148‑1 (MQ=255)                                                                                aGTGGTCGGTTAGCACTTTAGGGCCTCCTTGAGACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGggt                                                             <  6:267944/149‑1 (MQ=255)                                                                                aGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGggt                                                             <  8:112361/149‑1 (MQ=255)                                                                                    gACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGggtcggt                                                         <  7:67052/149‑1 (MQ=255)                                                                                        gTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                     <  4:115030/149‑1 (MQ=255)                                                                                        gTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                     <  4:80353/149‑1 (MQ=255)                                                                                         ttAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                     >  8:78761/1‑148 (MQ=255)                                                                                             ccATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAg                                                >  7:63674/1‑149 (MQ=255)                                                                                                   aGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATg                                                                          <  1:242362/117‑1 (MQ=255)                                                                                                   aGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATg                                                                          >  2:242362/1‑117 (MQ=255)                                                                                                     ggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGcggcgg                                          >  1:91049/1‑147 (MQ=255)                                                                                                       cctcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGggagg                                      <  2:36912/149‑1 (MQ=255)                                                                                                        ctcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGggagga                                     >  6:372545/1‑149 (MQ=255)                                                                                                         tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAt                                    >  6:17300/1‑149 (MQ=255)                                                                                                         tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAt                                    >  4:139576/1‑149 (MQ=255)                                                                                                              gATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGc                               >  4:357623/1‑149 (MQ=255)                                                                                                                   gggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTc                          <  6:360856/149‑1 (MQ=255)                                                                                                                      cTGGGGTGAGGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgc                       <  2:205940/149‑1 (MQ=255)                                                                                                                       tGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgct                      <  8:236526/149‑1 (MQ=255)                                                                                                                        ggggTGTTGGGGAGATCGCGCACGCGCTTGCCGTTCTGATGCGCGACTGGATTGGCGATTGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGTGGGGTCGGGCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                     <  4:247480/149‑1 (MQ=255)                                                                                                                        ggggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                     >  4:214281/1‑149 (MQ=255)                                                                                                                              aTGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCcaca                <  3:299688/148‑1 (MQ=255)                                                                                                                                      aTCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACAt       >  4:111035/1‑149 (MQ=255)                                                                                                                                       tCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATc      <  3:357623/149‑1 (MQ=255)                                                                                                                                          cgcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcg   <  3:214281/149‑1 (MQ=255)                                                                                                                                          cgcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcg   <  1:266846/149‑1 (MQ=255)                                                                                                                                          ccgtcgcgCATGCCATCGGCAAGCGCGTCGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcg   >  2:120821/5‑149 (MQ=255)                                                                                                                                           gcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGc                                                                         >  2:48801/1‑78 (MQ=255)                                                                                                                                           gcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGc                                                                         <  1:48801/78‑1 (MQ=255)                                                                                                                                           gcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcgc  >  3:310749/1‑149 (MQ=255)                                                                                                                                           cgtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCCCCACTTCGGCGATGCCCTTGGCGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcgc  >  3:253323/4‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                             TTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804972‑4805258

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