Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,581	T→A	8.2%	intergenic (+39/‑50)	iscX → / → ndk	cysteine desulfurase activity regulator/nucleoside diphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	983,581	0	T	A	8.2%	177.5 / 5.9	61	intergenic (+39/‑50)	iscX/ndk	cysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (3/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAA  >  NC_002947/983434‑983729                                                                                                                                                    |                                                                                                                                                        ccTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGta                                                                                                                                                        >  8:299792/1‑149 (MQ=255)             ggTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCgg‑g                                                                                                                                            <  6:199176/148‑1 (MQ=255)                       ttGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAAttttt                                                                                                                                  <  4:940/148‑1 (MQ=255)                            agagTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTaa                                                                                                                             >  2:199946/1‑148 (MQ=255)                                ttCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACt                                                                                                                         >  1:65236/1‑148 (MQ=255)                                             cGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCATGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGtt                                                                                                           >  6:154293/1‑149 (MQ=255)                                              gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGt                                                                                                                               >  4:15111/1‑128 (MQ=255)                                              gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGt                                                                                                                               <  3:15111/128‑1 (MQ=255)                                              gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGttt                                                                                                          >  5:42295/1‑149 (MQ=255)                                              gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGttt                                                                                                          >  5:36984/1‑149 (MQ=255)                                                   gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑t                                                                                                   <  8:111572/149‑1 (MQ=255)                                                         gTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGt                                                                                             <  6:36984/149‑1 (MQ=255)                                                         gTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGt                                                                                             <  6:42295/149‑1 (MQ=255)                                                            agaAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTc                                                                                           >  7:256095/1‑148 (MQ=255)                                                             gaAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                              <  6:131394/130‑1 (MQ=255)                                                             gaAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                              >  5:131394/1‑130 (MQ=255)                                                              aaaGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCaa                                                                                         <  7:299792/148‑1 (MQ=255)                                                                     tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTGCGTAACAGGCTAATTATGGAGTTTACCT‑‑TGGGTGTTAAACGTCCt                                                                                   >  4:250014/1‑147 (MQ=255)                                                                        aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCATGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTtctc                                                                              <  5:154293/149‑1 (MQ=255)                                                                        aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGGCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTtctc                                                                              <  5:350620/149‑1 (MQ=255)                                                                         ggCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTcc                                                                             <  8:33444/149‑1 (MQ=255)                                                                            aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCcatc                                                                          <  8:256095/149‑1 (MQ=255)                                                                                   gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGc                                                                   <  2:65236/149‑1 (MQ=255)                                                                                        cTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGAc                                                              >  4:295898/1‑149 (MQ=255)                                                                                             tCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTg                                                                                                                >  2:136265/1‑96 (MQ=255)                                                                                             tCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTg                                                                                                                <  1:136265/96‑1 (MQ=255)                                                                                               gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTg                                                       >  6:149003/1‑149 (MQ=255)                                                                                               gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTg                                                       <  1:78001/149‑1 (MQ=255)                                                                                                    aGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTg                                                          <  8:23895/141‑1 (MQ=255)                                                                                                    aGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTg                                                          >  7:23895/1‑141 (MQ=255)                                                                                                     gCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCaaa                                                  >  2:232385/1‑148 (MQ=255)                                                                                                           tAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCa                                           >  7:140009/1‑149 (MQ=255)                                                                                                            aaGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCAt                                          >  8:242741/1‑149 (MQ=255)                                                                                                             aGCGTCACTTGCAGGTTAGGCCATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCCTCAAGCCTGACGCGCATGCAAAAATCGtactc                                         >  6:250739/1‑145 (MQ=255)                                                                                                               cGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTtctc                                                                              <  2:292018/110‑1 (MQ=255)                                                                                                                 tCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGCGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                              <  1:37743/79‑1 (MQ=255)                                                                                                                 tCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                              >  2:37743/1‑78 (MQ=255)                                                                                                                   aCTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCaa                                    <  3:295898/148‑1 (MQ=255)                                                                                                                   aCTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAg                                   <  7:242741/149‑1 (MQ=255)                                                                                                                   aCTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCGCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGc                                      >  4:274576/1‑146 (MQ=255)                                                                                                                        cAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                              <  8:140009/149‑1 (MQ=255)                                                                                                                        cAGGTTAGGCAATCCCCCGGAAACCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATca                               >  8:160259/1‑148 (MQ=255)                                                                                                                              aGGCATTCCCCCGGAACCCGGGTATACTTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCATCGTACTTTCTTCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGCCAAGATCACCACTc                         <  5:250739/148‑1 (MQ=255)                                                                                                                                gCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTGATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCAGCAAGATCACCACTCGCt                      >  1:93899/1‑149 (MQ=255)                                                                                                                                  aaTCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑Tgg                                                                                                 >  2:227343/1‑72 (MQ=255)                                                                                                                                  aaTCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑Tgg                                                                                                 <  1:227343/72‑1 (MQ=255)                                                                                                                                  aaTCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTc                    >  3:344728/1‑149 (MQ=255)                                                                                                                                      ccccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCg                   >  5:161434/1‑146 (MQ=255)                                                                                                                                      ccccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGACATCGGCAAGATCACCACTCGCTTCgaa                 <  1:199946/148‑1 (MQ=255)                                                                                                                                       cccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGTAAGATCACCACTCGCTTCgaag                >  1:235615/1‑148 (MQ=255)                                                                                                                                       cccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag                >  1:341435/1‑148 (MQ=255)                                                                                                                                        ccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAAGGGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCg                   >  1:107508/1‑146 (MQ=255)                                                                                                                                        ccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaa                 >  6:164863/1‑146 (MQ=255)                                                                                                                                        ccGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaa                 <  5:164863/146‑1 (MQ=255)                                                                                                                                        ccGGAAACCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaagaa              <  7:160259/149‑1 (MQ=255)                                                                                                                                        attaaaCCCGCGAATTATACGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag                >  7:135696/5‑147 (MQ=255)                                                                                                                                         cGGAAACCGCCTTTAATTCGCCG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGaagaag             <  3:250014/149‑1 (MQ=255)                                                                                                                                           aaaCCCGCGAATTATACGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag                <  3:247423/143‑1 (MQ=255)                                                                                                                                           aaaCCCGCGAATTATACGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag                >  4:247423/2‑144 (MQ=255)                                                                                                                                             aCCCGCGAATTATACGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCgg         >  2:82530/1‑149 (MQ=255)                                                                                                                                              cccGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGc        >  5:290363/1‑149 (MQ=255)                                                                                                                                               ccGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGcc       >  7:54706/1‑149 (MQ=255)                                                                                                                                                    tataATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCTAAGAAGCCGGCCTGaaa  <  5:149003/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAA  >  NC_002947/983434‑983729

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