Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,1730TG23.1% 109.6 / 10.0 53G31G (GGT→GGGymdFconserved protein YmdF
Reads supporting (aligned to +/- strand):  ref base T (18/22);  new base G (12/0);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 5.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCG  >  NC_000913/1068034‑1068312
                                                                                                                                            |                                                                                                                                           
cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGTCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                                             <  2:1031648/140‑1 (MQ=255)
           tGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTGCGGGAATTTTGTAGAAGAGCGCGGAAAAGCAACCATAGCAAGTAAAAAAAGAGGGCAACACAGGCGGGGAAAAtttaacaa                                                                                                                                   >  2:513928/1‑135 (MQ=255)
                   gatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGc                                                                                                                          >  2:844959/1‑140 (MQ=255)
                   gatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGAGGACAGCACAGCGGGGGGAAATTCAAGAATGACCCGc                                                                                                                          >  2:1132655/1‑140 (MQ=255)
                     tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAg                                                                                                                        >  2:953121/1‑140 (MQ=255)
                     tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGATCCCGCCAATTTTTTGAAGACCACGCAAGAGCATCCGAAGGCGGGAAAAAAGGAGGACAGGACAAGGGGGGGAATTTTAAGGATGACCCGccg                                                                                                                        >  2:990821/1‑138 (MQ=255)
                       tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcg                                                                                                                      <  1:228750/140‑1 (MQ=255)
                       tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcg                                                                                                                      <  2:496705/140‑1 (MQ=255)
                       tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTACGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcg                                                                                                                      <  2:1157834/140‑1 (MQ=255)
                        cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgc                                                                                                                     <  2:136316/140‑1 (MQ=255)
                        cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgc                                                                                                                     <  2:730867/140‑1 (MQ=255)
                          ataCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAACTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAgcgcgc                                                                                                                   >  2:190088/1‑140 (MQ=255)
                            aCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCAt                                                                                                                 >  2:607126/1‑140 (MQ=255)
                              ccAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCt                                                                                                               >  2:314019/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                            <  2:885755/140‑1 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                            <  1:587689/140‑1 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                            >  1:308228/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                            >  1:467134/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                            >  1:687316/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGGTCCGGCAATTTTGCAGAAGACCGCGAAAGAACATCAGAAGCACGGAAAAAAGGTGGACAGCACAGCGGGGGCAATTTTAAGAAAAACCCCCAGCGCGCCTCTGaa                                                                                                            >  2:671929/1‑140 (MQ=255)
                                   cGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAAGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGGGCATCTGAAGc                                                                                                          >  2:408863/1‑140 (MQ=255)
                                    ggAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCa                                                                                                         <  1:65400/140‑1 (MQ=255)
                                         tCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaa                                                                                                    <  2:1052862/140‑1 (MQ=255)
                                            gTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGAGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCCGCGCGCATCTGAAGCAGGGaaaaa                                                                                                 >  1:1041821/1‑140 (MQ=255)
                                                  gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGGAGGGAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAggtgg                                                                                           >  1:632303/1‑140 (MQ=255)
                                                  gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGAACGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAggggg                                                                                           >  1:1151274/1‑140 (MQ=255)
                                                         aTCGAGGCGGTTCCGGCAATTTTGCAGAAAACCGCGAAAAAGCATCAGAAGCCGGGAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAAAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGGGGGAaaaga                                                                                    >  2:292113/1‑139 (MQ=255)
                                                            gAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGt                                                                                 <  2:864742/140‑1 (MQ=255)
                                                             aGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACAGCGGGGGTAATTTCAAGAAAGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGTGGGAAGAGCAGTc                                                                                >  2:955329/1‑140 (MQ=255)
                                                                cGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACg                                                                             >  1:427209/1‑140 (MQ=255)
                                                                   ttCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCa                                                                          <  2:308228/140‑1 (MQ=255)
                                                                    tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTTAAGGAGGGAAAAAAGGGGGGAAGAGCAGGCACGGCaa                                                                         >  2:158082/1‑140 (MQ=255)
                                                                            ttttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAACACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACa                                                                 >  2:577642/1‑140 (MQ=255)
                                                                                 cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                                            <  2:170724/140‑1 (MQ=255)
                                                                                 cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                                            <  1:1132655/140‑1 (MQ=255)
                                                                                 cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                                            <  1:314019/140‑1 (MQ=255)
                                                                                   gaagaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGcc                                                          <  1:440298/140‑1 (MQ=255)
                                                                                      gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                                       <  1:844959/140‑1 (MQ=255)
                                                                                         cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa                                                    >  2:647062/1‑140 (MQ=255)
                                                                                         cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa                                                    >  2:513943/1‑140 (MQ=255)
                                                                                                gagCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATg                                             >  2:517558/1‑140 (MQ=255)
                                                                                                  gCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGAc                                           <  1:190088/140‑1 (MQ=255)
                                                                                                    aTCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa                                         <  2:687316/140‑1 (MQ=255)
                                                                                                         aaGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa                                    <  1:607126/140‑1 (MQ=255)
                                                                                                         aaGCAGGTAAAAAAGGTGGACAACACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa                                    <  1:577642/140‑1 (MQ=255)
                                                                                                           gCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt                                  >  2:205955/1‑140 (MQ=255)
                                                                                                           gCAGGTAAAAAAGGTGGACAGCACAGCGGGGGAAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGCGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt                                  >  2:436399/1‑140 (MQ=255)
                                                                                                             aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCAACAACTAGCCGGGCTAATCAATGACGAATGCAtttt                                >  1:1045894/1‑140 (MQ=255)
                                                                                                              ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGCCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                               <  1:990821/140‑1 (MQ=255)
                                                                                                              ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                               >  2:849398/1‑140 (MQ=255)
                                                                                                               gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg                              >  1:17228/1‑140 (MQ=255)
                                                                                                               gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACACGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg                              >  1:736834/1‑140 (MQ=255)
                                                                                                                       ggTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc                      >  1:269148/1‑140 (MQ=255)
                                                                                                                         tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc                    >  2:648315/1‑140 (MQ=255)
                                                                                                                                cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaa             >  2:607229/1‑140 (MQ=255)
                                                                                                                                cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaa             >  1:131009/1‑140 (MQ=255)
                                                                                                                                       gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGAAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTAAAAAAGCCAt       <  2:269148/140‑1 (MQ=255)
                                                                                                                                           gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc   <  1:75964/140‑1 (MQ=255)
                                                                                                                                           gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc   >  2:219479/1‑140 (MQ=255)
                                                                                                                                            tAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAccg  <  1:648315/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCG  >  NC_000913/1068034‑1068312

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