Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG26.1% 95.1 / 20.0 46G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/33);  new base G (12/0);  total (13/33)
Fisher's exact test for biased strand distribution p-value = 3.34e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATA  >  NC_000913/2326315‑2326580
                                                                                                                                       |                                                                                                                                   
aCTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                <  1:1010937/140‑1 (MQ=255)
     tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                           <  2:656129/140‑1 (MQ=255)
            ggggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                    <  1:851523/139‑1 (MQ=255)
            cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                    <  1:1103553/140‑1 (MQ=255)
            cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTTGAAAATATggggt                                                                                                                    >  1:1104476/1‑136 (MQ=255)
            cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGtt                                                                                                                    >  2:1103553/1‑140 (MQ=255)
            cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGtt                                                                                                                    >  2:1094570/1‑140 (MQ=255)
            cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGtt                                                                                                                    >  1:983994/1‑140 (MQ=255)
              ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                  <  1:1014417/140‑1 (MQ=255)
                 tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                                                               <  2:983994/140‑1 (MQ=255)
                 tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                                                               <  1:361305/140‑1 (MQ=255)
                 gggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                                                               <  1:1047955/139‑1 (MQ=255)
                    cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                              >  2:38926/1‑138 (MQ=255)
                    cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                              <  1:38926/138‑1 (MQ=255)
                         aCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATTAGTTTAGCCCCCcactt                                                                                                       >  2:947239/1‑140 (MQ=255)
                            gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                    <  2:702523/140‑1 (MQ=255)
                             tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGCGGTTAACCCCCTCATTACTc                                                                                                   >  2:660776/1‑140 (MQ=255)
                             tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTc                                                                                                   >  1:164950/1‑140 (MQ=255)
                             tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATTAGGTTTGCCCCCCACTTtctc                                                                                                   >  2:1047955/1‑140 (MQ=255)
                               tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGa                                                                                                 <  1:89492/140‑1 (MQ=255)
                                        aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGc                                                                                        <  1:472490/140‑1 (MQ=255)
                                           gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                     <  1:237835/140‑1 (MQ=255)
                                           agacagaaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                     <  1:219864/134‑1 (MQ=255)
                                            tCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                                                                    <  1:354007/140‑1 (MQ=255)
                                              aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc                                                                                  <  2:546212/140‑1 (MQ=255)
                                              aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc                                                                                  >  2:940190/1‑140 (MQ=255)
                                               aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGTAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctct                                                                                 <  2:791089/140‑1 (MQ=255)
                                                 aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTgg                                                                               <  1:358457/140‑1 (MQ=255)
                                                      aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTTGCCCCCCACTTACTCGATGCAAAAGCAAGCTCTGGTTAtt                                                                          >  2:354007/1‑139 (MQ=255)
                                                      aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTgggtttc                                                                          >  2:383446/1‑135 (MQ=255)
                                                        gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                        >  1:1077450/1‑140 (MQ=255)
                                                        gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                        >  1:416921/1‑140 (MQ=255)
                                                            ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTGCCCCCCACTTACTCGATGCAAAAGAACGCCCTGGGTTTCGTCTTg                                                                    >  2:167297/1‑140 (MQ=255)
                                                              cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCtttgg                                                                  >  2:16924/1‑137 (MQ=255)
                                                               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                 <  1:539546/140‑1 (MQ=255)
                                                               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                 <  2:416921/140‑1 (MQ=255)
                                                               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                 <  1:153743/140‑1 (MQ=255)
                                                               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTAGTTTTGCCCCCCCCTTTCTCGATGCAAAAACACGCCCTGGTTATCGTCTTgggg                                                                 >  2:539546/1‑138 (MQ=255)
                                                               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGGCTTgggg                                                                 >  1:391531/1‑138 (MQ=255)
                                                                    ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGa                                                            <  1:940190/140‑1 (MQ=255)
                                                                     tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                                                           <  1:16924/140‑1 (MQ=255)
                                                                     tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGCTTGGGAAATATCAGGTTATCCCCCTCCTTACTCTACGGCAAAGCACGCTCTGGCTATTGTCTTGGAGGCCGAc                                                           >  2:131893/1‑140 (MQ=255)
                                                                             tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGAGAGAACAATGGCTTTGCCCCCCACCCACTCAACGGAAAAGCAACCCCTGGTTGTACTCTTGGTGAGGGACAggagtt                                                    >  2:1114682/1‑136 (MQ=255)
                                                                                      ccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATc                                          <  2:486970/140‑1 (MQ=255)
                                                                                        gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                        <  1:167297/140‑1 (MQ=255)
                                                                                           cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                     <  1:131893/140‑1 (MQ=255)
                                                                                           cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                     <  1:239392/140‑1 (MQ=255)
                                                                                             ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                   <  1:1094570/140‑1 (MQ=255)
                                                                                             ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                   <  1:602868/140‑1 (MQ=255)
                                                                                                   tCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGGCCCCCCCCTACTCGAATCAAAATCACGCTCTGGGTTTCGTCTTCGAGACGGACAAGGTTATTACGGAATCCCTCGCGATGGcc                             >  2:762820/1‑140 (MQ=255)
                                                                                                            tCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgcg                    <  1:365619/140‑1 (MQ=255)
                                                                                                                   cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCa             <  1:91015/140‑1 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATg           <  2:932501/140‑1 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGCAAAAAAACGCTCTGGTTATCGTCTTTGGGACGGGCAGGGTTTTGACGTAAACCTGCGCGGTGGGCTGATGTGCGCCCCCCCTg           >  2:1120873/1‑140 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATg           >  1:789557/1‑140 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCCCCCCTg           >  1:304070/1‑140 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGACAGGTTTATGAAGTAAACCTGCGCGGTGGGCTGATGTGCGCCACCCCTg           >  2:91015/1‑140 (MQ=255)
                                                                                                                     gagCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGACCCCTACTTACTCGGTGCAAAAACAACCTCTGGTTATCGTCTTGGAGGCGGAGAGGTTTATTAAGTAATCCTGCGCGATGGCCTGATGTGAGCCACCCATg           >  2:494640/1‑140 (MQ=255)
                                                                                                                      agCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTAGCCCCCCACTTACTCG‑‑‑‑‑‑‑AGCACGCTCTCGGTATCGTCTTGGAGGCGGACAGGGTTTTGACGGAATCCTGCGCGGTGGCCTGATGTGCGCCACCCCTGGGTATCat   >  2:841522/1‑140 (MQ=255)
                                                                                                                              ggCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata  <  2:61100/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATA  >  NC_000913/2326315‑2326580

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: