Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,511 C→G 10.5% intergenic (+22/‑29) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5110CG10.5% 141.7 / 7.2 57intergenic (+22/‑29)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base G (3/3);  total (28/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.20e-01

CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCG  >  NC_000913/3020363‑3020658
                                                                                                                                                    |                                                                                                                                                   
caGCGGGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGc                                                                                                                                                     <  1:389360/149‑1 (MQ=255)
         ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                            >  6:15078/1‑149 (MQ=255)
          ggcAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                           <  1:390015/147‑1 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAATGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                           <  8:15848/149‑1 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                           <  5:309140/149‑1 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                            <  7:219110/148‑1 (MQ=255)
            cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGt                                                                                                                                         >  5:341047/1‑149 (MQ=255)
             aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGtt                                                                                                                                        <  1:131004/149‑1 (MQ=255)
              aTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGtt                                                                                                                                        >  1:268634/1‑148 (MQ=255)
                 gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                                                                                                                                    <  2:268634/149‑1 (MQ=255)
                    tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATaa                                                                                                                                  <  8:364377/148‑1 (MQ=255)
                        ggTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAggg                                                                                                                                          >  7:343528/1‑136 (MQ=255)
                        ggTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAggg                                                                                                                                          <  8:343528/136‑1 (MQ=255)
                        ggTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGa                                                                                                                               >  7:168585/1‑147 (MQ=255)
                                 aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                    <  3:92892/149‑1 (MQ=255)
                                     cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                >  2:166275/1‑149 (MQ=255)
                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                           >  3:62780/1‑121 (MQ=255)
                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                           <  4:62780/121‑1 (MQ=255)
                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                >  3:26358/1‑148 (MQ=255)
                                           tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                          >  5:366456/1‑149 (MQ=255)
                                             tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                        >  1:44805/1‑149 (MQ=255)
                                             tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtttt                                                                                                         <  1:140960/148‑1 (MQ=255)
                                               aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                   <  4:54499/136‑1 (MQ=255)
                                               aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                       <  6:243172/148‑1 (MQ=255)
                                               aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                      <  5:15078/149‑1 (MQ=255)
                                                         cGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTCTGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGAct                                                                                             <  4:335344/148‑1 (MQ=255)
                                                           aTTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTg                                                                                          <  5:320482/149‑1 (MQ=255)
                                                              ccATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                        <  6:168021/148‑1 (MQ=255)
                                                                 tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                      <  4:269421/131‑1 (MQ=255)
                                                                 tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                      >  3:269421/1‑131 (MQ=255)
                                                                        acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTTGCATAAATGACAGATGCCGTCTTCGCGCAGGGAAGAGGGTTAAATAAGGATAGAAGATGATTTAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                             <  8:271815/149‑1 (MQ=255)
                                                                         caATCCAGAAACTCCGCTGCCAGTATGTGGCGGCGGATGGATGCGCTGGCATAAATGAAATATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAATAATATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                            <  6:41738/149‑1 (MQ=255)
                                                                          aatccAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCAGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAATATGGTTAAATTAGGAAAGAAGATGATTTAACAATATTTCAGGCCCGACACTGACGCACACTCGCTg                                                                            >  8:356757/4‑148 (MQ=255)
                                                                                 aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                    <  2:392843/149‑1 (MQ=255)
                                                                                   aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                  <  1:166275/149‑1 (MQ=255)
                                                                                       cGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgccg                                                              >  1:371100/1‑149 (MQ=255)
                                                                                        gCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgccgc                                                             >  6:299150/1‑149 (MQ=255)
                                                                                            ccAGTATGTGGCGTCGGATGGATGAGCTGGCATAAATGACAGATGCCCTCTTCCCGCACGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTAAGGCCCGACTCTGTCGAACAGGCGCTGTAAATGAACCGCCGCTAc                                                          >  7:117415/1‑148 (MQ=255)
                                                                                              aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                       >  8:217566/1‑149 (MQ=255)
                                                                                              aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                       >  6:57407/1‑149 (MQ=255)
                                                                                                    tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGa                                                  <  7:66942/148‑1 (MQ=255)
                                                                                                       cGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGcc                                              <  2:371100/149‑1 (MQ=255)
                                                                                                            gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTg                                         >  3:48614/1‑149 (MQ=255)
                                                                                                             atggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTgg                                        <  4:26358/149‑1 (MQ=255)
                                                                                                               ggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGtt                                      >  2:333956/1‑149 (MQ=255)
                                                                                                                           aTAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGCAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                          <  8:168585/149‑1 (MQ=255)
                                                                                                                            tAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                          >  6:216925/1‑148 (MQ=255)
                                                                                                                                   cAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGCAGATGATTGAACAATTTTTCAGGCCCGCCTCTGTCGACCAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCACACTc                   >  5:381783/1‑148 (MQ=255)
                                                                                                                                       tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc              >  2:222336/1‑149 (MQ=255)
                                                                                                                                        gCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCggg                               >  7:166954/1‑131 (MQ=255)
                                                                                                                                        gCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCggg                               <  8:166954/131‑1 (MQ=255)
                                                                                                                                          ccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTCCCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa            >  7:11168/1‑148 (MQ=255)
                                                                                                                                          ccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac           >  3:315930/1‑149 (MQ=255)
                                                                                                                                          ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa            >  7:242463/1‑148 (MQ=255)
                                                                                                                                          ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACT‑AATGCTaca          >  5:249676/1‑149 (MQ=255)
                                                                                                                                              ttCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac         <  2:44805/147‑1 (MQ=255)
                                                                                                                                                cccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCaa      <  8:344886/148‑1 (MQ=255)
                                                                                                                                                   gCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCg  >  6:247479/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCG  >  NC_000913/3020363‑3020658

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: ‑