Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5330TG54.0% 8.0 / 30.7 38V131G (GTG→GGG) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base T (16/1);  new base G (20/0);  total (37/1)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGA  >  NC_000913/417295‑417761
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
ctcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAACCACCTCAAGCGCGAGTCGATGACCCGGGATATTCTAGTGGTGCTTTTGACCGCCAGAGGGTAAGAAGAAGATCGCGTGCGCGGCCTTGAAATCGGCGCGGATGACTATACCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGAGGTAATGCGCCGTATTTTGCCAATGGCGGTGGGAGACGt                                                                                                                                                                                                                                              >  2:1416603/1‑239 (MQ=255)
               ccTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGGGGGGATTGAGATGCAggg                                                                                                                                                                                                                               >  2:123066/1‑239 (MQ=255)
                    cggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGATTGAGATGCAGGGATTaa                                                                                                                                                                                                                          >  1:1162066/1‑239 (MQ=255)
                                                cTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGGTTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGAt                                                                                                                                                                                              >  1:1158456/1‑239 (MQ=255)
                                                          aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCGGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGTCGAGGCGGTTTCGAATCAGGGATTAAGTTTTCAACCGCCAACTCACCCGGGGAAggc‑‑gggaggg                                                                                                                                                                                  >  2:132763/1‑235 (MQ=255)
                                                                                  tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGa                                                                                                                                                                                                   >  1:320331/1‑200 (MQ=255)
                                                                                      gATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCAGCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGAAGAGCCGCTGGGGAGGGGGGCGGCAGAAt                                                                                                                                                  >  1:1687889/1‑239 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGGAAAAGCGCCGGGGATGGGGGCGgaagaaat                                                                                                                                                 >  2:1827689/1‑237 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACCTCTCACCGGGTGATGGC‑‑GGGGGGAGAGACGCCGGGGATGGGGGCgccccaaat                                                                                                                                                  >  2:578637/1‑231 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGAAGAGCCGCTGGGGATgg                                                                                                                                                               >  1:1276564/1‑226 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGCCTCGACCCGACATCTCa                                                                                                                                                                                                       >  1:1684516/1‑191 (MQ=255)
                                                                                           tGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGTCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCCCACCGGGGGCGGGC‑‑GGGCGGAAAGCCGCCGGGGATGGGGGCCACAAAATTTaaa                                                                                                                                             >  1:1710324/1‑239 (MQ=255)
                                                                                                  cAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGC‑‑GGGCGAAGAGCCGCGGGGGAGGGGGGCGACAGAATTTtaaat                                                                                                                                           >  1:1122204/1‑232 (MQ=255)
                                                                                                       gggAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGGCGGGATTTAGTCTCGACCCGACCTCTCACCGGGGGATGGCGGGGGGGAAGAGCGGGGGGAGGGGGGCCGACGGGTTTTAAAatgcgcc                                                                                                                                      >  2:72979/1‑232 (MQ=255)
                                                                                                                  gaTCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGGGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTGACCGAGTGATGGC‑‑GGGGGGAGGAGGGCCGGGGGGGGGGCCGGCCGGATTTTAACTGCTGGACtttttttttgagcc                                                                                                                      >  2:797592/1‑231 (MQ=255)
                                                                                                                  gaTCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAGGCCCGACCCGACATCTCACCGGGTGGTGGCGGGGGGAAGAGCCGCGGGGGAGGGGGCGGAAAGAATTTAAATTGCGGGCATTTTTTTTGAGGc                                                                                                                      >  2:206532/1‑239 (MQ=255)
                                                                                                                       cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGC‑‑GGGGGGGGAGCCGCCGGGCGAGGGGGCGCACGAAATTTAAATGCTGCCCCTTTTTTTGTCggctcccc                                                                                                                >  2:1248761/1‑232 (MQ=255)
                                                                                                                                ccTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGGGATGGGGCCGACAGAATTTAAACTGCTGCCCTTTTTTTTGACGCATCCTTAGCGCgtgg                                                                                                       >  1:568422/1‑238 (MQ=255)
                                                                                                                                          ggCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGGAAAGGCGCCGGGGGAGGGGGCGGCAGAATTTTAAATGGCGCACTTTTTTTTg                                                                                                                          >  1:1112215/1‑210 (MQ=255)
                                                                                                                                            cgcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGCGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGGCATCTCACCGAGTGATGGG‑‑GGGGGGAGAGGCGCCGGGGGGGGGGGCCGACGGAATTATAACTCGTGCCCTTTTTTTTGGGCACACCGTGAGCGTGTGTagcccgccg                                                                                              >  2:465046/1‑230 (MQ=255)
                                                                                                                                                           cACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAAAGGCCGCGGGGGAGGGGGCGGCCGAAATTTAAATGCTGCCCTTTTTTTTTACGCAACCTGAGTGCGGGTACAACCCCGAGCACCTGTTTAAacac                                                                            >  2:162452/1‑238 (MQ=255)
                                                                                                                                                             ccAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAGATGCAGGGATTAAGTCTCGAACCGACATCTCCCCGAGTGATGGC‑‑GGGGGGAGAGGCGCTGGAGGTGGGGGCGGCAGAATTTAAACTGCTGCACTTTTTTTTTTAGCATCCTGAGCGCGTGTACAGCCCCGAGCAGCTGTTATAACACg                                                                           >  2:1307064/1‑237 (MQ=255)
                                                                                                                                                             ccAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGa                                                                                                                                                                                                   >  1:1173491/1‑125 (MQ=255)
                                                                                                                                                              cAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCa                                                                                                                                                                                                       >  1:383898/1‑120 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGGAGAGCCGCTGga                                                                                                                                                                      >  1:1347913/1‑148 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGGAAAGGCGCCGGGGGAGGGGGCGGCAAAAATTTaaa                                                                                                                                             >  2:994096/1‑171 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTTAGATGCAGGGGTTAAGCCTCGACCCGTCATCTCACCGAGTGATGGC‑‑GGGCGTAGAAACGCCGGGGAATGGGGCCAGAGGAATTTCAATTGCGCACTTTTTTTTTACGCATCCTGTGCGCGTGGACAGCCGCGCGCAGCCGCTAAATCACCGCTgg                                                                      >  2:330576/1‑239 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGGAGAGGCGCTGGGGAAGGGGGCGGCAGAAATTAAACCGCTGGACTTTTTTTTGACGCATCCCGAGCGGGGGGTCAACCGCGCGCAGCTGTTTAAACCCCgctgg                                                                       >  2:1844527/1‑235 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAAATGCTGCACTTTTTTTTGACGCATCCTGAGCGCgtgga                                                                                                      >  1:178338/1‑206 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGAGATGGGGCCCACAGAATTTTAAATGCTGCCCCTTTTTATGACGCCTCCTGAGCGGGGGTACAGCCCCGGGGAGCTGTTAAACCACGGCTgg                                                                      >  1:806298/1‑239 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGCGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGGGATGCAGGGCATAAGTCTCGACCCGa                                                                                                                                                                                                              >  1:327905/1‑111 (MQ=255)
                                                                                                                                                                                   gcggtgggCGCGAATCAAAGCGGGAACGCGCCTTATTTCGCCAATGGCGGTGGAAGAGGTGATTGGAATGCAGGGATTAAGTCTCGACCCGACATCTCCCCGAGTGATGGC‑‑GGGCGAAGCGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTAtgtg                                                    <  2:884096/232‑1 (MQ=255)
                                                                                                                                                                                                         gTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                                    >  2:304420/1‑48 (MQ=255)
                                                                                                                                                                                                                    tATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGGAAAGGCGCTGGGGATgg                                                                                                                                                               >  2:567341/1‑101 (MQ=255)
                                                                                                                                                                                                                    tATTTCGCCAATGGCGGGGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGGGAAAAGGCGCTGGGGATTGGGGCCGCAGAATTTTAAATGCTGCCATTTTTTTTTACGCCACCTGGGCGCGGGTACAGCcgcgg                                                                                             >  2:79167/1‑164 (MQ=255)
                                                                                                                                                                                                                               tGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGG‑‑GGGGGGAGAGGCGCCGGGGATGGGGGCGAAAGAAATTAAAATGCTGCACTTTTTTTTGAAGCATACTGAGGGCGGGTACAGAAGCGAGCAGCTGTTTAAACAAGACTGGGGGACTAACGTGGATGGGGGAGAACGCACGGGAGGTGTGCACAATTGTCGCACGCGTAAAga        >  2:190928/1‑237 (MQ=255)
                                                                                                                                                                                                                                    gTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGGAGAGCCGCCGGGGATGGGGGCGGCAGAATTTAAAATGCTGCCCTTTTTTTTTACGCCTCCTGGGCGCGGGTAAAACCCCCGGGAGCTGGTTAACCACCGCTGGGGGACTAACGGGTATGTGGAAGAACGCACGGTCGGTGTCCCCATTTGTCGCCTGCCTAAAGCACTGGa  >  2:37237/1‑239 (MQ=255)
                                                                                                                                                                                                                                            ggTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGaa                                                 >  1:982380/1‑184 (MQ=255)
                                                                                                                                                                                                                                             gTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGGAGAGCCGCCGGAGAAGGGGCCGGCAGGAATTTAAATCGTGCACTTTTTTTTGAGGCATCCCGAGCGCGTGTACAGCCGCGAGCCGCTGTTAAACCACGCCTGGGGGACTAACGTGTTtggtgtaga                                               >  2:1722252/1‑178 (MQ=255)
                                                                                                                                                                                                                                              tGATTGAGATGCAGGGGTTAAGt                                                                                                                                                                                                                        <  2:327905/23‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
CTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑‑GGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGA  >  NC_000913/417295‑417761

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