Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,144	G→C	8.8%	S267T (AGT→ACT)	cysNC →	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,491,144	0	G	C	8.8%	270.4 / 16.3	91	S267T (AGT→ACT)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (38/45);  new base C (4/4);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAT  >  NC_002947/1490998‑1491285                                                                                                                                                    |                                                                                                                                                 gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                               <  2:889068/150‑1 (MQ=255)  ccGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑c                                                                                                                                             >  2:1420504/1‑150 (MQ=255)    gTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑t                                                                                                                                          <  1:78273/150‑1 (MQ=255)     tGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tg                                                                                                                                         <  2:1737325/150‑1 (MQ=255)     cggTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCTTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tg                                                                                                                                         <  1:398263/149‑1 (MQ=255)        tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggtg                                                                                                                                      >  2:1388826/1‑150 (MQ=255)        tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑CGAtggtg                                                                                                                                      <  2:1237848/150‑1 (MQ=255)         aCGCCGGCCCGACGCTGATGGAAATCCTCGAACCCGTCGAGGTTTAGCCTGCCCGCAACCTCGCGCGA‑CAACTTCTCCCTGGGCTGTACCTCGACCTCGCTGTCCGGAATTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑CGATGGTGc                                                                                                                                     <  2:1161581/150‑1 (MQ=255)               gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAg                                                                                                                               >  1:749822/1‑150 (MQ=255)                 ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGt                                                                                                                             <  2:78691/150‑1 (MQ=255)                   tactcTGATTGAAATCCTCGCAACCGTCGAGGTTTCGGTTGCCCGCAACTTCAC‑CTTCCTGCTTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGa                                                                                                                           <  2:1260733/146‑1 (MQ=255)                      gCTGATGGAACTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACga                                                                                                                        <  1:442618/150‑1 (MQ=255)                      gCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAGCATCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACga                                                                                                                        >  2:1485528/1‑150 (MQ=255)                      gCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACga                                                                                                                        <  2:1169049/150‑1 (MQ=255)                       cTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACgag                                                                                                                       <  1:1420817/150‑1 (MQ=255)                            ggAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGt                                                                                                                  >  2:572309/1‑150 (MQ=255)                             gAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                 <  2:222416/150‑1 (MQ=255)                               aaTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGt                                                                                                               <  1:1896944/150‑1 (MQ=255)                               aaTCCTCGAAACCGTAGATGGTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGt                                                                                                               >  1:1894876/1‑150 (MQ=255)                                    tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                          >  1:661117/1‑150 (MQ=255)                                    tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                          >  1:63679/1‑150 (MQ=255)                                    tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTGCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTTAACTTCCGCGGCTTCGCGGGCACCAGCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                          <  1:739589/150‑1 (MQ=255)                                     cGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGcc                                                                                                         <  2:766068/150‑1 (MQ=255)                                      gAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCg                                                                                                        >  1:1036363/1‑150 (MQ=255)                                          ccGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCACCCTGAACGTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                    <  1:212708/150‑1 (MQ=255)                                          ccGTCGAGATTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGCACGTCAACCGTCCCACCCTGAACTCCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                    <  2:1039264/150‑1 (MQ=255)                                             tCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCa                                                                                                 >  2:1786296/1‑150 (MQ=255)                                             tCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCACCGCCAGTGG‑CG‑TGGTGca                                                                                                                                    <  1:217969/115‑1 (MQ=255)                                             tCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACATGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCACCGCCAGTGG‑CG‑TGGTGca                                                                                                                                    >  2:217969/1‑115 (MQ=255)                                                aGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAag‑a                                                                                             <  1:550639/150‑1 (MQ=255)                                                     tCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑CGAtggtg                                                                                                                                      >  1:427179/1‑105 (MQ=255)                                                     tCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑CGAtggtg                                                                                                                                      <  2:427179/105‑1 (MQ=255)                                                        gCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCgg                                                                                     <  1:127831/150‑1 (MQ=255)                                                          tGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       >  2:1647685/1‑99 (MQ=255)                                                          tGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       <  1:1647685/99‑1 (MQ=255)                                                              cGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCAGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAg                                                                               >  1:332835/1‑150 (MQ=255)                                                               gCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGt                                                                              <  2:63679/150‑1 (MQ=255)                                                                  aCTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCTCGGCTTCGCGGGCACCATCGCCCGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCa                                                                           <  1:1075233/150‑1 (MQ=255)                                                                   cTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTCCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCAt                                                                          >  2:579831/1‑150 (MQ=255)                                                                   cTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑Cga                                                                                                                                           <  2:146518/86‑2 (MQ=255)                                                                   cTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑Cga                                                                                                                                           >  1:146518/1‑85 (MQ=255)                                                                    ttCAC‑CGACCTGCTGTTCCCGGTGCAGTACAGCGACCGTCCCAACCTGAATTTCCGCGGGTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGATGAGATTGTCGTGCTGCCGGCGGGCAAG‑AGCAGCCGGGTCAAGTCCATc                                                                         <  2:850733/150‑1 (MQ=255)                                                                    ttCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCCGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATc                                                                         >  2:669059/1‑150 (MQ=255)                                                                    ttCAC‑CGACCTGCGTTTCCCGGCGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGTTGCCGTCGCGCAAG‑AGCAGCCGGGTCAAGTCCATc                                                                         >  2:1338638/1‑150 (MQ=255)                                                                     tCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCg                                                                        >  1:590404/1‑150 (MQ=255)                                                                     tCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCg                                                                        <  1:1420504/150‑1 (MQ=255)                                                                       aC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTc                                                                      >  1:121577/1‑150 (MQ=255)                                                                       aC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTc                                                                      <  1:341194/150‑1 (MQ=255)                                                                          cGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGTCG‑CGG‑GCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCAc                                                                    >  1:329175/1‑150 (MQ=255)                                                                          cGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCAc                                                                    >  2:248964/1‑150 (MQ=255)                                                                           gACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTcccc                                                                   <  1:1388826/150‑1 (MQ=255)                                                                                    ttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGTCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGt                                                          <  1:188793/150‑1 (MQ=255)                                                                                        cGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGCGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                      <  1:817171/150‑1 (MQ=255)                                                                                          gTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTg                                                    <  2:1598057/150‑1 (MQ=255)                                                                                           tGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTgg                                                   <  2:121577/150‑1 (MQ=255)                                                                                            gCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCGAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGa                                                  <  2:1036363/150‑1 (MQ=255)                                                                                                   gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGt                                                                                                                                                 >  2:285486/1‑51 (MQ=255)                                                                                                   gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGt                                                                                                                                                 <  1:285486/51‑1 (MQ=255)                                                                                                   gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCt                                           >  2:1526064/1‑150 (MQ=255)                                                                                                   gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGAAGACTTTGTCGTGCTGCCGTCGGGAAAGCACCGGCAGGGACAGGT‑GTTCGTGACCTACGAAGGTGAACTGGAAAACGCt                                           >  1:888317/1‑150 (MQ=255)                                                                                                    tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       <  1:1753696/58‑1 (MQ=255)                                                                                                    tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       >  2:1753696/1‑58 (MQ=255)                                                                                                    tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                          <  2:139454/150‑1 (MQ=255)                                                                                                     cAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                         <  1:1897897/150‑1 (MQ=255)                                                                                                      aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGc                                        >  2:1786371/1‑150 (MQ=255)                                                                                                      aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGc                                        <  1:849807/150‑1 (MQ=255)                                                                                                        ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACt                                                     <  1:517918/135‑1 (MQ=255)                                                                                                        ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACt                                                     >  2:517918/1‑135 (MQ=255)                                                                                                        ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaaa                                               >  1:1335908/1‑141 (MQ=255)                                                                                                        ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaaa                                               <  2:1335908/141‑1 (MQ=255)                                                                                                        ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAGGCTGGccc                                      >  2:1042284/1‑150 (MQ=255)                                                                                                           tCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccagg                                   >  1:225779/1‑150 (MQ=255)                                                                                                              cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggcca                                >  1:599503/1‑150 (MQ=255)                                                                                                                 ccTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCctgggccagcg                             <  2:1163426/150‑3 (MQ=255)                                                                                                                 ccTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCcaggccaggc                             <  1:1786296/150‑1 (MQ=255)                                                                                                                  cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCaggccaggcc                            <  2:329175/150‑1 (MQ=255)                                                                                                                  cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGCGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCaggccaggcc                            >  1:1604596/1‑150 (MQ=255)                                                                                                                  cTGAACTTCCGCGGCTTCGCGGGCACCACCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCaggccaggcc                            <  2:661117/150‑1 (MQ=255)                                                                                                                   tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCg                           <  2:767785/150‑1 (MQ=255)                                                                                                                    gAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                          >  2:1595026/1‑150 (MQ=255)                                                                                                                      aCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCATGCCAGGCCGTGa                        <  1:1073578/150‑1 (MQ=255)                                                                                                                        ttCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                  >  1:1081028/1‑138 (MQ=255)                                                                                                                        ttCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                  <  2:1081028/138‑1 (MQ=255)                                                                                                                        ttCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGAcc                      >  2:227636/1‑150 (MQ=255)                                                                                                                                cTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCAGAAAGGTGACGAGATTGTCGTGGTGCAGTCGGGCAAG‑AGCAGCAGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGACCCAGGCCAGGCCGTGACCCTGACCAt              >  2:927225/1‑150 (MQ=255)                                                                                                                                 ttCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACAGAGGGGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATg             >  1:115186/1‑150 (MQ=255)                                                                                                                                        gCACCA‑CGCCACTGG‑CGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACg      >  1:1581239/1‑150 (MQ=255)                                                                                                                                            cATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAt  >  1:850580/1‑150 (MQ=255)                                                                                                                                            cATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAt  <  1:1526064/150‑1 (MQ=255)                                                                                                                                            aaTCGCGAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTGGTGCTGGCGTAGGGCAAG‑AGCAGCCGTGTTTAGTCGCTAGTCACATACGAAGGTGAACTGGACAACGCTGGCCCAGGCCAGGCCGTGACCCTGAGCATGGAAGAGGCGAt  >  2:826185/2‑150 (MQ=255)                                                                                                                                             cacGCCACTGG‑CGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAt  >  1:1073406/3‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAT  >  NC_002947/1490998‑1491285

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