Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,149:1	+A	7.7%	coding (805/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,149	1	.	A	7.7%	286.8 / 18.9	78	coding (805/1902 nt)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base . (33/39);  new base A (3/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCC  >  NC_002947/1491001‑1491295                                                                                                                                                        |                                                                                                                                                    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:63679/1‑150 (MQ=255)                                 tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                                    >  1:661117/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                                                                                                                                                 >  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)                                                                                                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)                                                                                                                                          cacGCCACTGG‑CGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGAt            >  1:1073406/3‑150 (MQ=255)                                                                                                                                                   gg‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGGGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTcc  <  1:1629049/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCAC‑CGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGG‑CG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG‑AGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCC  >  NC_002947/1491001‑1491295

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