Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,634,133	A→C	6.4%	M22L (ATG→CTG)	rnc →	ribonuclease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,634,133	0	A	C	6.4%	191.5 / 3.2	62	M22L (ATG→CTG)	rnc	ribonuclease 3
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base C (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

ACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTCGCGGCTGC  >  NC_002947/1633989‑1634280                                                                                                                                                 |                                                                                                                                                    acaTCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCt                                                                                                                                                 >  3:70878/1‑149 (MQ=255)   aTCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTcc                                                                                                                                               >  4:47337/1‑149 (MQ=255)    tCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCt                                                                                                                                              <  4:159430/149‑1 (MQ=255)    tCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCt                                                                                                                                              <  3:47337/149‑1 (MQ=255)      aaTCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTg                                                                                                                                             <  2:218374/148‑1 (MQ=255)       aTCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATg                                                                                                                                                   <  1:186304/141‑1 (MQ=255)       aTCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATg                                                                                                                                                   >  2:186304/1‑141 (MQ=255)                  acaGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTAGACCTTCAAGAACCAGGAGCATCTGGTCCTggttc                                                                                                                                         >  8:130288/1‑137 (MQ=255)                  acaGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTAGACCTTCAAGAACCAGGAGCATCTGGTCCTggttc                                                                                                                                         <  7:130288/140‑4 (MQ=255)                          aCAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTgg                                                                                                                                            >  5:247894/1‑129 (MQ=255)                          aCAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTgg                                                                                                                                            <  6:247894/129‑1 (MQ=255)                          aCAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTAc                                                                                                                        >  2:274579/1‑149 (MQ=255)                              tCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCg                                                                                                                    <  4:108162/149‑1 (MQ=255)                                tttCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCggg                                                                                                                  >  4:181290/1‑149 (MQ=255)                                      gATGTTGATTTGATCAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGAccc                                                                                                                                   >  8:245609/1‑126 (MQ=255)                                      gATGTTGATTTGATCAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGAccc                                                                                                                                   <  7:245609/126‑1 (MQ=255)                                        tGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGcaacaa                                                                                                          >  2:159388/1‑149 (MQ=255)                                         gTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCa                                                                                                                             <  2:25231/129‑1 (MQ=255)                                         gTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCa                                                                                                                             >  1:25231/1‑129 (MQ=255)                                         gTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAAt                                                                                                         >  4:60410/1‑149 (MQ=255)                                           tGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAAt                                                                                                         >  5:313271/1‑147 (MQ=255)                                                 gAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCt                                                                                                 >  6:112339/1‑149 (MQ=255)                                                  aacaacGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTg                                                                                                >  4:319837/1‑149 (MQ=255)                                                    caacGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTggt                                                                                              <  6:191751/149‑2 (MQ=255)                                                     aacGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAg                                                                                             <  3:181290/149‑1 (MQ=255)                                                        gcgTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTc                                                                                          >  5:28614/1‑149 (MQ=255)                                                             gAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCg                                                                                      <  6:313271/148‑1 (MQ=255)                                                              aacaacGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTcc                                                                                         <  1:260403/144‑1 (MQ=255)                                                              aacaacGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTcc                                                                                         >  2:260403/1‑144 (MQ=255)                                                              aacaacGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGt                                                                                    <  3:84351/149‑1 (MQ=255)                                                                        gggTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCAtt                                                                           <  3:319837/148‑1 (MQ=255)                                                                         ggTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCTACTGGAGTTCCTCGGTTACGCCAttat                                                                         >  8:177892/1‑147 (MQ=255)                                                                             gCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCaa                                                                      <  1:159388/148‑1 (MQ=255)                                                                                gCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTc                                                                  >  5:79623/1‑149 (MQ=255)                                                                                gCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTc                                                                  >  6:88382/1‑149 (MQ=255)                                                                                gCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTc                                                                  >  2:119906/1‑149 (MQ=255)                                                                                 cATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCg                                                                 >  1:63383/1‑149 (MQ=255)                                                                                  aTGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGt                                                                >  6:278729/1‑149 (MQ=255)                                                                                     aGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGc                                                             <  3:4574/149‑1 (MQ=255)                                                                                      gTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGcc                                                            <  5:278729/149‑1 (MQ=255)                                                                                        gCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCgg                                                          >  1:64433/1‑149 (MQ=255)                                                                                         cTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATc                                                                                                                                >  2:200258/1‑78 (MQ=255)                                                                                         cTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATc                                                                                                                                <  1:200258/78‑1 (MQ=255)                                                                                                 gCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCg                                                                                      >  4:46781/1‑112 (MQ=255)                                                                                                 gCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCg                                                                                      <  3:46781/112‑1 (MQ=255)                                                                                                  ccGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCt                                                 <  2:63383/148‑1 (MQ=255)                                                                                                   cGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGt                                               <  1:274579/149‑1 (MQ=255)                                                                                                          gAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGc                                         <  5:88382/148‑1 (MQ=255)                                                                                                             cGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGttt                                     >  5:228249/1‑149 (MQ=255)                                                                                                               aaaGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTc                                   >  3:314215/1‑149 (MQ=255)                                                                                                                          acacCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCcg                         <  7:177892/148‑1 (MQ=255)                                                                                                                             ccTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGaa                     <  6:79623/149‑1 (MQ=255)                                                                                                                                    gAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAgg                             <  8:182485/134‑1 (MQ=255)                                                                                                                                    gAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAgg                             >  7:182485/1‑134 (MQ=255)                                                                                                                                      gccAGGAGCATCTGGTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTc                                             >  7:70661/2‑116 (MQ=255)                                                                                                                                      gccAGGAGCATCTGGTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTc                                             <  8:70661/115‑1 (MQ=255)                                                                                                                                        cAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTc          >  4:263136/1‑149 (MQ=255)                                                                                                                                          ggACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTcgc        >  6:312514/1‑149 (MQ=255)                                                                                                                                              cAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTcg         <  6:36774/144‑1 (MQ=255)                                                                                                                                              cAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTcg         >  5:36774/1‑144 (MQ=255)                                                                                                                                                gATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTCGCGGCtgc  <  6:28614/149‑1 (MQ=255)                                                                                                                                                gATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTCGCGGCagc  >  5:42209/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGTGCTTCCCTTGCCGGTCTGGAGCGAAAGCTCGGTTACACCTTCAAGAACCAGGACCAGATGCTCCTGGCACTGACCCATCGCAGTTACGCCGGGCGCAACAATGAGCGCCTGGAGTTCCTCGGTGACGCCATTCTCAACTTCGTTGCCGGCGAAGCGCTGTTCGAGCGTTTTCCGCAGGCCCGCGAAGGCCAGTTGTCGCGGCTGC  >  NC_002947/1633989‑1634280

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