Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,584	A→T	8.2%	intergenic (+42/‑47)	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,584	0	A	T	8.2%	176.9 / 5.4	61	intergenic (+42/‑47)	iscX/ndk	cysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base T (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.77e-01

GGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAA  >  NC_002947/983446‑983729                                                                                                                                           |                                                                                                                                                     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                                                                                                                <  1:136265/96‑1 (MQ=255)                                                                                 tCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTg                                                                                                                >  2:136265/1‑96 (MQ=255)                                                                                   gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTg                                                       <  1:78001/149‑1 (MQ=255)                                                                                   gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTg                                                       >  6:149003/1‑149 (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                                                                                                 <  1:227343/72‑1 (MQ=255)                                                                                                                      aaTCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑Tgg                                                                                                 >  2:227343/1‑72 (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                >  4:247423/2‑144 (MQ=255)                                                                                                                               aaaCCCGCGAATTATACGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag                <  3:247423/143‑1 (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)                                                                                                                                           |                                                                                                                                                     GGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGG‑GTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCA‑‑TGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAA  >  NC_002947/983446‑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: ‑