Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,981,565	T→C	6.9%	S69P (TCG→CCG)	ihfB →	integration host factor subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,981,565	0	T	C	6.9%	164.0 / 2.1	58	S69P (TCG→CCG)	ihfB	integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (30/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.36e-01

CAAGGATGTGGAGCTGGCCATCAAGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCTCAAGCGCGCCCTGGC  >  NC_002947/1981417‑1981708                                                                                                                                                     |                                                                                                                                                 cAAGGATGTGGAGCTGGCCATCAAGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGt                                                                                                                                                  <  1:49354/149‑1 (MQ=255)      aTGTGGAGCTGGCCATCAAGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTc                                                                                                                                             <  4:148503/149‑1 (MQ=255)                    aTCAAGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGc                                                                                                                                >  7:6299/1‑147 (MQ=255)                        aGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCaa                                                                                                                              >  5:333408/1‑145 (MQ=255)                        aGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCaa                                                                                                                              <  6:333408/145‑1 (MQ=255)                         gACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCg                                                                                                                         >  4:109176/1‑149 (MQ=255)                          aCCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGt                                                                                                                        <  3:109176/149‑1 (MQ=255)                          aCCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGt                                                                                                                        <  1:200069/149‑1 (MQ=255)                          aCCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGt                                                                                                                        >  3:81667/1‑149 (MQ=255)                            cATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGc                                                                                                                      >  6:50692/1‑149 (MQ=255)                                 ttGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCAc                                                                                                                 >  6:23053/1‑149 (MQ=255)                                 ttGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGGCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCAc                                                                                                                 >  4:192312/1‑149 (MQ=255)                                    aGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACtt                                                                                                               >  4:293838/1‑148 (MQ=255)                                         aTGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGcc                                                                                                         >  1:320469/1‑149 (MQ=255)                                              acaATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGc                                                                                                                                <  6:72045/121‑1 (MQ=255)                                              acaATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACt                                                                                                                >  8:270056/1‑137 (MQ=255)                                              acaATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACt                                                                                                                <  7:270056/137‑1 (MQ=255)                                              acaATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGc                                                                                                     >  6:153043/1‑148 (MQ=255)                                                 aTGCCTGGCTACCGGCGATCGCATCGAGAGTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAag                                                                                                 <  8:243072/149‑1 (MQ=255)                                                     cTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGtt                                                                                             >  3:165803/1‑149 (MQ=255)                                                       ggCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGc                                                                                           >  5:146404/1‑149 (MQ=255)                                                          tACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGt                                                                                         >  8:93943/1‑148 (MQ=255)                                                                cGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCggg                                                                                  <  5:63094/149‑1 (MQ=255)                                                                 gATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGt                                                                                 <  4:165803/149‑1 (MQ=255)                                                                     gCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAAt                                                                             <  8:6299/149‑1 (MQ=255)                                                                     gCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAAt                                                                             <  5:23053/149‑1 (MQ=255)                                                                      cATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATg                                                                            <  6:146404/149‑1 (MQ=255)                                                                       aTCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGTCCgg                                                                                   <  4:64150/141‑1 (MQ=255)                                                                       aTCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGGCCgg                                                                                   >  3:64150/1‑141 (MQ=255)                                                                        tCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGaa                                                                          <  5:50692/149‑1 (MQ=255)                                                                           agaTTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGAc                                                                       <  7:286958/149‑1 (MQ=255)                                                                             aTTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCaaa                                                                                                  >  6:3320/1‑120 (MQ=255)                                                                             aTTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCaaa                                                                                                  <  5:3320/120‑1 (MQ=255)                                                                             aTTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGAc                                                                       >  5:146297/1‑147 (MQ=255)                                                                             aTTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGa                                                                     <  6:51124/149‑1 (MQ=255)                                                                                 gcgGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCaa                                                                              <  3:238087/136‑1 (MQ=255)                                                                                 gcgGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCaa                                                                              >  4:238087/1‑136 (MQ=255)                                                                                 gcgGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCAt                                                                 <  3:140981/149‑1 (MQ=255)                                                                                  cgGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATg                                                                <  6:146297/149‑1 (MQ=255)                                                                                          gCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGcc                                                        <  3:293838/149‑1 (MQ=255)                                                                                                 ctcGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACAcc                                                  <  4:81667/148‑1 (MQ=255)                                                                                                    gCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTg                                                <  5:153043/147‑1 (MQ=255)                                                                                                     cTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGAtt                                             >  5:53571/1‑149 (MQ=255)                                                                                                        cACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGCCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTcc                                           >  4:110705/1‑148 (MQ=255)                                                                                                        cACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGCCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCca                                          >  4:115395/1‑149 (MQ=255)                                                                                                           tATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCcac                                         >  8:277492/1‑147 (MQ=255)                                                                                                           tATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACaa                                       <  7:93943/149‑1 (MQ=255)                                                                                                                    ccTCGCGTAGGCCGTAACCCGAAGCCCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATc                              >  7:60907/1‑149 (MQ=255)                                                                                                                      tCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCGAAAGGCAAGTTCGTGCCGCACTTCAAGCCCGGAAAAGAGTTGCGTGACCGGGTAAATGaa                                                                          >  1:238929/1‑104 (MQ=255)                                                                                                                      tCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGaa                                                                          <  2:238929/103‑1 (MQ=255)                                                                                                                             ggCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGt                      >  2:208685/1‑148 (MQ=255)                                                                                                                                cGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAAcc                  >  4:286063/1‑149 (MQ=255)                                                                                                                                 gTAACCCGAAGACCGGCCAGCCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCt                 <  3:115395/149‑1 (MQ=255)                                                                                                                                 gTAACCCGAAGACCGGCCAGCCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCt                 <  3:110705/149‑1 (MQ=255)                                                                                                                                    aCCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGAGTCCACAAGGAGCAATCTGATGCGTAACCTc                >  6:188890/1‑147 (MQ=255)                                                                                                                                                gCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCCCACCTGATTCCACAAGGCGCAATCTGAGGCGTAACCTCAAGCGCCCACTggc  >  8:291282/1‑149 (MQ=255)                                                                                                                                                gCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCTCAAGCGCGCCCTggc  >  1:200479/1‑149 (MQ=255)                                                                                                                                                   aGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCTCaa              <  3:128632/134‑1 (MQ=255)                                                                                                                                                   aGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCTCaa              >  4:128632/1‑134 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CAAGGATGTGGAGCTGGCCATCAAGACCATGCTTGAGCAGATGTCACAATGCCTGGCTACCGGCGATCGCATCGAGATTCGCGGTTTTGGCAGCTTCTCGCTGCACTATCGCGCCCCTCGCGTAGGCCGTAACCCGAAGACCGGCCAGTCGGTCAGCCTCG‑AAGGCAAGTTCGTGCCGCACTTCAAGCCCGGCAAAGAGTTGCGTGACCGGGTCAATGAAGACGAGCATGAAGAGGCCCACACCTGATTCCACAAGGAGCAATCTGATGCGTAACCTCAAGCGCGCCCTGGC  >  NC_002947/1981417‑1981708

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