Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,401,795	G→A	R86R (CGC→CGT)	nuoE ←	NADH:ubiquinone oxidoreductase, chain E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,401,795	0	G	A	100.0%	224.4 / NA	61	R86R (CGC→CGT)	nuoE	NADH:ubiquinone oxidoreductase, chain E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/25);  total (36/25)

CCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGGCGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAGC  >  NC_000913/2401648‑2401936                                                                                                                                                    |                                                                                                                                              ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGa                                                                                                                                            <  1:454810/151‑1 (MQ=255)   ccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGAcc                                                                                                                                          <  1:215969/151‑1 (MQ=255)       gcagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAAcc                                                                                                                                      >  2:598260/1‑151 (MQ=255)            aaGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTg                                                                                                                                 >  2:646247/1‑151 (MQ=255)            aaGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTg                                                                                                                                 >  1:525493/1‑151 (MQ=255)                 ggCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGa                                                                                                                            <  1:7402/151‑1 (MQ=255)                   cagcagCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGa                                                                                                                            >  2:454161/1‑149 (MQ=255)                       agcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCt                                                                                                                      >  1:747938/1‑151 (MQ=255)                         cGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGa                                                                                                                    <  2:700926/151‑1 (MQ=255)                          gTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGAc                                                                                                                   >  2:667984/1‑151 (MQ=255)                                 ggCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTTACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGa                                                                                                             <  2:38198/150‑1 (MQ=255)                                 ggCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGa                                                                                                             >  1:558452/1‑150 (MQ=255)                                    cATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGt                                                                                                         >  2:284090/1‑151 (MQ=255)                                     aTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGtt                                                                                                        <  1:509595/151‑1 (MQ=255)                                     aTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGtt                                                                                                        >  1:681429/1‑151 (MQ=255)                                       cAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGc                                                                                                      >  1:771624/1‑151 (MQ=255)                                        aaaTGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTTCTCGAGCGCCGCCTGAATCCCCTGATAACCGGTGATATGACAGAGCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGAAAAACGTTGc                                                                                                      <  1:569389/150‑1 (MQ=255)                                              cGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACAcc                                                                                                <  2:321585/150‑1 (MQ=255)                                                  tGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCg                                                                                            >  2:470309/1‑150 (MQ=255)                                                  tGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCg                                                                                            >  2:853020/1‑150 (MQ=255)                                                  tGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCg                                                                                            <  2:3622/150‑1 (MQ=255)                                                   gCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGAc                                                                                          <  1:338509/151‑1 (MQ=255)                                                        ggTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGc                                                                                     <  2:442115/151‑1 (MQ=255)                                                                   cAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAAt                                                                           >  1:470088/1‑150 (MQ=255)                                                                   cAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAAt                                                                           <  1:771821/150‑1 (MQ=255)                                                                      cTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATAccc                                                                       >  2:449338/1‑151 (MQ=255)                                                                         tttCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAg                                                                     <  2:624398/150‑1 (MQ=255)                                                                           tctcGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGca                                                                   <  2:786199/150‑1 (MQ=255)                                                                              cGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACAt                                                                <  1:852961/150‑1 (MQ=255)                                                                                  cgccgcCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCg                                                           <  1:667046/151‑1 (MQ=255)                                                                                         tGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGt                                                    >  2:320516/1‑151 (MQ=255)                                                                                           aaTACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGt                                                    <  2:764166/149‑1 (MQ=255)                                                                                            aTACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCa                                                                      <  1:307220/130‑1 (MQ=255)                                                                                            aTACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCa                                                                      >  2:307245/1‑130 (MQ=255)                                                                                              aCCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                             >  1:286189/1‑137 (MQ=255)                                                                                              aCCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                             <  2:286214/137‑1 (MQ=255)                                                                                                  tGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCAc                                           >  1:407415/1‑151 (MQ=255)                                                                                                  tGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCAc                                           >  1:465424/1‑151 (MQ=255)                                                                                                   gATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCa                                            <  1:704960/149‑1 (MQ=255)                                                                                                    aTAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCAcc                                          >  2:76458/1‑150 (MQ=255)                                                                                                      aaCCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATc                                       >  2:713904/1‑151 (MQ=255)                                                                                                       aCCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATcc                                      >  2:681873/1‑151 (MQ=255)                                                                                                          gTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGc                                   >  1:764563/1‑151 (MQ=255)                                                                                                              atatGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCAccc                               >  2:173606/1‑151 (MQ=255)                                                                                                               tatGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCAccc                               >  2:174384/1‑150 (MQ=255)                                                                                                                  gACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGcc                           >  2:302337/1‑151 (MQ=255)                                                                                                                        ccACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACgctg                     <  1:371805/151‑1 (MQ=255)                                                                                                                        ccACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACgctg                     >  2:232080/1‑151 (MQ=255)                                                                                                                         cACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACgctgc                    >  2:58673/1‑151 (MQ=255)                                                                                                                         cACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACgctg                     >  1:355773/1‑150 (MQ=255)                                                                                                                              tGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTg               >  2:244885/1‑151 (MQ=255)                                                                                                                              tGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTg               >  2:263509/1‑151 (MQ=255)                                                                                                                               gTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGa              <  2:411205/151‑1 (MQ=255)                                                                                                                                tCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCt                >  1:162466/1‑148 (MQ=255)                                                                                                                                  acaATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTTTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCGCCATCCGGCACCCAGCCACGCTGCTTCTGAACg           <  1:2090/151‑1 (MQ=255)                                                                                                                                  acaATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACg           <  2:452052/151‑1 (MQ=255)                                                                                                                                  acaATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACg           >  2:142181/1‑151 (MQ=255)                                                                                                                                   caATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACg           >  2:739622/1‑150 (MQ=255)                                                                                                                                     aTAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGAtt        >  2:402434/1‑151 (MQ=255)                                                                                                                                          ggATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAg   <  1:675540/151‑1 (MQ=255)                                                                                                                                           gATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAgc  >  2:483433/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                              CCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGGCGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAGC  >  NC_000913/2401648‑2401936

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