Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	88,807	G→T	5.1%	Q171K (CAG→AAG)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	88,807	0	G	T	5.1%	248.1 / 2.4	80	Q171K (CAG→AAG)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base T (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

AACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGT  >  NC_002947/88659‑88956                                                                                                                                                     |                                                                                                                                                      aaCTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTgg                                                                                                                                                      <  2:74226/150‑1 (MQ=255)    tCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCTAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgcgc                                                                                                                                                   >  2:785114/1‑150 (MQ=255)    tCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgcgc                                                                                                                                                   >  2:1641969/1‑150 (MQ=255)    tCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgcgc                                                                                                                                                   >  1:56637/1‑150 (MQ=255)               cATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAacac                                                                                                                                        <  2:841110/150‑1 (MQ=255)                  ggggATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAAcaccac                                                                                                                                     >  2:86394/1‑150 (MQ=255)                        atCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACctcctc                                                                                                                               <  1:1804911/150‑1 (MQ=255)                           cacaCCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgcgc                                                                                                                                                   <  1:1044888/127‑1 (MQ=255)                           cacaCCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgcgc                                                                                                                                                   >  2:1044888/1‑127 (MQ=255)                            acacCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCtcgtcg                                                                                                                           >  1:1052938/1‑150 (MQ=255)                                cTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAAt                                                                                                                       >  2:1004189/1‑150 (MQ=255)                                cTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGAGCGCCTTGAACACCACCTCCTCGTCGAAAt                                                                                                                       >  2:774967/1‑150 (MQ=255)                                 ttCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATc                                                                                                                      >  1:396778/1‑150 (MQ=255)                                 tacGTAGCGGACCCAGCCGCGAGTGGCAATGAAGTAGGGCACAAGCGGGAGGGTCATAGATACGGTCAGGCAAAATAGCCGGCCATCGTTAACGCGCACATGGTCGTACAGCTACTCGCGAGCATTGACCACCACTTCCGCGTCGaacac                                                                                                                      >  2:802667/3‑147 (MQ=255)                                   cGTAGCGGTCCCAGTAGCGTTTGGCAAGGGTGTAGGGGTCATGCTGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCa                                                                                                                    <  2:55575/150‑1 (MQ=255)                                     tAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAAt                                                                                                                       >  2:1448279/1‑145 (MQ=255)                                     tAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAAt                                                                                                                       <  1:1448279/145‑1 (MQ=255)                                     tAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGc                                                                                                                  <  1:414040/150‑1 (MQ=255)                                       gCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGCTTATGCTGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTg                                                                                                                >  1:1893084/1‑150 (MQ=255)                                        cGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGa                                                                                                               <  2:452816/150‑1 (MQ=255)                                                       ttGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCgcacgcac                                                                                                <  1:1641969/150‑1 (MQ=255)                                                       ttGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCgcacgcac                                                                                                <  1:1004189/150‑1 (MQ=255)                                                       ttGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCgcacgcac                                                                                                <  1:1338894/150‑1 (MQ=255)                                                          gCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGa                                                                                             <  1:929503/150‑1 (MQ=255)                                                            aaTGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTTGAAATCCAGCTGATTGGTGCGCACGCACGGAcc                                                                                           >  1:318511/1‑150 (MQ=255)                                                            aaTGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGAcc                                                                                           >  2:64388/1‑150 (MQ=255)                                                            aaTGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGAcc                                                                                           >  1:499868/1‑150 (MQ=255)                                                               cgtgcAGAGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACGTGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGc                                                                                        <  2:612115/149‑1 (MQ=255)                                                                gtgtAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                                                       >  2:1292932/1‑150 (MQ=255)                                                                   tAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgc                                                                                    >  1:233347/1‑150 (MQ=255)                                                                    tgcggTCATGCGCGGGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgca                                                                                   <  2:1047948/147‑1 (MQ=255)                                                                     ggggTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgcag                                                                                  >  2:1574342/1‑150 (MQ=255)                                                                      gggTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgcagc                                                                                 >  2:1257073/1‑150 (MQ=255)                                                                         tCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCAcc                                                                              <  1:1257058/150‑1 (MQ=255)                                                                          cATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCT‑CACGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGa                                                                            <  2:338585/150‑1 (MQ=255)                                                                           aTGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGa                                                                            <  1:86394/150‑1 (MQ=255)                                                                                     gtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGaa                                                                                                                                           >  2:432098/1‑77 (MQ=255)                                                                                     gtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGa                                                                            >  1:922974/1‑140 (MQ=255)                                                                                     gtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGa                                                                            <  2:922974/140‑1 (MQ=255)                                                                                     ctCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGaa                                                                                                                                           <  1:432098/76‑1 (MQ=255)                                                                                      tCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACAtgttgt                                                                 >  1:1724466/1‑150 (MQ=255)                                                                                      tCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACAtgttgt                                                                 >  2:1211338/1‑150 (MQ=255)                                                                                       cATCGAAACGGTCAGGCAAAATGGCCGGCGCTTCTTTCCTCTTACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                                                <  2:1893084/150‑1 (MQ=255)                                                                                       cATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                                                <  1:1574342/150‑1 (MQ=255)                                                                                        aTCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTgg                                                               <  2:634031/150‑1 (MQ=255)                                                                                        aTCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTgg                                                               >  2:461122/1‑150 (MQ=255)                                                                                              aCGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAg                                                         <  2:1029587/150‑1 (MQ=255)                                                                                                 gTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCtt                                                      <  2:396778/150‑1 (MQ=255)                                                                                                  tCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTAc                                                            >  2:57480/1‑143 (MQ=255)                                                                                                  tCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTAc                                                            <  1:57480/143‑1 (MQ=255)                                                                                                    aGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGa                                                   >  2:1830856/1‑150 (MQ=255)                                                                                                    aGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCCCCCCGTCGAAATCCGGGTGAGTGGTGCGCTCGGACGGGCGCGCCGTCGGCACCGAGGACAGGCTGTGGTACCAGCCTGGa                                                   >  2:1653698/1‑150 (MQ=255)                                                                                                    aGGCAAAAGGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGa                                                   >  2:1268826/1‑150 (MQ=255)                                                                                                     ggCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTgc                                                                                                        >  2:941129/1‑96 (MQ=255)                                                                                                     ggCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTgc                                                                                                        <  1:941129/96‑1 (MQ=255)                                                                                                      gCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACg                                                 >  2:289009/1‑150 (MQ=255)                                                                                                              gccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg                                         <  1:615850/150‑1 (MQ=255)                                                                                                              gccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg                                         <  2:233347/150‑1 (MQ=255)                                                                                                              gccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg                                         <  1:1268826/150‑1 (MQ=255)                                                                                                                  gccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCGGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTc                                     >  2:163706/1‑150 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTTGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTc                                  >  1:1758376/1‑150 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTTGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTc                                  >  1:1758292/1‑150 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgca                                                                                   >  2:717082/1‑101 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgca                                                                                   <  1:717082/101‑1 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAg                                                         >  2:1194607/1‑127 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAg                                                         <  1:1194607/127‑1 (MQ=255)                                                                                                                     aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTc                                  >  1:1157934/1‑150 (MQ=255)                                                                                                                      tCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTcc                                 <  2:1724466/150‑1 (MQ=255)                                                                                                                        gTTTTCGCGCACATGGTCGTACAGGTACTTGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCa                               <  2:1758292/150‑1 (MQ=255)                                                                                                                        gTTTTCGCGCACATGGTCGTACAGGTACTTGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCa                               <  2:1758376/150‑1 (MQ=255)                                                                                                                             cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                                                <  2:119767/112‑1 (MQ=255)                                                                                                                             cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                                                >  1:119767/1‑112 (MQ=255)                                                                                                                             cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCa                          <  1:1211338/150‑1 (MQ=255)                                                                                                                                gcACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGc                                                >  1:893496/1‑125 (MQ=255)                                                                                                                                gcACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGc                                                <  2:893496/125‑1 (MQ=255)                                                                                                                                gcACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCg                       >  1:218338/1‑150 (MQ=255)                                                                                                                                    aTGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCGGCGGATTGGTGCGCTCGCACGGACCCGCCTGCAGCACCGAGGACATGTTGGGGTACCAGCTTGGACGCAAATCGGGTTCGTCCCAGTTCACAGCCCa                   >  2:865956/1‑150 (MQ=255)                                                                                                                                          gTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTa             >  1:158010/1‑150 (MQ=255)                                                                                                                                              aGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCg         >  2:887050/1‑150 (MQ=255)                                                                                                                                                 tACTGGCGCGCCTTGAGCACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCAGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTGGTCGGcc      >  2:977906/1‑150 (MQ=255)                                                                                                                                                 tACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGcc      <  1:481679/150‑1 (MQ=255)                                                                                                                                                     ggCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGt  <  1:461122/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGT  >  NC_002947/88659‑88956

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