Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,491,957	A→C	5.7%	N154K (AAT→AAG)	PP_3096 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,491,957	0	A	C	5.7%	220.0 / 2.4	69	N154K (AAT→AAG)	PP_3096	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base C (2/2);  total (35/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 = 9.98e-01

GCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGG  >  NC_002947/3491813‑3492096                                                                                                                                                 |                                                                                                                                            gcgcGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTg                                                                                                                                          <  1:199364/148‑1 (MQ=255)   gcgGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTTAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCCTTgcgc                                                                                                                                       >  7:66353/1‑149 (MQ=255)   gcgGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCCTTgcgc                                                                                                                                       >  7:40017/1‑149 (MQ=255)   gcgGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTgcgc                                                                                                                                       >  3:137609/1‑149 (MQ=255)      gCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCg                                                                                                                                    >  6:328278/1‑149 (MQ=255)        ttCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGACGCCCGAGGTGCCATTGCGCTCGTg                                                                                                                                  >  6:251105/1‑149 (MQ=255)                gAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAACGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGAtgt                                                                                                                          >  4:126957/1‑149 (MQ=255)                      gCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAAtt                                                                                                                    <  4:309473/149‑1 (MQ=255)                         ggCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCAtt                                                                                                                                           >  4:30200/1‑123 (MQ=255)                           cAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTAt                                                                                                               <  1:65495/149‑1 (MQ=255)                           cAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTAt                                                                                                               <  8:316897/149‑1 (MQ=255)                              aaTAGCTGAACAGCTCACGACGCAACGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTAtttt                                                                                                            <  3:126957/149‑1 (MQ=255)                                   cTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGc                                                                                                         >  4:81306/1‑147 (MQ=255)                                    tGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAg                                                                                                      >  8:175893/1‑149 (MQ=255)                                      aaCAGCTTACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGc                                                                                                    >  2:181839/1‑149 (MQ=255)                                      aaCAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTg                                                                                                                                          <  3:150462/111‑1 (MQ=255)                                      aaCAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTg                                                                                                                                          >  4:150462/1‑111 (MQ=255)                                      aaCAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGc                                                                                                    >  6:105904/1‑149 (MQ=255)                                           cTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTc                                                                                                                                     >  5:334391/1‑111 (MQ=255)                                           cTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTc                                                                                                                                     <  6:334391/111‑1 (MQ=255)                                           cTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCg                                                                                                >  6:158919/1‑148 (MQ=255)                                            tCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCg                                                                                                >  5:208759/1‑147 (MQ=255)                                             cACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCg                                                                                                >  2:202255/1‑146 (MQ=255)                                              acgacgCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAg                                                                                             <  2:95075/148‑1 (MQ=255)                                                  cgCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGaacaa                                                                                        >  4:45702/1‑149 (MQ=255)                                                   gCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTg                                                                                                                  >  5:328813/1‑122 (MQ=255)                                                   gCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTg                                                                                                                  <  6:328813/122‑1 (MQ=255)                                                   gCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGaacaac                                                                                       <  7:36038/149‑1 (MQ=255)                                                       tGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTgg                                                                                   <  3:81306/149‑1 (MQ=255)                                                                 tGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTg                                                                         >  7:93472/1‑149 (MQ=255)                                                                   ctctTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGc                                                                        <  4:137609/148‑1 (MQ=255)                                                                     ctTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTAt                                                                                                               <  3:189881/107‑1 (MQ=255)                                                                     ctTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTAt                                                                                                               >  4:189881/1‑107 (MQ=255)                                                                      ttCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCgag                                                                     >  7:134500/1‑148 (MQ=255)                                                                      ttCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCga                                                                      <  5:328278/147‑1 (MQ=255)                                                                        cGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCgaga                                                                    <  8:93472/147‑1 (MQ=255)                                                                        cGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCgag                                                                     <  5:105904/146‑1 (MQ=255)                                                                        cGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGagaa                                                                   <  3:45702/148‑1 (MQ=255)                                                                        cGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGCAGAACAACGTGGTCATGCGGTGCgag                                                                     <  1:181839/146‑1 (MQ=255)                                                                           aCTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTg                                                                                    <  7:161469/128‑1 (MQ=255)                                                                           aCTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTg                                                                                    >  8:161469/1‑128 (MQ=255)                                                                           aCTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGagaa                                                                   <  4:161843/145‑1 (MQ=255)                                                                           aCTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGagaa                                                                   >  3:161843/1‑145 (MQ=255)                                                                            cTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCCTTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGt                                                              <  8:40017/149‑1 (MQ=255)                                                                            cTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCCTTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGt                                                              <  8:66353/149‑1 (MQ=255)                                                                             ttCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTc                                                             <  7:175893/149‑1 (MQ=255)                                                                             ttCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTc                                                             >  6:24732/1‑149 (MQ=255)                                                                               cAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGc                                                                        <  1:174695/136‑1 (MQ=255)                                                                               cAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGc                                                                        >  2:174695/1‑136 (MQ=255)                                                                                aGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCga                                                           <  5:24732/148‑1 (MQ=255)                                                                                            cgcgGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCTGTGCGAGAAGATGTCGAGAAACGCATGggc                                              <  5:158919/149‑1 (MQ=255)                                                                                                 cccGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATg                                                                              <  5:337694/112‑1 (MQ=255)                                                                                                 cccGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATg                                                                              >  6:337694/1‑112 (MQ=255)                                                                                                       gCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTcgc                                   <  8:134500/149‑1 (MQ=255)                                                                                                        ccAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCga                                                                      >  4:75753/1‑113 (MQ=255)                                                                                                        ccAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCga                                                                      <  3:75753/113‑1 (MQ=255)                                                                                                         cAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCga                                                           <  4:37995/123‑1 (MQ=255)                                                                                                         cAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCga                                                           >  3:37995/1‑123 (MQ=255)                                                                                                             cAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACt                             >  1:171289/1‑149 (MQ=255)                                                                                                                gTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGagaag                                                                  >  2:255572/1‑109 (MQ=255)                                                                                                                gTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGagaag                                                                  <  1:255572/109‑1 (MQ=255)                                                                                                                              cATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagc            <  6:208759/149‑1 (MQ=255)                                                                                                                                tCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGAATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcag          >  7:52920/1‑149 (MQ=255)                                                                                                                                  gAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCgag                                                          >  2:193447/1‑99 (MQ=255)                                                                                                                                  gAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCgag                                                          <  1:193447/99‑1 (MQ=255)                                                                                                                                  gAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcag          >  2:87822/1‑147 (MQ=255)                                                                                                                                        cGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCgg  >  5:73735/1‑149 (MQ=255)                                                                                                                                               ccATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcagc         >  3:324689/1‑135 (MQ=255)                                                                                                                                               ccATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcagc         <  4:324689/135‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            GCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGTTGCTCTTCGAACTTCAGCTTCTGGGCGCGCGGCCCGAGGCCAACCAGGTTGAGCAGGTAGCCATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGG  >  NC_002947/3491813‑3492096

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