Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,905,499	C→T	7.0%	intergenic (‑93/+35)	PP_4313 ← / ← PP_4314	peptidylprolyl isomerase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,905,499	0	C	T	7.0%	163.7 / 2.2	57	intergenic (‑93/+35)	PP_4313/PP_4314	peptidylprolyl isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

GCTTTGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAGGCTTG  >  NC_002947/4905357‑4905647                                                                                                                                                |                                                                                                                                                     gCTTTGCCGTCTTCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACg                                                                                                                                                 >  8:20476/1‑148 (MQ=255)     tGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAAt                                                                                                                                             <  1:145192/148‑1 (MQ=255)       ccGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGc                                                                                                                                           <  4:169907/148‑1 (MQ=255)         gTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTc                                                                                                                                        <  5:19947/149‑1 (MQ=255)              ccTCGACCAGGTCGATGATCTGCAGTCCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCCGTCGGC‑GAATGCCTCCg                                                                                                                                      <  5:287152/146‑1 (MQ=255)                          cGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGTCGACGGAATGCCTcc                                                                                                                                       <  8:76043/133‑1 (MQ=255)                          cGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGTCGACGGAATGCCTcc                                                                                                                                       >  7:76043/1‑133 (MQ=255)                               atCTGCGGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAAt                                                                                                                   <  7:295925/148‑1 (MQ=255)                                  tGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CTCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGc                                                                                                               >  1:241381/1‑149 (MQ=255)                                  tGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACAt                                                                                                                             >  8:224378/1‑135 (MQ=255)                                  tGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACAt                                                                                                                             <  7:224378/135‑1 (MQ=255)                                  tGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑AGCCGACGGAATGCCTCCGTACAGCCATGACGTTGAATCAGc                                                                                                               >  3:241495/1‑149 (MQ=255)                                    cAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGc                                                                                                               <  4:205848/147‑1 (MQ=255)                                          ttGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACAGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGc                                                                                                       >  2:158406/1‑149 (MQ=255)                                             gCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGGAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTc                                                                                                     <  3:213364/148‑1 (MQ=255)                                                 aCGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGaa                                                                                                  <  1:158406/147‑1 (MQ=255)                                                         acTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTgg                                                                                        >  4:41117/1‑149 (MQ=255)                                                           tCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGcc                                                                                      >  7:276281/1‑149 (MQ=255)                                                             cAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCgg                                                                                    <  1:210978/149‑1 (MQ=255)                                                             cAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCg                                                                                     >  4:99483/1‑148 (MQ=255)                                                               aaCGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGa                                                                                   >  3:183114/1‑148 (MQ=255)                                                                aCGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACa                                                                                 <  4:309808/149‑1 (MQ=255)                                                                    tGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCta                                                                             >  2:281625/1‑148 (MQ=255)                                                                     gAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTcc                                                                            >  5:92314/1‑149 (MQ=255)                                                                      aaaaaCAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCt                                                                           >  1:149981/1‑149 (MQ=255)                                                                      aaaaaCAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTAGAACACCTCTTGGCCGGACAGCTCCt                                                                           >  5:197202/1‑149 (MQ=255)                                                                      aaaaaCAAAACCGCCCTTTTCTCAGGAACGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTcc                                                                            <  5:288441/148‑1 (MQ=255)                                                                        aaacaaaacCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTg                                                                         >  8:33688/1‑149 (MQ=255)                                                                          acaaaacCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGt                                                                       <  3:41117/149‑1 (MQ=255)                                                                               acCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGAt                                                                  <  2:241381/149‑1 (MQ=255)                                                                                  gCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTgg                                                               <  3:99483/149‑1 (MQ=255)                                                                                      ttttCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGAAGTTGAATAAGCTGTTATGATCGAAAAACTCTTGGCCGGACAGAACCATGGTATTATTGGCGa                                                            >  8:200131/1‑148 (MQ=255)                                                                                        ttCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAg                                                          <  6:92314/148‑1 (MQ=255)                                                                                             aGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCt                                                                                                      <  2:158547/99‑1 (MQ=255)                                                                                             aGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTAc                                                    >  3:139534/1‑149 (MQ=255)                                                                                               gAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCa                                                        >  8:62197/1‑143 (MQ=255)                                                                                               gAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCa                                                        <  7:62197/143‑1 (MQ=255)                                                                                                 aTGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCAcg                                                <  4:139534/149‑1 (MQ=255)                                                                                                  tGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCAcgc                                               >  1:10522/1‑149 (MQ=255)                                                                                                     aCCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCggg                                            >  4:276751/1‑149 (MQ=255)                                                                                                     aCCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCggg                                            >  6:307169/1‑149 (MQ=255)                                                                                                      ccGTTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGt                                           <  4:80354/149‑1 (MQ=255)                                                                                                       cGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAg                                                          >  6:37470/1‑133 (MQ=255)                                                                                                       cGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAg                                                          <  5:37470/133‑1 (MQ=255)                                                                                                         gTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTtc                                        >  6:14364/1‑149 (MQ=255)                                                                                                            tGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTtctc                                      <  7:20476/148‑1 (MQ=255)                                                                                                             ggCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCtt                                    >  5:45528/1‑149 (MQ=255)                                                                                                              gCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCtt                                    <  8:276281/148‑1 (MQ=255)                                                                                                                aaGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTcga                                 <  3:276751/149‑1 (MQ=255)                                                                                                                 aGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTcgac                                <  5:191181/149‑1 (MQ=255)                                                                                                                       cGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTcgacg                               >  5:214685/1‑144 (MQ=255)                                                                                                                           gAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGc                      <  7:200131/149‑1 (MQ=255)                                                                                                                                   cGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGaa              <  6:154999/149‑1 (MQ=255)                                                                                                                                    ggAGGCATTCCGTCGGC‑GAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGaaa             >  8:193111/1‑149 (MQ=255)                                                                                                                                         cagt‑CGCCGACGGAATGCCTCCGTACAGACATGACCTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCa        <  3:38305/146‑1 (MQ=255)                                                                                                                                         cATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCa        <  6:197202/149‑1 (MQ=255)                                                                                                                                          aTT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAg       <  5:307169/149‑1 (MQ=255)                                                                                                                                                ccGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAGGCTTg  <  4:183114/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GCTTTGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATT‑CGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAGGCTTG  >  NC_002947/4905357‑4905647

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