Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,612,542 C→T 6.8% R122C (CGC→TGC)  PP_4932 → D‑arabinose 1‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,612,5420CT 6.8% 171.8 / 2.7 59R122C (CGC→TGC) PP_4932D‑arabinose 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (2/2);  total (25/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 = 7.82e-01

CAGCGAAGAACGCCTGGGCCCGCTGCTGCGCGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCGGGCTGTCC  >  NC_002947/5612397‑5612682
                                                                                                                                                  |                                                                                                                                            
caGCGAATAACGCCTGGGCCCGCTGCTGCGCGGCCAACGTGATGAATGGGTAATAGTCAGCAAGGTCGGTGAAGAGTTCGAGGAATGCCTGTCTCATTTCGACTTCAGCTCTGCCCACACCCGCCGCACGTTGGAGCGAAGCCTGCGcc                                                                                                                                           >  4:243224/1‑149 (MQ=255)
   cGAAGAACGCCTGGGCCCGCTGCTGCGCGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGc                                                                                                                                        <  3:155663/149‑1 (MQ=255)
                    cGCTGCTGCGAGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATc                                                                                                                       <  3:243224/149‑1 (MQ=255)
                       tgctgcGCGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGTGCCGGCTCGAAACCGACCGCATCGaa                                                                                                                    >  4:222375/1‑149 (MQ=255)
                       tgctgcGCGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGTGCCGGCTCGAAACCGACCGCATCGaa                                                                                                                    >  4:218313/1‑149 (MQ=255)
                            cgcgGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActggt                                                                                                               <  3:50940/149‑1 (MQ=255)
                                gCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActggtgctg                                                                                                           >  2:181470/1‑149 (MQ=255)
                                  cAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTGGGTGGAGCGCAGCCTGCGCCGGCTCGAAAGCGACCGCATCGAActggtgctggt                                                                                                         <  5:296283/149‑1 (MQ=255)
                                  cAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActggtgctggt                                                                                                         <  7:101088/149‑1 (MQ=255)
                                  cAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTACGACTTCAGCGCTGCCCACACCCGCCGCTCTGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActggtgctggt                                                                                                         >  5:104373/1‑149 (MQ=255)
                                       tgatgaATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCAtt                                                                                                    <  6:206499/149‑1 (MQ=255)
                                        gatgaATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTc                                                                                                   >  6:15221/1‑149 (MQ=255)
                                             aTGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACg                                                                                              <  6:121618/149‑1 (MQ=255)
                                                       gTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActg                                                                                                                 <  7:54114/120‑1 (MQ=255)
                                                       gTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActg                                                                                                                 >  8:54114/1‑120 (MQ=255)
                                                       gTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCt                                                                                    >  6:105866/1‑149 (MQ=255)
                                                          aGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGGGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTgg                                                                                  <  7:123939/148‑1 (MQ=255)
                                                             aaGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCAt                                                                              <  1:168754/149‑1 (MQ=255)
                                                              aGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATc                                                                             <  5:15221/149‑1 (MQ=255)
                                                                     tGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActg                                                                                                                 <  7:3514/106‑1 (MQ=255)
                                                                     tGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAActg                                                                                                                 >  8:3514/1‑106 (MQ=255)
                                                                     tGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGaa                                                                       >  6:109958/1‑148 (MQ=255)
                                                                     tGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGaa                                                                       >  1:12210/1‑148 (MQ=255)
                                                                         gagTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAAcagca                                                                  <  1:103889/149‑1 (MQ=255)
                                                                            ttCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAAcagcag                                                                 >  7:175611/1‑147 (MQ=255)
                                                                            ttCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAAcagcag                                                                 >  6:245284/1‑147 (MQ=255)
                                                                              ctaGGATGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACAGCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGAATTCCGACGGCAACGACCTGGCCATCCTCGAAGAGCaggag                                                              <  6:104373/146‑1 (MQ=255)
                                                                              cGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCaggag                                                              <  5:109958/148‑1 (MQ=255)
                                                                                aGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTc                                                           >  8:100680/1‑149 (MQ=255)
                                                                                 ggAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCt                                                          <  7:220302/149‑1 (MQ=255)
                                                                                   aaGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTAc                                                        <  2:45805/149‑1 (MQ=255)
                                                                                        cTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGAc                                                   >  7:153478/1‑149 (MQ=255)
                                                                                        cTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGACTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGAc                                                   >  3:149598/1‑149 (MQ=255)
                                                                                             gCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACg                                                                                              >  1:218630/1‑102 (MQ=255)
                                                                                             gCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACg                                                                                              <  2:218630/101‑1 (MQ=255)
                                                                                                   cGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGaa                                         <  3:252234/148‑1 (MQ=255)
                                                                                                   cGACTTCAGCGCTGCCAACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCAATCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGaa                                         >  4:252234/1‑148 (MQ=255)
                                                                                                      cTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCa                                     <  1:181470/149‑1 (MQ=255)
                                                                                                      cTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCa                                     <  1:268954/149‑1 (MQ=255)
                                                                                                         cAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCTGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCaa                                    <  1:14257/147‑1 (MQ=255)
                                                                                                          aGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGTGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCa                                 <  3:218313/149‑1 (MQ=255)
                                                                                                          aGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGTGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCa                                 <  3:222375/149‑1 (MQ=255)
                                                                                                             gcTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGACTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCagga                              <  4:149598/149‑1 (MQ=255)
                                                                                                                   cacaCCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAAc                                        >  6:76515/1‑133 (MQ=255)
                                                                                                                   cacaCCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAAc                                        <  5:76515/133‑1 (MQ=255)
                                                                                                                    acacCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAAc                                        >  2:172426/1‑132 (MQ=255)
                                                                                                                    acacCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACACCAGGAGGTCTACCAGACCCTGGCCGAAc                                        <  1:172426/132‑1 (MQ=255)
                                                                                                                      acCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGc                          >  4:57200/1‑144 (MQ=255)
                                                                                                                            cgcTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCgg               <  7:100680/149‑1 (MQ=255)
                                                                                                                                  gTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCgg         <  5:245284/149‑1 (MQ=255)
                                                                                                                                  gTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCgg         <  2:12210/149‑1 (MQ=255)
                                                                                                                                    ggAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCaggagg                            <  5:292447/128‑1 (MQ=255)
                                                                                                                                    ggAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCaggagg                            >  6:292447/1‑128 (MQ=255)
                                                                                                                                       gcgcAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCGGGCTGt    >  7:16200/1‑149 (MQ=255)
                                                                                                                                         gcAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCGGGCTGTcc  <  8:153478/149‑1 (MQ=255)
                                                                                                                                              ccTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCg                <  6:227307/131‑1 (MQ=255)
                                                                                                                                              ccTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCg                >  5:227307/1‑131 (MQ=255)
                                                                                                                                              ccTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCg                >  5:223081/1‑131 (MQ=255)
                                                                                                                                              ccTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCg                <  6:223081/131‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CAGCGAAGAACGCCTGGGCCCGCTGCTGCGCGGCCAACGTGATGAATGGGTAATCGTCAGCAAGGTCGGTGAAGAGTTCGAGGAAGGCCTGTCGCATTTCGACTTCAGCGCTGCCCACACCCGCCGCTCGGTGGAGCGCAGCCTGCGCCGGCTCGAAACCGACCGCATCGAACTGGTGCTGGTGCATTCCGACGGCAACGACCTGGCCATCCTCGAACAGCAGGAGGTCTACCAGACCCTGGCCGAACTCAAGCAGGAGGGCAAGATCCTCGGTTTCGGGCTGTCC  >  NC_002947/5612397‑5612682

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