Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	88,657	T→A	7.5%	S221C (AGC→TGC)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	88,657	0	T	A	7.5%	274.7 / 11.0	93	S221C (AGC→TGC)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (34/52);  new base A (3/4);  total (37/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 = 7.10e-01

CAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACT  >  NC_002947/88508‑88806                                                                                                                                                       |                                                                                                                                                      cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCt                                                                                                                                                       >  2:705599/1‑150 (MQ=255) cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCt                                                                                                                                                       <  1:1196466/150‑1 (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: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                                                                                                                              >  1:1279914/1‑140 (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                                                                                                                                                 >  2:1600280/1‑81 (MQ=38)                                                                            cacaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGATCGAGTTCAGCCAG‑GCAGAACTc                                                                                                                                                 <  1:1600280/81‑1 (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                                                                                       <  1:1840569/101‑1 (MQ=255)                                                                                                                  gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                       >  2:1840569/1‑101 (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:1569148/150‑1 (MQ=255)                                                                                                                               aTGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTc                         <  1:1568844/150‑1 (MQ=255)                                                                                                                                 gCGGGTCGAGTTTTGCCTG‑CGTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGCGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc                       <  2:502259/150‑1 (MQ=255)                                                                                                                                   gggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTGCTGCAATGGTGAAAGGGTCCAGAATGGGTGCATTTCAATCCGTCAGGCAAAATGGCCGGCCATCGTGTTCgtgc                     >  1:1868174/1‑150 (MQ=255)                                                                                                                                    ggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                                   >  1:1541302/1‑119 (MQ=255)                                                                                                                                    ggTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                                   <  2:1541302/119‑1 (MQ=255)                                                                                                                                     gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                                   >  2:813158/1‑102 (MQ=255)                                                                                                                                     gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                                   <  1:813158/102‑1 (MQ=255)                                                                                                                                     gTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGcac                   >  1:1440144/1‑150 (MQ=255)                                                                                                                                      ccGAGTTCAGCCAG‑GCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                                  <  2:97390/117‑1 (MQ=255)                                                                                                                                      ccGAGTTCAGCCAG‑GCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                                  >  1:97390/2‑118 (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)                                                                                                                                                   g‑gCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGt     <  2:120816/150‑1 (MQ=255)                                                                                                                                                     gCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTTATCGAACCGATCCGGCAACATGGCCGGCCATCGTTTTCGGGCACATGGTCGTACAGGTa    <  1:1476780/150‑1 (MQ=255)                                                                                                                                                       tGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAGTGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACt  >  2:481679/1‑150 (MQ=255)                                                                                                                                                       tGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACt  >  2:1076324/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTG‑GCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACT  >  NC_002947/88508‑88806

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