Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	88,653	T→A	6.1%	Q222L (CAG→CTG)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	88,653	0	T	A	6.1%	292.1 / 9.0	98	Q222L (CAG→CTG)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (39/53);  new base A (3/3);  total (42/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

TCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATG  >  NC_002947/88504‑88792                                                                                                                                                      |                                                                                                                                             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:1676067/150‑1 (MQ=255)                          aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATggg                                                                                                                    <  2:59425/150‑1 (MQ=255)                           tCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATgggg                                                                                                                   >  2:1030134/1‑150 (MQ=255)                              atatAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  1:1612103/150‑1 (MQ=255)                              atatAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  1:1612018/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:1882321/150‑1 (MQ=255)                                     tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                         <  2:1882391/150‑1 (MQ=255)                                     tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCCGCACGGGGCATGGGGATATCCacac                                                                                                         <  1:1203114/150‑1 (MQ=255)                                        aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                >  1:1279914/1‑140 (MQ=255)                                        aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  2:1279914/140‑1 (MQ=255)                                        aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGata                                                                                                                <  2:1280051/140‑1 (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                                                                                                  >  2:1902141/1‑150 (MQ=255)                                            gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                  <  1:1149364/150‑1 (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                                                                                                     >  2:1428376/1‑98 (MQ=255)                                                                                             ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTc                                                                                                     <  1:1428376/98‑1 (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)                                                                                                                                                     cTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCgg                                                                                              <  2:1064018/49‑1 (MQ=255)                                                                                                                                                     cTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCgg                                                                                              >  1:1064018/1‑49 (MQ=255)                                                                                                                                                      |                                                                                                                                             TCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATG  >  NC_002947/88504‑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: ‑