Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,473,8840TG24.6% 87.3 / 10.3 61L94R (CTG→CGG) yfdPCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base G (15/0);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGT  >  NC_000913/2473742‑2474020
                                                                                                                                              |                                                                                                                                        
ctcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCt                                                                                                                                          <  1:15115/143‑1 (MQ=255)
             tCGAACCATTTGATGAGAGCCGGCTGAAACTCAAACACAAACTCTCATAGGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAg                                                                                                                             <  1:158627/143‑1 (MQ=255)
             tCGAACCATTTGAGGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAg                                                                                                                             <  1:783215/143‑1 (MQ=255)
               gAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGa                                                                                                                           >  1:975694/1‑143 (MQ=255)
               gAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGAATAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTTACTATTCTTCAGGa                                                                                                                           >  2:900170/1‑143 (MQ=255)
                  tgatgtgaagaGAAGCGTCTGAAGATCGAACACAAACGCGCAGATGGTCGCCCGACTAACCGTGGGACGATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGa                                                                                                                        <  1:1053420/138‑1 (MQ=255)
                       tgatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAgg                                                                                                                   <  1:992112/143‑1 (MQ=255)
                        gatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGACGAAGACACCACCGAAACGCCGATGTATGGCAATCGCGGTGGTCGCCTGACCATTCTGCAGGAAGACCAggt                                                                                                                  >  2:953908/1‑142 (MQ=255)
                        gatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCTGACTATTATGCAAGAAGACCAggg                                                                                                                  >  2:468644/1‑143 (MQ=255)
                         atgaGAACCGTCGGAAAGTCAACCACAAACTTCGATATGTGCTCCCGTGAAACCGTGGGAAAATTTCCGACGCAGACACCACCGAACCGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGa                                                                                                                 <  2:753500/143‑1 (MQ=255)
                           gagaACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTAAAAGAAGACACGGGACa                                                                                                               >  2:738304/1‑143 (MQ=255)
                                     tGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTGCAAAAAAACCAGGCACAATTACTGCct                                                                                                     >  2:545480/1‑143 (MQ=255)
                                          aTCAAACACAAACTCTCGTACGTTCGCCCGACTGACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGc                                                                                                <  2:910736/143‑1 (MQ=255)
                                          aTCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGTACTATTCTGCAGTAAGACCAGCGACAATTACTGACTCTTGc                                                                                                >  2:444022/1‑143 (MQ=255)
                                            cAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGACGACCAGGGACAATTACTGACTCTTGCCg                                                                                              >  1:1099465/1‑143 (MQ=255)
                                            cAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCg                                                                                              <  1:130999/143‑1 (MQ=255)
                                              aaCACAACCTCTCATTTGTCCGCCCGACTAACCGTGGGACAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTTTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGt                                                                                            <  1:607635/143‑1 (MQ=255)
                                              aaCACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTATACAGGAAGACAAGCGACAATTCCTGACTCTTGCCGGt                                                                                            >  2:838986/1‑143 (MQ=255)
                                                cacaAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGa                                                                                          <  1:949980/143‑1 (MQ=255)
                                                   aaaCTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGCGcc                                                                                       >  1:172881/1‑143 (MQ=255)
                                                         tcATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCAGCAAGAAGACCAGCGACACTTACTGACTCTTGCCGGGGAGCCTGCCgg                                                                                 >  2:727072/1‑143 (MQ=255)
                                                         tcATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCACGGTGGTCGCCTGACTATTCTGCACGAAGACCAGGGACAATTACTGACTCTTGCCGGCGAGCATGACgg                                                                                 >  1:451332/1‑143 (MQ=255)
                                                           ttttGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaa                                                                               <  1:1157497/140‑1 (MQ=255)
                                                           cgatGTGCGGCAGACGAGCCGTGGGTAAATTTCGGAAGCAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaa                                                                               <  2:795124/141‑1 (MQ=255)
                                                                     ccGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAAGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGACAACCTCCgcgc                                                                     >  2:43309/1‑143 (MQ=255)
                                                                               gTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGACAACTCCGCGCCGCAGGTcgt                                                           >  1:1137006/1‑141 (MQ=255)
                                                                                   gAAAATTTCCGAAGAATACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAat                                                       <  1:549363/143‑1 (MQ=255)
                                                                                   gAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGCAAGAAAACCAGCGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCCCAGCCGGTCATTact                                                       >  2:863151/1‑141 (MQ=255)
                                                                                    aaaaTTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAata                                                      >  1:354937/1‑143 (MQ=255)
                                                                                      aaTTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATc                                                    <  1:1146312/143‑1 (MQ=255)
                                                                                       aTTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCAGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGCAAACTCCGCGCAGCAGGTCATTAATATcc                                                   >  2:548566/1‑142 (MQ=255)
                                                                                        tttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTGCAGGAAGACCAGGGGCAATTACTGACTCTTGCCGGTGAGCCTGACGGAACACTCCGCGCAGCAGGTCAGTAATATCAt                                                  >  2:846379/1‑143 (MQ=255)
                                                                                        tttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCAACAAAAAGACCAGGGACAATTACTGACTCATGCCGCTGAGCCTGCCGGAAAACTCCGCGCAGCCGGCCCTTAATATCac                                                  >  2:957021/1‑142 (MQ=255)
                                                                                         ttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAtt                                                 <  2:1051544/143‑1 (MQ=255)
                                                                                         tcccGAACAAGACACCACCGAATCGCCGATGTATGTCAATCGCGGTGGCCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACCCCGCGCAGCAGGTCATTAATATCAtt                                                 <  1:805212/141‑1 (MQ=255)
                                                                                             gaagaaGACACCACCGAAACGCCGATGTATGTTAATCGCGGTGGTCGCCTGACCATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCtta                                             >  2:893825/1‑143 (MQ=255)
                                                                                                gcggcaaccaccGGTGCGCCGATGTGTGCCAATCGCGGTGGCCGCCTGACTAGTCTGCAGGAAGACCAGGGACAATCACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaatt                                          <  2:663938/137‑1 (MQ=255)
                                                                                                gaagaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaatt                                          >  2:910395/1‑143 (MQ=255)
                                                                                                gaagaCACCACCGAAACGCCGATGTATGTCAAGCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaatt                                          <  1:953908/143‑1 (MQ=255)
                                                                                                   gaCACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGACCATTCTACAGAAAAACCAGAGACAATTACTGACTCTTGCCGGTGAACCAGACGGAAAACTCAGCGCAGCAGGTCATTAATATCATTCTTAATTAAc                                       >  2:598485/1‑143 (MQ=255)
                                                                                                   gaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAATTTCGCGCAGCAGGTCATTAATATCATTCTTAATTAAc                                       >  1:98329/1‑143 (MQ=255)
                                                                                                        caccGAAACGCCGATGTATGTCAATCTCGGTGATCGCCGTACCATTCAACACAAAGACCAGGGACAATTACGGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCAGTAATATCATGCTTAATTAACTaaat                                  >  2:615475/1‑141 (MQ=255)
                                                                                                         accGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTATGCAAGAAGACCAGGGACAATTACTGACTCTTGCAGGTGAGCCTGGCGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATAAAATAatta                                 >  2:56684/1‑143 (MQ=255)
                                                                                                         accAAAACGCCGACGTCTGTCAATGGCGATAGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAACACTCCGCGCAGCAGGTCATTACTATCATTCTTCCTTAAATAatta                                 >  1:548430/1‑143 (MQ=255)
                                                                                                          ccGAAGCGCCGATGTATGTCAAGCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAattat                                <  1:304274/143‑1 (MQ=255)
                                                                                                                  ccGATGTATGTCAATCTCGGTGGTCGCCGTATCATTAAACAAAAAAACCAGCGACAATTACTGACTCTTGCCGGTGAGCCTGACGAAAAACTCCGCGCACAAGGACATTAACATCATCCACAATTAACTAAttatttatacca                        >  2:161962/1‑139 (MQ=255)
                                                                                                                   cGATGTATGTCAATCGCGGTGGTCGCCTGACGATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtcat                       <  1:530172/143‑1 (MQ=255)
                                                                                                                    gATGTATGTCAATCTCGGTGGTCGCCGGACTATTCAACAAGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTACGAGCAGCAGGGCATTCACAACATACTTAATTAACTAATTAGTTAGCgcatc                      >  2:217876/1‑143 (MQ=255)
                                                                                                                       gtatgtCAATCGCGGTGGTCGCCTGATTATTAAAAAGAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTAATCAca                   >  2:536872/1‑142 (MQ=255)
                                                                                                                        tatgtCAATCTCGGTGGTCGCCGGATCATTAAAAAAAAAAACCAGGAACAATTACTGACTCTTGCCGGTGAACCTGACGGAAAACTCCGCGCAGCAAGACATTAATATCAATCTTAAATAACTAATTAATTatcacatcacaa                  >  2:345216/1‑141 (MQ=255)
                                                                                                                         atgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGAGGAGCTTGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGa                 >  1:1166206/1‑143 (MQ=255)
                                                                                                                          tgtCAATCGCGGTGGTCGCCGTACCATTCTACAAAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTACCTCATCACTGaa                >  2:1003800/1‑143 (MQ=255)
                                                                                                                          cgtCACTCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTGCTGACTCTTGCCGGTGAGCCTGACGGAAATCTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGaa                <  2:578565/142‑1 (MQ=255)
                                                                                                                            tCAATCGCGGTGGTCGCCTCACTATTAAGAAGGAAGACCTGGGACAATTACTGACTCTTACAGGTAGGCCTGACGGAAAACTCCGCGCAGCAGGTGAATAATATCATTCTTAAGTAAACAATTATCTATATcaccaccaaatc              >  2:528012/1‑137 (MQ=255)
                                                                                                                            tCAATCGCGGTGGTCGCCGGATTATTAAGAAAAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTAAGCGCAGCAGGTCATTAATATCATTCTTAAATAACTAATTAATTAACTCATCACTGAata              >  2:1109125/1‑143 (MQ=255)
                                                                                                                             cAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAatat             <  1:1034446/143‑1 (MQ=255)
                                                                                                                               aTCGCGGTGAACGCCTGACTATTCTGCAGGAAGAACAGGGACACTTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCt           <  2:745730/143‑1 (MQ=255)
                                                                                                                                  gcgGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTaa        <  2:12539/143‑1 (MQ=255)
                                                                                                                                   cgGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTCGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAat       <  1:1066984/143‑1 (MQ=255)
                                                                                                                                        gtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCCTCGCGGAATATCTTAATATAGt  >  1:280524/1‑143 (MQ=255)
                                                                                                                                        gtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGt  >  1:81163/1‑143 (MQ=255)
                                                                                                                                        gtCGCCTGACTATTCGGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGt  <  1:43309/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
CTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGT  >  NC_000913/2473742‑2474020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: