Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,371,919	A→G	T299A (ACC→GCC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,371,919	0	A	G	98.2%	174.0 / NA	56	T299A (ACC→GCC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATC  >  NC_000913/4371784‑4372050                                                                                                                                        |                                                                                                                                    gTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGt                                                                                                                                 <  1:148138/140‑1 (MQ=255) gTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGt                                                                                                                                 <  2:785760/140‑1 (MQ=255) gTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGt                                                                                                                                 <  2:646729/140‑1 (MQ=255)    gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGAt                                                                                                                              >  2:816389/1‑140 (MQ=255)      aGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTTAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGAtct                                                                                                                            >  1:947709/1‑140 (MQ=255)      aGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGAtct                                                                                                                            >  1:62325/1‑140 (MQ=255)       ggCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGAtctc                                                                                                                           >  2:416603/1‑140 (MQ=255)         cGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTg                                                                                                                         >  2:814474/1‑140 (MQ=255)         cGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGAACGTCGTAAAGCTATGCTGCAGGATATCGCAACACTGACTGGCGGTAGCGTGGACTCTg                                                                                                                         >  2:65980/1‑140 (MQ=255)             gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTGAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAaga                                                                                                                     <  2:400630/140‑1 (MQ=255)             gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAaga                                                                                                                     <  2:489651/140‑1 (MQ=255)              cgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAagag                                                                                                                    <  2:713589/140‑1 (MQ=255)                cTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTTAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGAt                                                                                                                  <  2:947709/140‑1 (MQ=255)                cTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGAt                                                                                                                  <  1:384246/140‑1 (MQ=255)                cTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGAt                                                                                                                  <  1:903790/140‑1 (MQ=255)                       ctctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATg                                                                                                           <  2:24143/140‑1 (MQ=255)                           ggTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGc                                                                                                       >  2:508826/1‑140 (MQ=255)                               gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGa                                                                                                   <  1:712140/140‑1 (MQ=255)                                  aaCACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaa                                                                                                >  2:726232/1‑140 (MQ=255)                                      ccATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCa                                                                                            >  2:79015/1‑140 (MQ=255)                                      ccATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCa                                                                                            >  2:750879/1‑140 (MQ=255)                                         tGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAcc                                                                                         >  1:308869/1‑140 (MQ=255)                                         tGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAcc                                                                                         >  1:861833/1‑140 (MQ=255)                                          gCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAccc                                                                                        >  1:598230/1‑140 (MQ=255)                                          gCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAccc                                                                                        >  1:970309/1‑140 (MQ=255)                                          gCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAccc                                                                                        >  1:584000/1‑140 (MQ=255)                                                   cGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGAcc                                                                               >  2:601185/1‑140 (MQ=255)                                                   cGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGAcc                                                                               >  2:388119/1‑140 (MQ=255)                                                     tGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGGTCTCTGAAGAGATCGGGATGGAGCTGGAAAAAACAACCCTGGAAAACCTg                                                                             >  2:242812/1‑140 (MQ=255)                                                        aaGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGt                                                                          <  1:985127/140‑1 (MQ=255)                                                        aaGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGt                                                                          <  2:584000/140‑1 (MQ=255)                                                          gTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAcc                                                                                         >  1:549845/1‑123 (MQ=255)                                                          gTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAcc                                                                                         <  2:549845/123‑1 (MQ=255)                                                              cTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCt                                                                    >  2:6531/1‑140 (MQ=255)                                                                gCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaa                                                                  <  1:816389/140‑1 (MQ=255)                                                                gCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaa                                                                  <  2:178990/140‑1 (MQ=255)                                                                gCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaa                                                                  <  1:814474/140‑1 (MQ=255)                                                                gCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaa                                                                  <  2:321923/140‑1 (MQ=255)                                                                   gTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGCGATCTCTGAAGAGATCGGTATGGAGCTGGGAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACg                                                               >  2:563080/1‑140 (MQ=255)                                                                            ccGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACgtgatcacaa                                                      <  1:130802/140‑8 (MQ=255)                                                                              gggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                           >  1:300259/1‑133 (MQ=255)                                                                              gggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                           <  2:300259/133‑1 (MQ=255)                                                                                     ggCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGa                                             <  2:861833/140‑1 (MQ=255)                                                                                      gCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGac                                            <  1:472696/140‑1 (MQ=255)                                                                                        gATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGacac                                          <  1:79015/140‑1 (MQ=255)                                                                                           cgtcgtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccac                                       <  1:6531/140‑1 (MQ=255)                                                                                             tcgtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccacca                                     <  1:508826/140‑1 (MQ=255)                                                                                             tcgtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccacca                                     <  2:62325/140‑1 (MQ=255)                                                                                                     cTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTatcatc                             <  1:416603/140‑1 (MQ=255)                                                                                                      tATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGa                                             <  1:38312/123‑1 (MQ=255)                                                                                                      tATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGa                                             >  2:38312/1‑123 (MQ=255)                                                                                                             cAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGt                     >  2:555535/1‑140 (MQ=255)                                                                                                                 atatCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAAGGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGt                 <  2:308869/140‑1 (MQ=255)                                                                                                                 atatCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGt                 <  1:726232/140‑1 (MQ=255)                                                                                                                  tatCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAAACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTg                >  2:608870/1‑140 (MQ=255)                                                                                                                     cGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTgaag             >  1:525300/1‑140 (MQ=255)                                                                                                                       cAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTgaagaa           >  2:191034/1‑140 (MQ=255)                                                                                                                          cccTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCt        >  2:165078/1‑140 (MQ=255)                                                                                                                           ccTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTg       <  1:388119/140‑1 (MQ=255)                                                                                                                             tGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCa     >  2:594668/1‑140 (MQ=255)                                                                                                                                cTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATc  <  1:594668/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    GTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATC  >  NC_000913/4371784‑4372050

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