Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,138	Δ1 bp	7.6%	coding (794/1902 nt)	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,138	0	T	.	7.6%	246.9 / 17.2	92	coding (794/1902 nt)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base . (3/4);  total (43/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.72e-01

AATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAT  >  NC_002947/1490994‑1491285                                                                                                                                                  |                                                                                                                                                       aaTCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCa                                                                                                                                                   >  2:1420817/1‑150 (MQ=255) aaTCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCa                                                                                                                                                   >  1:139454/1‑150 (MQ=255)     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                                                                                                                                      <  2:427179/105‑1 (MQ=255)                                                         tCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGG‑CGAtggtg                                                                                                                                      >  1:427179/1‑105 (MQ=255)                                                            gCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCgg                                                                                     <  1:127831/150‑1 (MQ=255)                                                              tGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       <  1:1647685/99‑1 (MQ=255)                                                              tGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       >  2:1647685/1‑99 (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:341194/150‑1 (MQ=255)                                                                           aC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTc                                                                      >  1:121577/1‑150 (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                                                                                                                                                 <  1:285486/51‑1 (MQ=255)                                                                                                       gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGt                                                                                                                                                 >  2:285486/1‑51 (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                                                                                                                                       >  2:1753696/1‑58 (MQ=255)                                                                                                        tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑Cg‑tggt                                                                                                                                       <  1:1753696/58‑1 (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                                  <  2:1081028/138‑1 (MQ=255)                                                                                                                            ttCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                  >  1:1081028/1‑138 (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)                                                                                                                                                  |                                                                                                                                                       AATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAT  >  NC_002947/1490994‑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: ‑