Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	40,676	G→A	T419I (ACC→ATC)	caiT ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,676	0	G	A	98.5%	227.7 / ‑4.2	66	T419I (ACC→ATC)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (32/33);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTTATCT  >  NC_000913/40546‑40810                                                                                                                                   |                                                                                                                                       aaCAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATa                                                                                                                                <  1:75780/139‑1 (MQ=255)           aTAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGa                                                                                                                      >  2:470988/1‑139 (MQ=255)           aTAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGa                                                                                                                      >  2:11777/1‑139 (MQ=255)             aaTGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATg                                                                                                                    >  2:1581737/1‑139 (MQ=255)                gCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATgaag                                                                                                                 >  1:425058/1‑139 (MQ=255)                          aTTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCcaca                                                                                                       >  2:1064808/1‑139 (MQ=255)                                cAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtggtgg                                                                                                 <  1:823302/139‑1 (MQ=255)                                   ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtgg                                                                                                    <  2:986426/133‑1 (MQ=255)                                   ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtgg                                                                                                    >  1:986426/1‑133 (MQ=255)                                   ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              <  2:1462766/139‑1 (MQ=255)                                   ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              >  2:743652/1‑139 (MQ=255)                                       tACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCt                                                                                          <  1:147914/139‑1 (MQ=255)                                       tACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCt                                                                                          <  1:1082014/139‑1 (MQ=255)                                        aCGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTg                                                                                         >  1:597754/1‑139 (MQ=255)                                            aCCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTg                                                                                     >  2:499142/1‑139 (MQ=255)                                             ccAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTgg                                                                                    >  2:565928/1‑139 (MQ=255)                                                gcagAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCag                                                                                 <  1:270133/139‑1 (MQ=255)                                                gcagAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCag                                                                                 >  1:1318330/1‑139 (MQ=255)                                                gcagAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCag                                                                                 >  1:1034389/1‑139 (MQ=255)                                                 cagAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCaga                                                                                >  2:1236335/1‑139 (MQ=255)                                                  agagGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCagag                                                                               <  2:425058/139‑1 (MQ=255)                                                      gtggtTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGc                                                                           >  1:1256048/1‑139 (MQ=255)                                                         gtTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCa                                                                        >  2:823694/1‑139 (MQ=255)                                                         gtTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCa                                                                        <  2:478078/139‑1 (MQ=255)                                                           tcttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTCTGAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCg                                                                                               <  1:961264/114‑1 (MQ=255)                                                           tcttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCag                                                                                                                          >  1:584034/1‑87 (MQ=255)                                                           tcttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCag                                                                                                                          <  2:584034/87‑1 (MQ=255)                                                           tcttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCg                                                                                               >  2:961264/1‑114 (MQ=255)                                                                       cgTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATgg                                                          >  2:740858/1‑139 (MQ=255)                                                                         tACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCGATGATGgcg                                                        <  2:605240/139‑1 (MQ=255)                                                                         tACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcg                                                        <  2:597754/139‑1 (MQ=255)                                                                           cTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcg                                                      >  1:1649837/1‑139 (MQ=255)                                                                             tCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcgcg                                                    <  2:523381/139‑1 (MQ=255)                                                                             tCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcgcg                                                    <  1:11777/139‑1 (MQ=255)                                                                             tCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcgcg                                                    <  1:746827/139‑1 (MQ=255)                                                                             tCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcgcg                                                    <  2:167186/139‑1 (MQ=255)                                                                             tCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcgcg                                                    <  2:1624190/139‑1 (MQ=255)                                                                                 ggCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCag                                                                                 >  1:373261/1‑106 (MQ=255)                                                                                 ggCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCag                                                                                 <  2:373261/106‑1 (MQ=255)                                                                                  gCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCcaca                                               <  2:985359/139‑1 (MQ=255)                                                                                  gCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCcaca                                               <  1:1064808/139‑1 (MQ=255)                                                                                         gACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACt                                        <  2:1649837/139‑1 (MQ=255)                                                                                               gCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGt                                                                     <  2:1351004/104‑1 (MQ=255)                                                                                               gCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGt                                                                     >  1:1351004/1‑104 (MQ=255)                                                                                                ccAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGAt                                 <  2:319537/139‑1 (MQ=255)                                                                                                        taAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCa                                                                        <  2:1096817/92‑1 (MQ=255)                                                                                                        taAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCa                                                                        >  1:1096817/1‑92 (MQ=255)                                                                                                          agagCAGGCGTTAACCAGCGTAACGATGGCAATACAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGa                       <  2:1057130/139‑1 (MQ=255)                                                                                                          agagCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGa                       <  1:823694/139‑1 (MQ=255)                                                                                                              cAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcg                                                        <  2:1586895/102‑1 (MQ=255)                                                                                                              cAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcg                                                        >  1:1586895/1‑102 (MQ=255)                                                                                                              cAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGt                   >  1:883663/1‑139 (MQ=255)                                                                                                               aGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATgaa                                                                                                                  >  1:1062071/1‑43 (MQ=255)                                                                                                               aGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATgaa                                                                                                                  <  2:1062071/43‑1 (MQ=255)                                                                                                                ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtg                 >  2:1628300/1‑139 (MQ=255)                                                                                                                ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtg                 >  2:657378/1‑139 (MQ=255)                                                                                                                ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtg                 >  1:97303/1‑139 (MQ=255)                                                                                                                      aaCCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtgatga             <  1:1431756/137‑1 (MQ=255)                                                                                                                      aaCCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtgatga             >  2:1431756/1‑137 (MQ=255)                                                                                                                           gCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGttttt      >  1:1291997/1‑139 (MQ=255)                                                                                                                            cGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTta     >  2:1056068/1‑139 (MQ=255)                                                                                                                             gTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTtat    <  1:743652/139‑1 (MQ=255)                                                                                                                             gTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTtat    <  2:97303/139‑1 (MQ=255)                                                                                                                              tAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCt                                                                                          >  1:1431810/1‑52 (MQ=255)                                                                                                                              tAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCt                                                                                          <  2:1431810/52‑1 (MQ=255)                                                                                                                               aaCGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTtatct  <  1:452106/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       AACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTTATCT  >  NC_000913/40546‑40810

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