Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,530,076:1	(A)7→8	100%	coding (155/483 nt)	ydcD →	putative immunity protein for RhsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,530,069	1	.	A	94.4%	233.9 / 6.7	72	coding (148/483 nt)	ydcD	putative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base A (34/34);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACG  >  NC_000913/1529970‑1530161                                                                                                     |                                                                                               tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                 <  2:1976727/100‑1 (MQ=255) tGTTTATTTGTTATTTACCTTGGTTTTGATAACAATTCAAATTCTGACATCTTATTTTATGGACATAAAACACCATAGATCGTTGAGATATATCTTTCTg                                                                                                 <  2:1869458/100‑1 (MQ=255)     tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aa                                                                                          >  1:344344/1‑100 (MQ=255)      aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaa                                                                                         >  2:910519/1‑100 (MQ=255)       tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaca                                                                                        <  1:2226365/100‑3 (MQ=255)       tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaa                                                                                        >  1:1191501/1‑100 (MQ=255)        ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaaa                                                                                       >  1:247113/1‑100 (MQ=255)        ttgttgTGTACCTTAGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAGGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaa                                                                                        <  1:1020076/99‑1 (MQ=255)          gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaaaat                                                                                     <  1:40689/100‑3 (MQ=255)           ttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaaaata                                                                                    >  1:1781744/1‑97 (MQ=255)            tgttTACCTTGGGTTTGATAACAATTAAAATTCTGACATCGTATTTTATGGACATAAAACACAAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaaaatat                                                                                   >  2:468902/1‑96 (MQ=255)             gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa‑aaaaaatatt                                                                                  >  2:2304853/1‑95 (MQ=255)             gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AA‑AAAAAtatta                                                                                  >  2:1262405/1‑100 (MQ=255)              tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAAtatta                                                                                  >  1:1546327/1‑100 (MQ=255)              tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGAAATAAAACATCAAAGAGAGTTGTAATATATCTTTCTT‑‑aa‑aaaaaatatt                                                                                  >  2:721385/1‑94 (MQ=255)               ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAAtattat                                                                                 <  1:1045384/100‑1 (MQ=255)               ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAAtattat                                                                                 <  1:1512716/100‑1 (MQ=255)               ttACCTTGAGTTTGATAACAATTCAAATTCTGACATCGTATTTTATTGACATAAAACACCAAATAGCTTTGATATATATCTTTCTGA‑AA‑AAAAAtattat                                                                                 <  2:1671003/100‑1 (MQ=255)                  ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AA‑AAAAATAttatttat                                                                             <  1:93081/100‑1 (MQ=255)                   nttGGGTTTGATAACAAATCAAATTCTGACATCGTATTTTATTGACATAAAACACCAAAGAGCATTGAGATATATGTTTCTGA‑AA‑AAAAATAttattta                                                                              >  2:487923/2‑99 (MQ=255)                        gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATtaaa                                                                         >  2:1322267/1‑99 (MQ=255)                           tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAA‑AAAAATATTATTTATAAaataa                                                                      >  2:1696585/1‑100 (MQ=255)                            gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaatt                                                                    <  1:1897480/100‑1 (MQ=255)                            gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaatt                                                                    >  1:240858/1‑100 (MQ=255)                            gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaatt                                                                    >  2:505406/1‑100 (MQ=255)                             aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaatta                                                                   <  2:739420/100‑1 (MQ=255)                             aTAAAAATTCAAATTCTGACATAGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaatt                                                                    <  1:1090406/99‑1 (MQ=255)                              tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaattaa                                                                  >  1:1299149/1‑100 (MQ=255)                               aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaattaat                                                                 >  2:2118911/1‑100 (MQ=255)                               aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAAtaattaat                                                                 >  2:1051629/1‑100 (MQ=255)                                aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATg                                                                <  2:1773098/100‑1 (MQ=255)                                 ccaTTCAAATTCTGACATCGTATTTTATGGACATAAACCACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AACAAAAATATTCTTTATAAGAtaattaa                                                                  >  2:278557/3‑97 (MQ=255)                                 cAATTCAAATTCTGACATAGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTTATAaataatatttt                                                                <  1:1445937/99‑6 (MQ=255)                                  aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGTGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGAc                                                              <  2:1100829/100‑1 (MQ=255)                                  aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGa                                                               >  1:1756993/1‑99 (MQ=255)                                   aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGAcc                                                             >  2:423978/1‑100 (MQ=255)                                    ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCa                                                            >  1:940176/1‑100 (MQ=255)                                       aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCaaaa                                                         >  1:1085200/1‑100 (MQ=255)                                       aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGAACaaaa                                                         >  2:1251708/1‑100 (MQ=255)                                             tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                   <  1:2304880/100‑1 (MQ=255)                                             tGACATCATATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                   <  2:1860855/100‑1 (MQ=255)                                              gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                  <  1:1613694/100‑1 (MQ=255)                                                cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                >  2:1221757/1‑100 (MQ=255)                                                 aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTAAAAAATAATTAATGACCAAAAAATTAGTaga                                               >  2:1423855/1‑100 (MQ=255)                                                    gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                               >  1:705510/1‑97 (MQ=255)                                                     tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATAAATATAA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                            <  2:1863511/99‑1 (MQ=255)                                                            tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                    >  2:994677/1‑100 (MQ=255)                                                            tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                    <  1:505409/100‑1 (MQ=255)                                                             ggACCTACAACACCAAAAAACGATGAGATATATCTTTCTGAAAA‑AAAAATATTATTTATAAAATAATTAATGAACAAAAAATTAGTAAATGGAATGATCa                                    >  2:1077746/1‑100 (MQ=255)                                                             ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAAAGGGAATGGTCa                                    >  1:215994/1‑100 (MQ=255)                                                             ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                   >  2:743810/1‑100 (MQ=255)                                                             ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGAGAATGGTCAt                                   >  2:1002707/1‑100 (MQ=255)                                                              gACATAAAACACCAAAGAGCGATGAGATATAACTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAATAGAGGGAAtaataa                                    >  2:495508/1‑93 (MQ=255)                                                               aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                  >  1:1302194/1‑99 (MQ=255)                                                                 aTAAAACACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTAGTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAGGGTCATTTta                               >  1:1915041/1‑100 (MQ=255)                                                                  tAAAACACTAAAGAGCGTTGAGATATAAAAAAAAAA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                <  2:1893112/98‑1 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAGAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                              >  2:1282120/1‑100 (MQ=255)                                                                   aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                             <  1:271299/100‑1 (MQ=255)                                                                   aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                             <  1:57588/100‑1 (MQ=255)                                                                   aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                             <  2:994743/100‑1 (MQ=255)                                                                     aaCACGAAAGAGCGTTGAGATATATTTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                           <  1:1753923/100‑1 (MQ=255)                                                                     aaCACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                           <  2:681265/100‑1 (MQ=255)                                                                     aaCACCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                           <  1:457300/100‑1 (MQ=255)                                                                     aaCACCAAAGAGCGTTGAGATATAAATAAATGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                           <  2:1592579/100‑1 (MQ=255)                                                                      acacCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                          <  1:160889/100‑1 (MQ=255)                                                                       cacCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                         <  1:1200864/100‑1 (MQ=255)                                                                       cacCAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                         <  1:515249/100‑1 (MQ=255)                                                                        tccAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                         <  2:1428779/98‑1 (MQ=255)                                                                        acCAAAGAGCGTTGAGATATAAAAATAAAA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                          <  2:72121/98‑1 (MQ=255)                                                                         ccAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                       <  1:311341/100‑1 (MQ=255)                                                                         ccAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtaag                       >  1:1428794/1‑98 (MQ=255)                                                                          cAAAGAGCGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                      >  1:531047/1‑100 (MQ=255)                                                                                 cGTTGAGATATATCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa               >  1:332941/1‑100 (MQ=255)                                                                                      agaTATATCTTTATGA‑AA‑AAAAATATTATTTATAAAATAAATAATGAACAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata             <  2:71831/97‑1 (MQ=255)                                                                                      aaatatatTTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGATCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat              <  2:1546402/94‑1 (MQ=255)                                                                                        atatatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt          <  1:269746/98‑1 (MQ=255)                                                                                        atatatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc        <  1:881791/100‑1 (MQ=255)                                                                                         tatatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa       <  2:347686/100‑1 (MQ=255)                                                                                         tatatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa       >  1:1581004/1‑100 (MQ=255)                                                                                           tatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAgg     >  2:2318049/1‑100 (MQ=255)                                                                                           tatCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAgg     >  1:1283193/1‑100 (MQ=255)                                                                                            atCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTATATGGGATGGTTGTTTTATAAGTATAATGGTTAATAATTACAGGa    >  2:1680957/1‑100 (MQ=255)                                                                                             tCTTTCTGAAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGAAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGa    >  2:1456898/1‑100 (MQ=255)                                                                                             tCTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc   >  2:1117938/1‑100 (MQ=255)                                                                                              cTTTCTGA‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg  <  2:939836/100‑1 (MQ=255)                                                                                                     |                                                                                               TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACG  >  NC_000913/1529970‑1530161

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