Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	88,649	T→A	7.1%	A223A (GCA→GCT)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	88,649	0	T	A	7.1%	293.8 / 11.4	98	A223A (GCA→GCT)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (38/53);  new base A (3/4);  total (41/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

AGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATG  >  NC_002947/88500‑88792                                                                                                                                                      |                                                                                                                                                 aGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCt                                                                                                                                                  <  2:1238888/150‑1 (MQ=255)   ggTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGc                                                                                                                                                <  1:1515316/150‑1 (MQ=255)     tCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                              <  2:1100096/150‑1 (MQ=255)      cTGCAGCAATTGGCCGATGTTGTCATCGATATGGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTg                                                                                                                                             <  2:379466/150‑1 (MQ=255)       tGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTg‑g                                                                                                                                           >  1:59425/1‑150 (MQ=255)        gcagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑Gc                                                                                                                                          >  1:634031/1‑150 (MQ=255)        gcagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑Gc                                                                                                                                          >  1:1882391/1‑150 (MQ=255)        gcagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑Gc                                                                                                                                          >  1:1882321/1‑150 (MQ=255)         cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCt                                                                                                                                         <  1:1196466/150‑1 (MQ=255)         cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCt                                                                                                                                         >  2:705599/1‑150 (MQ=255)                 ggCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGCCGAGTTCAGCCAG‑GCAGAACTCga                                                                                                                                 <  1:1739817/150‑2 (MQ=255)                  gCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGc                                                                                                                                >  1:74226/1‑150 (MQ=255)                  gCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCAGCCTG‑GCAGAACTCgac                                                                                                                                <  1:807789/150‑3 (MQ=255)                     gATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTACCAAAAAATCTCTTGTGTTGAATGCGGGGCGAGTTCTGCCTG‑GCTGAACTCGGCACg                                                                                                                             >  2:49668/1‑150 (MQ=255)                      aTGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACgg                                                                                                                            <  1:602463/150‑1 (MQ=255)                            tCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATg                                                                                                                      >  2:1241379/1‑150 (MQ=255)                              aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATggg                                                                                                                    <  2:59425/150‑1 (MQ=255)                              aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATggg                                                                                                                    <  2:1676067/150‑1 (MQ=255)                               tCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATgggg                                                                                                                   >  2:1030134/1‑150 (MQ=255)                                  atatAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  1:1612018/150‑1 (MQ=255)                                  atatAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  1:1612103/150‑1 (MQ=255)                                    ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATc                                                                                                              <  1:1134846/150‑1 (MQ=255)                                    ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCATGTGCT‑ATCTCGGCACGGGGCATGGGGATATc                                                                                                              <  1:49668/150‑1 (MQ=255)                                      aGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCca                                                                                                            >  2:343576/1‑150 (MQ=255)                                         tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                         <  2:1882391/150‑1 (MQ=255)                                         tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                         <  2:1882321/150‑1 (MQ=255)                                         tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCCGCACGGGGCATGGGGATATCCacac                                                                                                         <  1:1203114/150‑1 (MQ=255)                                            aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  2:1279914/140‑1 (MQ=255)                                            aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  2:1280051/140‑1 (MQ=255)                                            aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                >  1:1279914/1‑140 (MQ=255)                                            aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCtt                                                                                                      >  1:179669/1‑150 (MQ=255)                                            aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATACACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                >  1:1280051/1‑140 (MQ=255)                                             cgcgCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTc                                                                                                     <  1:1084774/150‑1 (MQ=255)                                                gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                  <  1:1149364/150‑1 (MQ=255)                                                gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                  >  2:1902141/1‑150 (MQ=255)                                                  cGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGc                                                                                                >  2:929503/1‑150 (MQ=255)                                                  cGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTCCGGGAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGc                                                                                                >  1:612115/1‑150 (MQ=255)                                                   gAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGACCTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCg                                                                                               <  1:1082104/150‑1 (MQ=255)                                                     aGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCGTGGGGATATCCACACCTTCGTAGCGGt                                                                                             <  2:293274/150‑1 (MQ=255)                                                       tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                           >  2:51995/1‑150 (MQ=255)                                                           cgcgGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGt                                                                                       <  2:260091/150‑1 (MQ=255)                                                            gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                      <  1:1529272/150‑1 (MQ=255)                                                                  ggCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGttt                                                                                <  2:381859/150‑1 (MQ=255)                                                                   gCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                               <  2:578607/150‑1 (MQ=255)                                                                        gcgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAAt                                                                          <  1:51995/150‑1 (MQ=255)                                                                        gcgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGGGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAAt                                                                          >  1:338585/1‑150 (MQ=255)                                                                         gatctcTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                         <  1:967596/145‑1 (MQ=255)                                                                         cgAACCTTATCCACCGGCAGTGTCTTGTGCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                         <  2:957577/150‑1 (MQ=255)                                                                           aaCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                       <  2:68845/150‑1 (MQ=255)                                                                               ttATCCACATGCAGGGGCGTGTGCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGCCATGGGGATATCCACACCTTCGTAGCGGTCGCAATAGCATTTGGCAATggagatc                                                                   >  2:967596/1‑145 (MQ=255)                                                                               ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAg                                                                   <  2:278916/150‑1 (MQ=255)                                                                                    cacaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTc                                                              >  1:452816/1‑150 (MQ=255)                                                                                    cacaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGATCGAGTTCAGCCAG‑GCAGAACTc                                                                                                                                   <  1:1600280/81‑1 (MQ=38)                                                                                    cacaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGATCGAGTTCAGCCAG‑GCAGAACTc                                                                                                                                   >  2:1600280/1‑81 (MQ=38)                                                                                     acaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCa                                                             <  1:1241379/150‑1 (MQ=255)                                                                                     acaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATACACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCa                                                             >  1:356917/1‑150 (MQ=255)                                                                                      caGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCAt                                                            <  1:445264/150‑1 (MQ=255)                                                                                      caGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGAGTCGAGTTATGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCATACCGGTCCCAGTGACGTTTGTCAATGATGTAGGTGTCAt                                                            >  2:309347/1‑150 (MQ=255)                                                                                        ggCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGACCTCGGCACGGGGCATGGGGAGATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGc                                                          <  1:705599/150‑1 (MQ=255)                                                                                          cAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgg                                                        <  2:93554/150‑1 (MQ=255)                                                                                             tGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGTCAATGGTGTAGGGGTCATGCgggtg                                                     >  1:1512/1‑150 (MQ=255)                                                                                              ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                    <  2:179669/150‑1 (MQ=255)                                                                                                cTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGCGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGGTATCCGCACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgggt                                                  >  1:1244466/1‑150 (MQ=255)                                                                                                cTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgggt                                                  >  1:692519/1‑150 (MQ=255)                                                                                                 ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTc                                                                                                     <  1:1428376/98‑1 (MQ=255)                                                                                                 ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTc                                                                                                     >  2:1428376/1‑98 (MQ=255)                                                                                                    ttCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATc                                              >  2:1257058/1‑150 (MQ=255)                                                                                                        aacggTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                          <  1:1844075/147‑1 (MQ=255)                                                                                                        aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                          >  2:414040/1‑150 (MQ=255)                                                                                                         aaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCATG‑GCTGAACTTTTACCGAAATTATATGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAAc                                         <  1:1468654/150‑1 (MQ=255)                                                                                                           ggTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACgg                                       >  1:120816/1‑150 (MQ=255)                                                                                                                tAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                               >  1:153468/1‑105 (MQ=255)                                                                                                                tAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                               <  2:153468/105‑1 (MQ=255)                                                                                                                    ccTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     >  2:1638905/1‑127 (MQ=255)                                                                                                                    ccTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     <  1:1638905/127‑1 (MQ=255)                                                                                                                     cTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGATCATGCGGGTGGGTCATCGAAACGGTCAGGCaaaa                             >  2:1338894/1‑150 (MQ=255)                                                                                                                          gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                         >  2:1840569/1‑101 (MQ=255)                                                                                                                          gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                         <  1:1840569/101‑1 (MQ=255)                                                                                                                          gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggcc                        <  1:1902141/150‑1 (MQ=255)                                                                                                                           cagGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggccg                       <  2:692519/150‑1 (MQ=255)                                                                                                                                 cTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATc                 <  1:1030134/150‑1 (MQ=255)                                                                                                                                 cTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATc                 <  1:144941/150‑1 (MQ=255)                                                                                                                                     gAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGttt             <  1:343576/150‑1 (MQ=255)                                                                                                                                       aTGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTc           <  2:1075580/150‑1 (MQ=255)                                                                                                                                       aTGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTc           <  1:1568844/150‑1 (MQ=255)                                                                                                                                       aTGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTc           <  1:1569148/150‑1 (MQ=255)                                                                                                                                         gCGGGTCGAGTTTTGCCTG‑CGTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGCGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc         <  2:502259/150‑1 (MQ=255)                                                                                                                                           gggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTGCTGCAATGGTGAAAGGGTCCAGAATGGGTGCATTTCAATCCGTCAGGCAAAATGGCCGGCCATCGTGTTCgtgc       >  1:1868174/1‑150 (MQ=255)                                                                                                                                            ggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                     <  2:1541302/119‑1 (MQ=255)                                                                                                                                            ggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                     >  1:1541302/1‑119 (MQ=255)                                                                                                                                             gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     <  1:813158/102‑1 (MQ=255)                                                                                                                                             gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     >  2:813158/1‑102 (MQ=255)                                                                                                                                             gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGcac     >  1:1440144/1‑150 (MQ=255)                                                                                                                                              ccGAGTTCAGCCAG‑GCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                    >  1:97390/2‑118 (MQ=255)                                                                                                                                              ccGAGTTCAGCCAG‑GCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                    <  2:97390/117‑1 (MQ=255)                                                                                                                                               cGAGTTCAGCCAG‑GCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACAt   >  1:1029587/1‑150 (MQ=255)                                                                                                                                                gAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     >  1:208935/1‑99 (MQ=255)                                                                                                                                                gAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                     <  2:208935/99‑1 (MQ=255)                                                                                                                                                gAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATg  <  2:1512/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATG  >  NC_002947/88500‑88792

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑