Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,335,862	A→G	D11G (GAC→GGC)	acnA →	aconitate hydratase 1, aconitase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,335,862	0	A	G	100.0%	234.3 / NA	66	D11G (GAC→GGC)	acnA	aconitate hydratase 1, aconitase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/27);  total (39/27)

TCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATGGT  >  NC_000913/1335728‑1335995                                                                                                                                       |                                                                                                                                      tCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACg                                                                                                                                   >  1:1352815/1‑139 (MQ=255)    aTTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTg                                                                                                                                <  2:106522/139‑1 (MQ=255)     tttGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTGACCCTGTAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTCCGAGAAGCCAGTGAGGGCACGTTGc                                                                                                                               >  1:1261183/1‑139 (MQ=255)     tttGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGc                                                                                                                               <  1:142432/139‑1 (MQ=255)                aTCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGAt                                                                                                                    >  1:1069806/1‑139 (MQ=255)                    aaTCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATaaaa                                                                                                                >  1:116523/1‑139 (MQ=255)                           aCGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCa                                                                                                         <  2:1660508/139‑1 (MQ=255)                           aCGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCa                                                                                                         <  1:1386930/139‑1 (MQ=255)                                           aTCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCg                                                                                         >  2:1176805/1‑139 (MQ=255)                                           aTCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCg                                                                                         >  2:607673/1‑139 (MQ=255)                                                tAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgc                                                                                    <  1:607673/139‑1 (MQ=255)                                                 aaTATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgct                                                                                   >  2:414158/1‑139 (MQ=255)                                                 aaTATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgct                                                                                   >  1:943945/1‑139 (MQ=255)                                                      tCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTaa                                                                              >  1:751329/1‑139 (MQ=255)                                                        aCCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAAt                                                                            >  2:571571/1‑139 (MQ=255)                                                        aCCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAAt                                                                            >  2:296180/1‑139 (MQ=255)                                                          ccTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTAcc                                                                                                          <  1:1362663/107‑1 (MQ=255)                                                          ccTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTAcc                                                                                                          >  2:1362663/1‑107 (MQ=255)                                                               agagAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGgca                                                                                                                                     >  1:488283/1‑75 (MQ=255)                                                               agagAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGgca                                                                                                                                     <  2:488283/75‑1 (MQ=255)                                                                   aaTCAGGCCTTCGCAACCCTGTCAATAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGat                                                                 >  1:1011996/1‑139 (MQ=255)                                                                    aTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTggg                                                                     >  2:897427/1‑134 (MQ=255)                                                                    aTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTggg                                                                     <  1:897427/134‑1 (MQ=255)                                                                    aTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGata                                                                >  1:1189615/1‑139 (MQ=255)                                                                    aTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGata                                                                >  1:963551/1‑139 (MQ=255)                                                                        gggCTTCGCAACCCTGTCATTTAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTCGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATCTCAc                                                            <  2:963551/139‑1 (MQ=255)                                                                        gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGc                                                                                               >  1:1259246/1‑104 (MQ=255)                                                                        gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGc                                                                                               <  2:1259246/104‑1 (MQ=255)                                                                        gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                                            >  1:344338/1‑139 (MQ=255)                                                                        gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                                            >  2:1029/1‑139 (MQ=255)                                                                           cTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCa                                                                                                         <  2:987103/91‑1 (MQ=255)                                                                           cTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCa                                                                                                         >  1:987103/1‑91 (MQ=255)                                                                            ttCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGt                                                        <  1:1029/139‑1 (MQ=255)                                                                            ttCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGt                                                        <  2:344338/139‑1 (MQ=255)                                                                             tCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTAcc                                                                                                          >  1:125220/1‑88 (MQ=255)                                                                             tCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTAcc                                                                                                          <  2:125220/88‑1 (MQ=255)                                                                                     cTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAg                                               <  1:571571/139‑1 (MQ=255)                                                                                        tCATTAAGGAGGAGCTATGTCGTCAACACTACAAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGGCTACCCAAGTCa                                            >  2:986419/1‑139 (MQ=255)                                                                                             aaGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCTCTGGGCGATATCAccc                                                          <  1:646051/120‑1 (MQ=255)                                                                                             aaGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAccc                                                          >  2:646051/1‑120 (MQ=255)                                                                                                  ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgctg                                                                                  <  1:693281/91‑1 (MQ=255)                                                                                                  ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgctg                                                                                  >  2:693281/1‑91 (MQ=255)                                                                                                  ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt                                  >  1:899830/1‑139 (MQ=255)                                                                                                  ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt                                  >  1:450959/1‑139 (MQ=255)                                                                                                  ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt                                  >  1:1352624/1‑139 (MQ=255)                                                                                                    agCTATGTCGTCAACCCTACGAGAAGCCAGAAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTg                                >  1:1565897/1‑139 (MQ=255)                                                                                                              tcAACCCTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAAcc                      >  2:878513/1‑139 (MQ=255)                                                                                                                   ccTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTgg                                                                      >  1:1011679/1‑86 (MQ=255)                                                                                                                   ccTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTgg                                                                      <  2:1011679/86‑1 (MQ=255)                                                                                                                   ccTACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACctgctg                 >  2:1487046/1‑139 (MQ=255)                                                                                                                     tACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCtgctgc                >  2:1635466/1‑138 (MQ=255)                                                                                                                     tACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCtgctgc                <  1:1635466/138‑1 (MQ=255)                                                                                                                     tACGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTgcg               >  2:1498587/1‑139 (MQ=255)                                                                                                                      aCGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTgcgc              <  1:1487046/139‑1 (MQ=255)                                                                                                                       cGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACtt                                                                                                             >  1:900049/1‑43 (MQ=255)                                                                                                                       cGAGAAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACtt                                                                                                             <  2:900049/43‑1 (MQ=255)                                                                                                                        gagaAGCCAGTAAGGGCACGTTGCAGGCCAAAGATa                                                                                                                   <  2:182446/36‑1 (MQ=38)                                                                                                                        gagaAGCCAGTAAGGGCACGTTGCAGGCCAAAGATa                                                                                                                   >  1:182446/1‑36 (MQ=255)                                                                                                                         agaAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTgg           <  2:1011996/139‑1 (MQ=255)                                                                                                                          gaAGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGc          <  2:1189615/139‑1 (MQ=255)                                                                                                                            aGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAg        >  2:28602/1‑139 (MQ=255)                                                                                                                            aGCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAg        >  2:326249/1‑139 (MQ=255)                                                                                                                             gCCAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAgg       >  2:127470/1‑139 (MQ=255)                                                                                                                              ccAGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGa      <  1:28602/139‑1 (MQ=255)                                                                                                                                aGTAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATg    >  1:720719/1‑139 (MQ=255)                                                                                                                                  tAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATGGt  <  2:899830/139‑1 (MQ=255)                                                                                                                                  tAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATGGt  <  2:1352624/139‑1 (MQ=255)                                                                                                                                  tAAGGGCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATGGt  <  1:414158/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCTGCTGCGCTGGCAGGATGGT  >  NC_000913/1335728‑1335995

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