Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,811 (T)7→6 intergenic (+48/+705) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8050T.93.2% 168.5 / 6.6 44intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base . (21/20);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGT  >  NC_002947/4740658‑4740945
                                                                                                                                                      |                                                                                                                                            
gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                              <  1:429292/148‑1 (MQ=255)
      gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTATACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                       <  8:80613/149‑3 (MQ=255)
      gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                       <  2:229876/149‑3 (MQ=255)
      gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                       >  5:57107/1‑147 (MQ=255)
       aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                      <  4:411306/149‑4 (MQ=255)
        aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                     <  3:288627/149‑5 (MQ=255)
        aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCAGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGACCACGGGCAGCCtattttattt                                                                                                                                     >  3:52564/1‑145 (MQ=255)
         agaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                                    >  7:112002/1‑144 (MQ=255)
                          cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                  >  5:281698/1‑149 (MQ=255)
                             aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAACTGGGGGGGGACGAAGGAGGCAAGAAGGCTGCCATCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                  >  1:140138/1‑149 (MQ=255)
                              cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                              <  8:112002/149‑1 (MQ=255)
                              cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                              <  6:317096/149‑1 (MQ=255)
                                 gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCCCTTGGCgggc                                                                                                           >  5:293997/1‑147 (MQ=255)
                                       ccGGCCGCGCCAGCTGTACTCCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                     <  6:323069/149‑1 (MQ=255)
                                          tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                   >  1:138777/1‑148 (MQ=255)
                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                <  3:76119/149‑1 (MQ=255)
                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                 <  6:114010/148‑1 (MQ=255)
                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                 >  2:390938/1‑148 (MQ=255)
                                              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                              <  5:47333/149‑1 (MQ=38)
                                               gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                             <  3:61422/149‑1 (MQ=255)
                                                ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAgg                                                                                            <  6:281698/149‑1 (MQ=38)
                                                  aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCTCTTAGCGGCCTTGCGCGCGTAGCGGcc                                                                                          >  4:415425/1‑149 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  1:390938/149‑1 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  5:243306/149‑1 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  6:166698/149‑1 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  2:138777/149‑1 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCCAAGAGGCCg                                                                                         >  2:117637/1‑149 (MQ=255)
                                                          acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                             >  4:320016/1‑138 (MQ=38)
                                                          acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                             <  3:320016/138‑1 (MQ=38)
                                                                        gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                          >  1:386346/1‑80 (MQ=38)
                                                                        gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                          <  2:386346/80‑1 (MQ=38)
                                                                                 cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGGCCGCGAAGAGGCCGGCACAGCCACCCCCAAACCGCAGATATtc                                                            >  4:92756/1‑148 (MQ=255)
                                                                                  gTTGCCCTGAAAGACCGCATGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                           <  1:117637/148‑1 (MQ=255)
                                                                                       ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCGCAGTTTTTATacgccc                                                     >  4:230782/1‑149 (MQ=255)
                                                                                        cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCAGAGTTTTTATacgcccc                                                    >  7:351385/1‑148 (MQ=255)
                                                                                         tcaaaGACCGCAAGTAAGCTTCAGTACCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCAGGCCTCATCGAGGACTTGAACGCGAATAGGCAGGCCCAGACACCAACAAACCTCAGTATATCAacgccca                                                    >  6:151225/3‑147 (MQ=255)
                                                                                                   cAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                         >  8:109988/1‑149 (MQ=255)
                                                                                                   cAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTACGCCCGCTCACGCAGcc                                         >  1:138168/1‑149 (MQ=255)
                                                                                                          gCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGGCCGCGGATAGGCCGGCACAGAACCCACAAAACCGAAGGTTTTCTCCGCCCGCGCACCCAGCCCGTTCAg                                  >  7:363249/1‑149 (MQ=255)
                                                                                                           cTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                 <  7:109988/149‑1 (MQ=255)
                                                                                                           cTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                 <  7:336897/149‑1 (MQ=255)
                                                                                                            tGCAGTAGCCGA‑‑‑AAGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                >  2:119203/1‑149 (MQ=255)
                                                                                                                    ccGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCGCAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTTAAc                        >  3:37898/1‑149 (MQ=255)
                                                                                                                               aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTGTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                >  1:343986/1‑149 (MQ=255)
                                                                                                                                 gatggCTGCCGTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca              <  2:138168/146‑1 (MQ=255)
                                                                                                                                  aagccggCCCTCGGGCAGCC‑TTTTTTATTTGGCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac             <  4:37898/143‑1 (MQ=255)
                                                                                                                                  aagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac             >  5:106519/1‑149 (MQ=255)
                                                                                                                                    ggcggcCCTCGGGCAGCC‑TTTTTAATTGGCCCTGTGCCGCCCTCTGCGCGCCGTGCCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca           <  3:415425/149‑1 (MQ=255)
                                                                                                                                    ggCTGCCGTCGGCCAGCTTTTTTTTTTTGCCCTTG‑GCCGCCCTCTTGCCCGCCTTCCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca           <  4:219462/149‑1 (MQ=255)
                                                                                                                                        gCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAAGaa       >  1:329460/1‑149 (MQ=255)
                                                                                                                                         cccTCGGGCAGCTTTTTTTTTTTGTCCTGG‑GCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc      <  3:145983/149‑1 (MQ=255)
                                                                                                                                              gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                            <  7:246489/122‑1 (MQ=39)
                                                                                                                                              gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                            >  8:246489/1‑122 (MQ=39)
                                                                                                                                              gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt  <  2:140138/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGT  >  NC_002947/4740658‑4740945

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