Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,628,160	G→A	6.2%	G118D (GGC→GAC)	era →	GTPase Era

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,628,160	0	G	A	6.2%	208.3 / 4.1	65	G118D (GGC→GAC)	era	GTPase Era
Reads supporting (aligned to +/- strand):  ref base G (40/20);  new base A (2/2);  total (42/23)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGCCCTGAACCGTTACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAGGGGCACAA  >  ACB143/1628026‑1628292                                                                                                                                       |                                                                                                                                     ggCCCTGAACCGTTACATGAACCGCAATGCCTCGGCGGCCCTCTAGGACGTCGACGTGGTTATCTTCGTGGTTGTCCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGc                                                                                                                                     <  2:347482/136‑1 (MQ=255)    ccTGAACCGTTACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCg                                                                                                                                  >  5:7863/1‑136 (MQ=255)         aCCGTTACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGat                                                                                                                             >  5:392945/1‑136 (MQ=255)           cGTTACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGatca                                                                                                                           >  6:128954/1‑136 (MQ=255)            gTTACATGAACCGCAATGCCTCGGCGGCCCTCTAGGACGTCGACGT‑GTTTTCTTCGTGGTTGACCGTACCAAGTTGACCGACGAGGACCAACTGGTTCTGGAGCTCGTGCAGTATGTCACCGGCCCGCTGatcatc                                                                                                                         <  7:22713/136‑1 (MQ=255)             ttACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGatcatc                                                                                                                         <  8:198440/136‑1 (MQ=255)             gaaCATGAACCGCCATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGGGGTTGACCGTACCAAGAGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGatcatc                                                                                                                         <  8:136865/134‑1 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGCCGAGGACCCACTGCTTCTGGAGCGCGTGCCGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  4:79827/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  1:129306/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  5:298211/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  3:330377/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  2:276196/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  7:358561/1‑136 (MQ=255)                cATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATcgcg                                                                                                                      >  1:331090/1‑136 (MQ=255)                        gCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGCCGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGACCCGCTGATCATCGCGGTcaaca                                                                                                               >  3:279121/1‑135 (MQ=255)                        gCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGACCCGCTGATCATCGCGGTcaaca                                                                                                               >  3:273662/1‑135 (MQ=255)                        gCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTAGTGGTTGACCGTACCAAGTGGACCGACGAGGAAAAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTcaaca                                                                                                               >  3:253333/1‑135 (MQ=255)                         cAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAg                                                                                                             <  3:274083/136‑1 (MQ=255)                           aTGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgac                                                                                                           <  3:202367/136‑1 (MQ=255)                            tGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgacc                                                                                                          >  3:186240/1‑136 (MQ=255)                            tGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTAACAAGAGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATAGCTGTCAACACgacc                                                                                                          >  4:85994/1‑136 (MQ=255)                                    cggcCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCAt                                                                                                  >  7:330682/1‑136 (MQ=255)                                       cccTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGa                                                                                               <  1:341924/136‑1 (MQ=255)                                            aaGGACGTCGACGTGGTTATCTTAGTGGTTGAACGTACCAAGTGGCCCGACGAGGACCCACTGGTTGTGGATCGCGTGAAGTATGTCACCGGCCCGCTGATCATAGCGGTCAAAAATAACGAACGCATTGAagcg                                                                                           >  3:36980/1‑133 (MQ=255)                                            aaGGACGTCGACGTGGTTATATTCCTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCCTGCAGTATTTCACCGGCCCGCTGATCATAACGGTCAACAACACCGACCCCAAGGAagac                                                                                           >  2:150966/1‑134 (MQ=255)                                                    cGACGTGGTTATCTTCGTGGTTGACAGTACCAAGTGGACCGACGAGGCCCAACTGGGTCTCGAGCGCGTGCAGAATGTCACCCGCACCCTGATCATCACGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAg                                                                                  >  2:295748/1‑136 (MQ=255)                                                     gACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGc                                                                                 >  7:380782/1‑136 (MQ=255)                                                      aCGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGCCCGACCGCATGGAAGAGAAAGCCGAGCt                                                                                >  2:243584/1‑136 (MQ=255)                                                      aCGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCt                                                                                >  3:16467/1‑136 (MQ=255)                                                        gTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGa                                                                              <  4:330377/136‑1 (MQ=255)                                                                    gTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTg                                                                   >  2:336387/1‑135 (MQ=255)                                                                           aCCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGc                                                            <  2:235911/135‑1 (MQ=255)                                                                             cGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTg                                                          >  3:155155/1‑135 (MQ=255)                                                                                  cAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAagag                                                    <  3:316908/136‑1 (MQ=255)                                                                                         aCCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGc                                             >  6:384901/1‑136 (MQ=255)                                                                                         aCCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGc                                             >  8:54709/1‑136 (MQ=255)                                                                                           cgacgaGGACAAAATGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAAAAAGACCGACCGCATGGAAGAGAAAGACGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGAAACTGcc                                            >  3:74198/1‑135 (MQ=255)                                                                                             acgaGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGCGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGaa                                         >  5:177918/1‑136 (MQ=255)                                                                                                  gACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAGCAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCg                                    >  7:360279/1‑136 (MQ=255)                                                                                                          ggTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCCAGAGCAACTGCCGAACGCCGACGGGATg                            >  1:146921/1‑136 (MQ=255)                                                                                                          ggTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATg                            >  5:8299/1‑136 (MQ=255)                                                                                                               tGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGAt                       >  4:24782/1‑136 (MQ=255)                                                                                                                 gAGCGCGTGCAGTTTGTCACCGGCCCGCTGATCAACGCGGTAAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGAttt                     <  8:330682/136‑1 (MQ=255)                                                                                                                    cgcgTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCg                  <  7:273152/136‑1 (MQ=255)                                                                                                                    cgcgTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCg                  >  6:278743/1‑136 (MQ=255)                                                                                                                    cgcgTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCg                  <  7:267852/136‑1 (MQ=255)                                                                                                                       gTGCATTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCGCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGcac               <  3:144809/136‑1 (MQ=255)                                                                                                                        tGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgaccgacc                                                                                                      >  1:379429/1‑48 (MQ=255)                                                                                                                        tGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgaccgacc                                                                                                      <  2:379429/48‑1 (MQ=255)                                                                                                                        tGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGcaca              <  2:15601/136‑1 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgac                                                                                                           >  1:198887/1‑42 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAgac                                                                                                           <  2:198887/42‑1 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGATCGCCGAGGTGATGCCGATTTCCGCAcag             >  7:10198/1‑136 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCAcag             >  4:314516/1‑136 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCAcag             >  6:52172/1‑136 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCAcag             >  5:3747/1‑136 (MQ=255)                                                                                                                         gCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCAcag             >  4:379739/1‑136 (MQ=255)                                                                                                                            gTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCAcagcag          >  1:331400/1‑136 (MQ=255)                                                                                                                              aTGTCACCGGCCCGCTTATCATCGAGTTCAACAATACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGATGTGATGCCGATTTCCGCACAGCAggg        <  1:185604/136‑1 (MQ=255)                                                                                                                               tatCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAggg        <  5:366915/133‑1 (MQ=255)                                                                                                                               tGTCACCGTCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAggg        <  4:16467/135‑1 (MQ=255)                                                                                                                               tGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAggg        <  3:57000/135‑1 (MQ=255)                                                                                                                               tGTCACCGACCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAggg        <  4:273662/135‑1 (MQ=255)                                                                                                                               tGTCACCGACCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAggg        <  4:279121/135‑1 (MQ=255)                                                                                                                                    ccGGCCCGCTGATCATCGCGGTCAACAAGCCCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCCCCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAGGGGcacac  >  4:245661/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                     GGCCCTGAACCGTTACATGAACCGCAATGCCTCGGCGGCCCTCAAGGACGTCGACGTGGTTATCTTCGTGGTTGACCGTACCAAGTGGACCGACGAGGACCAACTGGTTCTGGAGCGCGTGCAGTATGTCACCGGCCCGCTGATCATCGCGGTCAACAAGACCGACCGCATGGAAGAGAAAGCCGAGCTGATCCCGCACCTGCAGTGGCTGCAAGAGCAACTGCCGAACGCCGAGGTGATGCCGATTTCCGCACAGCAGGGGCACAA  >  ACB143/1628026‑1628292

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