Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,581	T→A	6.9%	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	6.9%	174.3 / 3.0	58	intergenic (+39/‑50)	iscX/ndk	cysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (32/22);  new base A (2/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

ACCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGC  >  NC_002947/983433‑983723                                                                                                                                                      |                                                                                                                                               aCCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGt                                                                                                                                                <  8:284903/149‑1 (MQ=255)   cTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGta                                                                                                                                               >  5:122686/1‑148 (MQ=255)     gCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATaa                                                                                                                                            >  6:137180/1‑149 (MQ=255)      cACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGtata                                                                                                                                             <  3:310564/147‑1 (MQ=255)      cACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGaatta                                                                                                                                            >  3:252489/1‑146 (MQ=255)              ggTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                     <  7:198376/147‑1 (MQ=255)               gTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGtt                                                                                                                                  >  2:235325/1‑149 (MQ=255)                tGCTGGCACTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttc                                                                                                                                 <  7:251140/149‑2 (MQ=255)                  cTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTaa                                                                                                                               <  6:196906/149‑1 (MQ=255)                     gCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                              <  4:252489/147‑5 (MQ=255)                       aTTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAAttttt                                                                                                                          >  3:152646/1‑149 (MQ=255)                        ttGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                         <  1:235325/149‑1 (MQ=255)                        ttGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                         <  2:11886/149‑1 (MQ=255)                           ccAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCTGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACATGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGt                                                                                                                       <  7:100572/148‑1 (MQ=255)                                    aGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCata                                                                                                             >  3:312202/1‑149 (MQ=255)                                    aGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCata                                                                                                             >  1:324710/1‑149 (MQ=255)                                     gCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCatat                                                                                                            >  7:189007/1‑149 (MQ=255)                                       gacgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTc                                                                                                          <  6:122686/149‑1 (MQ=255)                                         cgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTg                                                                                                        >  7:105174/1‑149 (MQ=255)                                           aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                      <  5:137180/149‑1 (MQ=255)                                                      ggcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                           <  8:105174/149‑1 (MQ=255)                                                      ggcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                           <  4:152646/149‑1 (MQ=255)                                                                aaGTTCTCGATGCAATCCAGGCCGGGTGGATCGACGAAGCCGACTAAGCGTCACTTGGAGGTT‑AGGCAATCCCCCGGAACCCGCGTATTATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAAc                                                                                  <  2:324710/148‑1 (MQ=255)                                                                      tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACtt                                                                            >  2:327032/1‑148 (MQ=255)                                                                        gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTGAGGCATTCCCCCGGAACCCGCGTTAAATAGCCGGGTTTAATTTTTGGTAACAATCTTATTATGTAGTTTACCATGGATGTTAAAGGTCATTtc                                                                          >  1:267748/1‑149 (MQ=255)                                                                        gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                           >  6:109944/1‑131 (MQ=255)                                                                        gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                           <  5:109944/131‑1 (MQ=255)                                                                             aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                     >  4:65540/1‑148 (MQ=255)                                                                             aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTCAATTTTTCGTAACACTCATAtt                                                                                                           <  8:280361/110‑1 (MQ=255)                                                                             aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTCAATTTTTCGTAACACTCATAtt                                                                                                           >  7:280361/1‑110 (MQ=255)                                                                                 cAGGCCGCCTGGATCGACGAAGCCGAATAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                 >  2:211447/1‑148 (MQ=255)                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                 >  7:208253/1‑147 (MQ=255)                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                 <  2:73546/147‑1 (MQ=255)                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCaa                                                               >  1:73546/1‑149 (MQ=255)                                                                                    gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGc                                                             >  8:277681/1‑149 (MQ=255)                                                                                                gacgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                  <  1:164824/148‑1 (MQ=255)                                                                                                gacgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                  >  2:164824/1‑148 (MQ=255)                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCTGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTTTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGTTGTTGc                                                >  5:243673/1‑149 (MQ=255)                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                >  3:199259/1‑149 (MQ=255)                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                >  4:324408/1‑149 (MQ=255)                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                <  4:199259/149‑1 (MQ=255)                                                                                                          aCTAAGCGTCACTTGCAGGTT‑ATGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                        >  7:213874/1‑148 (MQ=255)                                                                                                              aGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACTCTGTTGCCAAAAACGTCAt                                    <  1:211447/148‑1 (MQ=255)                                                                                                              aGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCAt                                    <  3:65540/148‑1 (MQ=255)                                                                                                                cGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCgg                                 >  7:81567/1‑149 (MQ=255)                                                                                                                 gTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGCTGCCAAAAACGTCATCgg                                 >  1:285879/1‑148 (MQ=255)                                                                                                                      ttGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGa                            >  6:249067/1‑148 (MQ=255)                                                                                                                        gCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCgg                                 >  6:10468/1‑141 (MQ=255)                                                                                                                        gCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCgg                                 <  5:10468/141‑1 (MQ=255)                                                                                                                          aGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                        >  1:18149/1‑148 (MQ=255)                                                                                                                                aGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTc                   <  2:267748/148‑1 (MQ=255)                                                                                                                                  gCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                         <  6:284318/92‑1 (MQ=255)                                                                                                                                  gCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                         >  5:284318/1‑92 (MQ=255)                                                                                                                                  gCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCt                >  6:354277/1‑149 (MQ=255)                                                                                                                                  gCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCt                >  5:236979/1‑149 (MQ=255)                                                                                                                                    aaTCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTc              <  5:249067/149‑1 (MQ=255)                                                                                                                                      tCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCg             >  1:140761/1‑148 (MQ=255)                                                                                                                                         cccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                        <  7:194938/134‑1 (MQ=255)                                                                                                                                         cccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                        >  8:194938/1‑134 (MQ=255)                                                                                                                                          ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag          >  2:82359/1‑147 (MQ=255)                                                                                                                                             aaaCCCGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGcc     >  8:150194/2‑149 (MQ=255)                                                                                                                                                cccGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATAACAACTCGCTTCGAAGAAGCCGGc  >  7:98756/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               ACCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTT‑AGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGC  >  NC_002947/983433‑983723

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