Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,493,473	T→G	A302A (GCT→GCG)	tsgA →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,493,473	0	T	G	98.3%	208.5 / ‑6.8	59	A302A (GCT→GCG)	tsgA	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (28/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACATGGTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCTCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCC  >  NC_000913/3493326‑3493615                                                                                                                                                    |                                                                                                                                               aTACATGGTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTg                                                                                                                                             >  2:494098/1‑151 (MQ=255)     aTGGTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGct                                                                                                                                         <  2:806652/151‑1 (MQ=255)      tGGTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGctt                                                                                                                                        >  2:255699/1‑151 (MQ=255)       ggTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGcttc                                                                                                                                       <  1:49143/151‑1 (MQ=255)       ggTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGctt                                                                                                                                        <  2:626749/150‑1 (MQ=255)           ggCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCttcttc                                                                                                                                    >  2:625708/1‑150 (MQ=255)           ggCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCttcttc                                                                                                                                    >  2:785886/1‑150 (MQ=255)            gCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTtctc                                                                                                                                  >  2:904171/1‑151 (MQ=255)             cATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTcc                                                                                                                                 >  1:262397/1‑151 (MQ=255)                  gggCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAgcgc                                                                                                                            <  1:195953/151‑1 (MQ=255)                   ggCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGcgcg                                                                                                                           <  2:668674/151‑1 (MQ=255)                      gTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGGCAATTTTCGCTCTGGGCTTCTTCTCCAGcgcggt                                                                                                                         >  2:770176/1‑148 (MQ=255)                       ttCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCt                                                                                                                       >  1:826083/1‑151 (MQ=255)                              ttATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATaccac                                                                                                                >  1:105006/1‑151 (MQ=255)                               tATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATaccac                                                                                                                <  1:368348/150‑1 (MQ=255)                                aTTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATAccacca                                                                                                              <  1:880869/151‑1 (MQ=255)                                  tcttcGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACcat                                                                                                             <  1:572156/150‑1 (MQ=255)                                     tcGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACcatcat                                                                                                          <  1:806346/150‑1 (MQ=255)                                     tcGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACcatca                                                                                                           <  2:61096/149‑1 (MQ=255)                                         ttcttTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCAccc                                                                                                     >  1:505939/1‑151 (MQ=255)                                          tcttTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCAccc                                                                                                     >  1:197058/1‑150 (MQ=255)                                             ttGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTgg                                                                                                  <  2:387277/150‑1 (MQ=255)                                             ttGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTg                                                                                                   <  1:868044/149‑1 (MQ=255)                                                  ttGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTcac                                                                                            >  2:836423/1‑151 (MQ=255)                                                    gCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTcaca                                                                                           <  2:748668/150‑1 (MQ=255)                                                     cAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCAcagc                                                                                         <  1:128349/151‑1 (MQ=255)                                                     cAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCCGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCAcagc                                                                                         >  2:56532/1‑151 (MQ=255)                                                       aCGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCAcagcag                                                                                       <  2:857491/151‑1 (MQ=255)                                                       aCGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCAcagcag                                                                                       <  1:682907/151‑1 (MQ=255)                                                                gACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGt                                                                               >  2:203822/1‑150 (MQ=255)                                                                  ccGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTAcc                                                                            <  2:121955/151‑1 (MQ=255)                                                                      aCTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTc                                                                         <  2:368418/150‑1 (MQ=255)                                                                       ctggctggTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGc                                                                       >  2:415049/1‑151 (MQ=255)                                                                        tggctggTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTc                                                                         <  1:359983/148‑1 (MQ=255)                                                                           ctggTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCaa                                                                    >  1:539427/1‑150 (MQ=255)                                                                                   gCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTa                                                           <  2:352985/151‑1 (MQ=255)                                                                                    cTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTa                                                           >  1:434105/1‑150 (MQ=255)                                                                                      gCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACt                                                        <  1:778575/151‑1 (MQ=255)                                                                                        gATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACttt                                                      <  1:481823/151‑1 (MQ=255)                                                                                           tCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTc                                                   >  1:674066/1‑151 (MQ=255)                                                                                               aTGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGa                                               <  1:494111/151‑1 (MQ=255)                                                                                                     gTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCg                                         >  2:139531/1‑151 (MQ=255)                                                                                                       cTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAACACTGGTTAACTTTGTCCTGACCTGCgg                                        >  1:766068/1‑150 (MQ=255)                                                                                                       cTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCgg                                        >  1:454386/1‑150 (MQ=255)                                                                                                               ccGGAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCgg                               >  1:641962/1‑151 (MQ=255)                                                                                                                 ggAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTa                             <  2:795445/151‑1 (MQ=255)                                                                                                                  gAACACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTAc                            <  2:607288/151‑1 (MQ=255)                                                                                                                   aaCACCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTAc                            >  2:631638/1‑150 (MQ=255)                                                                                                                    acacCAGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACt                           >  1:308790/1‑150 (MQ=255)                                                                                                                         aGCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTg                     <  1:785904/151‑1 (MQ=255)                                                                                                                          gCACATATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGa                    <  2:500404/151‑1 (MQ=255)                                                                                                                            acaTATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGAcc                  <  1:904191/151‑1 (MQ=255)                                                                                                                            acaTATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGAcc                  >  1:654579/1‑151 (MQ=255)                                                                                                                             caTATGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCt                 >  2:370687/1‑151 (MQ=255)                                                                                                                              atatGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCtt                >  1:372366/1‑151 (MQ=255)                                                                                                                                atGGCGTGGTCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCTCCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCttt               <  1:378208/150‑1 (MQ=255)                                                                                                                                         tCAATTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCg     >  2:634829/1‑151 (MQ=255)                                                                                                                                           aaTTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGc   <  2:127571/151‑1 (MQ=255)                                                                                                                                            aTTCTCGCGCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGcc  >  2:600172/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                               ATACATGGTCGGCATGTGGGCGTTCAGCTTTATTCTTCGCTTCTTTGATTTGCAACGCATTCTGACCGTACTGGCTGGTCTGGCTGCGATTCTGATGTACGTCTTTAACACCGGAACACCAGCACATATGGCGTGGTCAATTCTCGCTCTGGGCTTCTTCTCCAGCGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCC  >  NC_000913/3493326‑3493615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑