Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,223,378	G→C	10.0%	P390R (CCG→CGG)	PP_1962 ←	site‑specific recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,223,378	0	G	C	10.0%	160.1 / 6.6	58	P390R (CCG→CGG)	PP_1962	site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base C (3/3);  total (28/30)
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.42e-01

AGCAGACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAGAGGGAGGCGTTCAGTC  >  NC_002947/2223240‑2223518                                                                                                                                              |                                                                                                                                              agcagACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCCGCAGCTGCCGAacacaatc                                                                                                                                     <  1:242643/149‑5 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTg                                                                                                                                    <  5:81259/147‑1 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTg                                                                                                                                    >  6:81259/1‑147 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCCGCAGCTGCCGAacacaatcatt                                                                                                                                  <  6:180442/149‑8 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCCGCAGCTGCCGAacacaatca                                                                                                                                    <  1:115936/147‑6 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCCGCAGCTGCCGAacacaatca                                                                                                                                    >  2:115936/1‑142 (MQ=255)    agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGAGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCCGCAGCTGCCGAacacaatcatt                                                                                                                                  >  5:180442/1‑142 (MQ=255)     gACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGt                                                                                                                                  >  2:193761/1‑148 (MQ=255)      aCGGATCAGTTCAATATTTCTTTCGCTCCGATTGATCAGGATTTCTTCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTa                                                                                                                                <  2:340669/149‑1 (MQ=255)         gATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCt                                                                                                                             <  3:173765/149‑1 (MQ=255)               gTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCTATGCGGATGGAGGCGCTCGATCACCGTCAC‑‑TTCTTAGAGATAATCAATGA‑TTGTGTTCGGAAGCTGCGGAAAACCTTTGGTTAACTTGATc                                                                                                                        >  5:317414/1‑148 (MQ=255)               gTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGAAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATc                                                                                                                        >  6:380650/1‑148 (MQ=255)                tCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACATTTGGTTACCTTGATCga                                                                                                                      >  4:360974/1‑149 (MQ=255)                      tttCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGatgatg                                                                                                                 >  3:42649/1‑148 (MQ=255)                      cttcgttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGatgatg                                                                                                                 >  3:21479/2‑148 (MQ=255)                         cgttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGACTTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGaa                                                                                                               <  1:94186/148‑1 (MQ=255)                           ttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCAc                                                                                                           >  4:91527/1‑149 (MQ=255)                                    aTTGATCAGGATTTGTGCAATCGTTTTCCGAAAGCGTTGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTAGAGATAATCAATGA‑TTGTTTTCGGCAGCTGCGGAACACCTTTGTTTACCTTGATCGATGATGAAGCACT‑CTGaaaaa                                                                                                 <  3:360974/149‑1 (MQ=255)                                    aTTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGaaaaa                                                                                                 <  2:300651/149‑1 (MQ=255)                                    aTTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGaaaaa                                                                                                 <  8:62070/149‑1 (MQ=255)                                         tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCACAGGTCTTCGAGAAAAACAATGA‑TTTTTTACGGCAGATGCGGAACACATTTGTTTACCTTGATCCATAATGAAGAAAT‑ATGAAAAAtc                                                                                               >  1:109049/1‑147 (MQ=255)                                         tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTGGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCAGT‑CTGAAAAATGaa                                                                                             <  4:42649/148‑1 (MQ=255)                                         tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATg                                                                                               >  8:243807/1‑146 (MQ=255)                                         tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGaa                                                                                             <  7:243807/148‑1 (MQ=255)                                         tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAg                                                                                            >  1:112349/1‑149 (MQ=255)                                          gatgATTTCTGCAATTGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGTTCCCCGTCTC‑‑TTCTTCGAGATTATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATTGATGATGTAGCACT‑CTGAAAAATGAAg                                                                                            <  2:109049/145‑1 (MQ=255)                                           aGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGcc                                                                                          >  2:62301/1‑149 (MQ=255)                                             gATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGccct                                                                                        >  1:242720/1‑149 (MQ=255)                                              aTTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGccctc                                                                                       <  2:242720/149‑1 (MQ=255)                                                   tGCAATCGTTTTCCGAAAGCGATGCGTATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGATATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGt                                                                                  <  3:91527/149‑1 (MQ=255)                                                     cAATCGTTTTCCGAAAGCGATGCCGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCTAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTa                                                                                 <  6:364697/148‑1 (MQ=255)                                                      aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGaaaaa                                                                                                 <  4:153671/131‑1 (MQ=255)                                                      aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGaaaaa                                                                                                 >  3:153671/1‑131 (MQ=255)                                                      aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAg                                                                                >  3:128637/1‑148 (MQ=255)                                                        tCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAg                                                                             <  5:291623/149‑1 (MQ=255)                                                           ttttCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTc                                                                          >  3:307806/1‑149 (MQ=255)                                                           ttttCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTc                                                                          >  8:347894/1‑149 (MQ=255)                                                                 gAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACg                                                                      >  6:295433/1‑147 (MQ=255)                                                                  aaaGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTa                                                                    <  1:62301/149‑1 (MQ=255)                                                                  aaaGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACg                                                                      >  2:162069/1‑146 (MQ=255)                                                                           cGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAgagcggagcg                                                          <  3:295199/149‑1 (MQ=255)                                                                                ggAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTtct                                                      <  5:295433/148‑1 (MQ=255)                                                                                        tCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGaa                                                                                             >  4:221567/1‑101 (MQ=255)                                                                                        tCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGaa                                                                                             <  3:221567/101‑1 (MQ=255)                                                                                          gATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTaa                                           >  3:236412/1‑149 (MQ=255)                                                                                                    aC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATTAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCAc                                  <  5:380650/148‑1 (MQ=255)                                                                                                    aC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCAc                                  <  6:317414/148‑1 (MQ=255)                                                                                                                gaTAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTa                       <  4:21479/149‑1 (MQ=255)                                                                                                                 aTAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAg                      <  4:128637/149‑1 (MQ=255)                                                                                                                   aaTCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAgg                     >  6:58903/1‑148 (MQ=255)                                                                                                                    aTCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTtct                                                      <  4:286856/114‑1 (MQ=255)                                                                                                                    aTCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTa                   <  2:112349/149‑1 (MQ=255)                                                                                                                          gA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg               >  7:250136/1‑147 (MQ=255)                                                                                                                          gA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg               >  2:301552/1‑147 (MQ=255)                                                                                                                          gA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg               >  7:225470/1‑147 (MQ=255)                                                                                                                           aaaggtgtTCCGCAGCTGCCGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTCTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg               >  2:288916/5‑147 (MQ=255)                                                                                                                             ttGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAgagcgga                                                             <  6:362560/99‑1 (MQ=255)                                                                                                                                      gCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAGAGGGAGGCGTTCAGtc  <  1:193761/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              AGCAGACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGA‑TTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACT‑CTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAGAGGGAGGCGTTCAGTC  >  NC_002947/2223240‑2223518

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