Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC42.1% 96.3 / 20.1 57G384G (GGT→GGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/1);  new base C (0/24);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 3.32e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585
                                                                                                                                                      |                                                                                                                                                
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCacac                                                                                                                                                 >  2:978742/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCGATGGTCAATG‑CCACAcc                                                                                                                                                >  2:1086597/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc                                                                                                                                                >  2:1232977/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc                                                                                                                                                >  2:1129929/1‑151 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc                                                                                                                                                >  1:511604/1‑150 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCt                                                                                                                                              >  2:846074/1‑150 (MQ=255)
              aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACCCCCTGGCCCccgccg                                                                                                                                   <  1:862474/151‑1 (MQ=255)
                    aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACa                                                                                                                              >  2:655139/1‑150 (MQ=255)
                      aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGCTACaa                                                                                                                             >  1:989496/1‑149 (MQ=255)
                          aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                       >  2:796002/1‑151 (MQ=255)
                          aGGGAGCCAAGAAAATGTATGTCGACTTGCCAGCCATAATTATTCAACAGATTAATTACGGCAGATATTTATCGGTGAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACACTAGCCCCCGCCGATACAAAGCg                                                                                                                         >  1:512065/1‑149 (MQ=255)
                           gggccccaaaaaaaTTTATGTCCACTTCCCACCCATTATTATTCAAAAGATAAATGCCGGAAAATTTTTAGCGGAAAGGAGGAAATGCAGATCTTTTCACCACGTTAAATGTTAAACC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGtt                                                                                                                       <  1:889891/140‑1 (MQ=255)
                           gggAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGTTATTCGTCCCGTTAAATGGTCAATC‑CCACCCCCGGGCCCCCGCCGATACAAAGCGTTg                                                                                                                      <  2:1093220/151‑1 (MQ=255)
                             taGCCCAGAAAATGTATGTCCACCTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                                                                     >  1:578277/2‑150 (MQ=255)
                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                              >  1:369792/1‑151 (MQ=255)
                                     aaaaTGTATGTCGATTTCCCAGCAATTATTATTAAAAAGATTAATGCGGGCAGATTTTTATGGGAAGGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTAAAGG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                <  2:734372/147‑1 (MQ=255)
                                       aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                          >  2:1009642/1‑151 (MQ=255)
                                                 gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                >  1:156786/1‑151 (MQ=255)
                                                    ttGCCAGCCATCTTTATACAACAGATAAAGGACGGCAGATTTTTAGGGGTAGGGAGAAAAGGGAGAGCTTTTCGCCACGTTAAATGGTAAATG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                  <  1:194366/146‑1 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            >  2:392824/1‑151 (MQ=255)
                                                      gccagcctTCATTATCCAACAGATTAATGACGGCAGATATTTTTCGGTAGGGAGGAAAGGCAGAGTTTTTCGTCACGTT‑AATGGGTAATGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:576358/150‑1 (MQ=255)
                                                       ccagcaaTCATTATAAAACAAATTAATGGCGGGAGATATTTATCGGAAGGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAGTC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:334255/149‑1 (MQ=255)
                                                       ccacccaTAATTATCCAACAGATTAAGGACGGCATATTTTTATCGGTAGGGAAGAAAGCCAGCGCTATTCGTCCCGTTCAATGGTCATTC‑CCCCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                           <  2:25043/150‑1 (MQ=255)
                                                           ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                      >  2:131485/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  2:794048/1‑151 (MQ=255)
                                                                 ttattaAACAGATAAATGAGGGAAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAATGGTCAACG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                    <  2:703999/147‑1 (MQ=255)
                                                                   attAAACAGATTAATGACGGAAAATTTTTAGGGGTAAGGAGAAAATGCAGAGTTTTTCG‑CCCCTTTAAAGGGAAATGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                               <  2:1030670/150‑1 (MQ=255)
                                                                       aaaagttaaaggagggCAAATTTTTAGCGGAAAGGAGAAAAGGAAGAGTTTTTCTTCACTTTAAATGGTAAACC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                <  1:899407/131‑1 (MQ=255)
                                                                        aCAGATAAAGGACGGAAATTATTTACGTGTAGGGAGGAAATGCAGAGCTTTTCGTCCCGTTAAATGGTCAATG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                  <  1:1301399/142‑1 (MQ=255)
                                                                            aTTAATGCCGGAAAATTTTTATCGGAAGGGAGGAAAGGCAGTGCTTTTCGTCACGTTAAATGGTCATTC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                      <  2:149064/150‑1 (MQ=255)
                                                                             ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                    >  1:1236285/1‑151 (MQ=255)
                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                 >  1:223945/1‑150 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                               >  1:239441/1‑151 (MQ=255)
                                                                                    cGGCAGTTATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGAAAAGC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                             <  2:263278/151‑1 (MQ=255)
                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                          >  2:343215/1‑151 (MQ=255)
                                                                                        atatttttATAGGAAAGGAGAAAAGCCAGATTTTTTCGTCCCTTCCATTGGTCATGC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                            <  1:1159899/143‑1 (MQ=255)
                                                                                        aGATTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTCGATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                           <  2:1005767/149‑1 (MQ=255)
                                                                                           tatTTAGCGAAAGGGAGGAAATGAGGACTTTTTCTCCACGTTAAATGGTAAAGG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                           <  1:75787/146‑1 (MQ=255)
                                                                                              ttATCGGTAAGGAGAAAATGCAGAGCTTTTCGCCACTTTAAATGGTAAATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                    <  2:875736/150‑1 (MQ=255)
                                                                                               tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                   >  2:1052659/1‑150 (MQ=255)
                                                                                                      aaGGAGGAAATGGATATTTTTTCGTCACGTTCAATGTTCAATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                           <  1:1214760/151‑1 (MQ=255)
                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                            >  2:1147736/1‑150 (MQ=255)
                                                                                                        ggaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCTGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                         <  1:639200/151‑1 (MQ=255)
                                                                                                         gagaaaggcaagtgCTTTTCGTCACGTTAAAGGGTCAAGC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                        <  1:585544/138‑1 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                       >  2:569433/1‑151 (MQ=255)
                                                                                                                  cAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                >  2:1688/1‑150 (MQ=255)
                                                                                                                   agggCTTTTCGTCCCTTTAAATGGTCAACC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa                               <  2:476770/147‑1 (MQ=255)
                                                                                                                     agCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            >  2:1110838/1‑151 (MQ=255)
                                                                                                                        tATTCTTCCCTTCCAATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            <  1:795972/148‑1 (MQ=255)
                                                                                                                         aTTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                         >  2:202883/1‑150 (MQ=255)
                                                                                                                            cGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     >  2:618805/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                  >  2:597617/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     <  2:1033228/148‑1 (MQ=255)
                                                                                                                                aCGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                  >  1:743130/1‑150 (MQ=255)
                                                                                                                                  gTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc               >  1:216441/1‑151 (MQ=255)
                                                                                                                                        tGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac         >  2:405982/1‑151 (MQ=255)
                                                                                                                                        tGGTCAATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga          <  2:592331/150‑1 (MQ=255)
                                                                                                                                               tG‑CCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc  >  1:363170/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585

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