Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,664,592	A→T	11.5%	E29D (GAA→GAT) ‡	dmsD →	twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,664,592	0	A	T	11.5%	190.1 / 13.5	61	E29D (GAA→GAT) ‡	dmsD	twin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (4/3);  total (32/29)
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.83e-01

AATGGCTATTGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTC  >  NC_000913/1664445‑1664736                                                                                                                                                    |                                                                                                                                                  aaTGGCTATTGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAg                                                                                                                                                  <  1:1116957/149‑1 (MQ=255)               gTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCActctct                                                                                                                                       >  1:1524136/1‑140 (MQ=255)                        ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                               <  2:1643171/144‑1 (MQ=255)                        ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                               >  1:1643171/1‑144 (MQ=255)                              cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                 >  1:685642/1‑136 (MQ=255)                              cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                 <  2:685642/136‑1 (MQ=255)                                       ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                    <  3:184090/124‑2 (MQ=255)                                       ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                    >  4:184090/1‑123 (MQ=255)                                          cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                    >  4:352476/1‑120 (MQ=255)                                          cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                    <  3:352476/121‑2 (MQ=255)                                                          aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAc                                                                                        >  4:548376/1‑149 (MQ=255)                                                            gAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCa                                                                                      >  4:20780/1‑149 (MQ=255)                                                                  cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                 <  1:152908/132‑1 (MQ=255)                                                                  cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                 >  2:152908/1‑132 (MQ=255)                                                                   ccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGt                                                                               >  3:1554352/1‑149 (MQ=255)                                                                       ttttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAtt                                                                           <  3:548376/149‑1 (MQ=255)                                                                         ttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAg                                                                         <  4:1554352/149‑1 (MQ=255)                                                                                  aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                            <  1:1089503/137‑1 (MQ=255)                                                                                  aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                            >  2:1089503/1‑137 (MQ=255)                                                                                   gATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGt                                                               <  3:20780/149‑1 (MQ=255)                                                                                             cTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                     >  3:1507179/1‑149 (MQ=255)                                                                                                            tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                 <  2:1387855/106‑1 (MQ=255)                                                                                                            tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                 >  1:1387855/1‑106 (MQ=255)                                                                                                              ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                 <  4:19363/72‑1 (MQ=255)                                                                                                              ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                 >  3:19363/1‑72 (MQ=255)                                                                                                                gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                  <  1:1540282/85‑1 (MQ=255)                                                                                                                gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                  >  2:1540282/1‑85 (MQ=255)                                                                                                                  gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                      <  3:175138/63‑1 (MQ=255)                                                                                                                  gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                      >  4:175138/1‑63 (MQ=255)                                                                                                                    gcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCa                              >  3:1894094/1‑149 (MQ=255)                                                                                                                      gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                            >  1:1549931/1‑149 (MQ=255)                                                                                                                      gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                            >  1:76170/1‑149 (MQ=255)                                                                                                                         gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                         >  2:665110/1‑149 (MQ=255)                                                                                                                         gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                         >  4:564414/1‑149 (MQ=255)                                                                                                                         gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                         >  4:252659/1‑149 (MQ=255)                                                                                                                               ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                <  1:1617916/72‑1 (MQ=255)                                                                                                                               ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                >  2:1617916/1‑72 (MQ=255)                                                                                                                                taCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTg                  <  4:1507179/149‑1 (MQ=255)                                                                                                                                    aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                       <  2:955469/122‑1 (MQ=38)                                                                                                                                    aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                       >  1:955469/3‑124 (MQ=38)                                                                                                                                     tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc             >  1:1517468/1‑149 (MQ=255)                                                                                                                                     tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc             >  1:1018841/1‑149 (MQ=255)                                                                                                                                     tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc             <  1:665110/149‑1 (MQ=255)                                                                                                                                     tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc             <  3:564414/149‑1 (MQ=255)                                                                                                                                       cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                             <  3:942896/99‑1 (MQ=255)                                                                                                                                       cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                             >  4:942896/1‑99 (MQ=255)                                                                                                                                        agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                 >  2:530828/1‑110 (MQ=255)                                                                                                                                        agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                 <  1:530828/110‑1 (MQ=255)                                                                                                                                           gagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt       <  3:252659/149‑1 (MQ=255)                                                                                                                                            agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt      >  3:1487747/1‑149 (MQ=255)                                                                                                                                            agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt      <  4:1894094/149‑1 (MQ=255)                                                                                                                                            agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt      <  2:1549931/149‑1 (MQ=255)                                                                                                                                             gTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgttt     <  2:1517468/149‑1 (MQ=255)                                                                                                                                              tGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGtttgtttg    >  1:194566/1‑149 (MQ=255)                                                                                                                                               gCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTttgtttgt   <  2:76170/149‑1 (MQ=255)                                                                                                                                                cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc  >  3:1579536/1‑149 (MQ=255)                                                                                                                                                cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc  <  2:194566/149‑1 (MQ=255)                                                                                                                                                 ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                     >  2:812168/1‑97 (MQ=255)                                                                                                                                                 ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtgtt                                                     <  1:812168/97‑1 (MQ=255)                                                                                                                                                  gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt       <  4:1062035/142‑1 (MQ=255)                                                                                                                                                  gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt       >  3:1062035/1‑142 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AATGGCTATTGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTC  >  NC_000913/1664445‑1664736

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