Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5330TG47.3% 39.0 / 17.8 56V131G (GTG→GGG) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base T (20/9);  new base G (26/0);  total (47/9)
Fisher's exact test for biased strand distribution p-value = 2.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GCCGC‑TGGAGATGGGGCCGACAG‑AATTTAAAC‑TGCTGCA‑CTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAG‑CTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACAT  >  NC_000913/417295‑417737
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                   
ctcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGt                                                                                                                                                                                                                     <  1:1784512/239‑1 (MQ=255)
        gATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTgaga                                                                                                                                                                                                             >  1:1823323/1‑239 (MQ=255)
               ccTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAAGCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGTGAATGAGATgacggg                                                                                                                                                                                                      >  2:1010686/1‑234 (MQ=255)
                  ggcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGAtt                                                                                                                                                                                                   >  1:964982/1‑239 (MQ=255)
                        tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGGGGGGATTTGGAGTCAGGGATTAAGtc                                                                                                                                                                                              >  2:2741712/1‑238 (MQ=255)
                        tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGCGGTGGTTGGGATGCAGGGGTTAAgac                                                                                                                                                                                              >  2:1740784/1‑236 (MQ=255)
                        tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTAGCCAATGGCGGGGGAAGAGGTGATTGAGATGCGGGGCGTAAGtc                                                                                                                                                                                              >  2:2726757/1‑238 (MQ=255)
                            gTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGATTTGGATGGCGGGGTTAAGTCTCGa                                                                                                                                                                                          >  1:1535513/1‑238 (MQ=255)
                                    ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAGTCTCGACCCGACa                                                                                                                                                                                   >  1:2729522/1‑237 (MQ=255)
                                      catcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCACTGGCGGTGGAAGAGGGGATTGAGGTGCAGGGATTAAGTCTCGACCCGaactccc                                                                                                                                                                               >  2:2642130/1‑233 (MQ=255)
                                       atcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAggg                                                                                                                                                                                                      >  1:156671/1‑215 (MQ=255)
                                       atcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACAtctc                                                                                                                                                                               >  1:362639/1‑238 (MQ=255)
                                       atcaAACACCGCAAGCGGGAGGCGATGAACCGATAGATGCAAGTGGTGAGGGTGGCCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGTCGCGGAGGACTATATCTCCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGGAATGCTCCGTATTTCGCCAATGTCGGTGGAAGAGGTTATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCa                                                                                                                                                                              <  2:2671955/239‑1 (MQ=255)
                                         caAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGGCAtcccccc                                                                                                                                                                            >  1:2327524/1‑234 (MQ=255)
                                              acCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAGGAGGGGATTGAGAGGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTg                                                                                                                                                                       >  2:2505670/1‑239 (MQ=255)
                                                cTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGGGTGAt                                                                                                                                                                     >  1:866257/1‑239 (MQ=255)
                                                 tCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCCGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCTGTATTTCGCCAATGGCGGCGGAAGAGGCGGTTGAGTGGCAGGGGTTATGTACCGACCCGACACCTCACGGAGTGTTg                                                                                                                                                                    >  2:29574/1‑239 (MQ=255)
                                                  cAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGGGGTGATTGAGATGACGGGGGTAAGTCTCGACCCGACATCTCACCGAGTGagtg                                                                                                                                                                   >  2:537007/1‑236 (MQ=255)
                                                  cAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATg                                                                                                                                                                    >  1:1237251/1‑238 (MQ=255)
                                                          aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTTAGCTGCAGGGATTTAGTCTTGACCCGCCATCTTACCGTGGGGTggc‑gggtgga                                                                                                                                                          >  2:631716/1‑235 (MQ=255)
                                                                                  tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAAAGGTTATTGAGATGCAGGGGGTAAGTTTCGAGCCGGCCTCTCACCGGGGTGTGGC‑GGGGGGAGG‑GACGG‑TGGCGATgggtgcgga                                                                                                                                 >  2:294786/1‑232 (MQ=255)
                                                                                  tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGACTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GGCCC‑CGGGGCTGGGGCCGGCa                                                                                                                                >  1:1889385/1‑239 (MQ=255)
                                                                                    gtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGAGAA‑GCCCC‑CGGGGATGGGGCCgcccga                                                                                                                              >  1:2655472/1‑234 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAggg                                                                                                                                                                                                      >  1:137516/1‑167 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGAAGA‑GCCGC‑TGGGGGTGGGGCCGGCAG‑AAtt                                                                                                                          >  1:39852/1‑239 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGC‑GGGGGGAGA‑GCCGC‑CGGGGATgggg                                                                                                                                      >  1:669710/1‑228 (MQ=255)
                                                                                       aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGGGGGGATTGGGATGCAGGGATTAAGTCTCGATCCGCAATCTCAACGAGGGGAGGC‑GGGGGGAGA‑GGCGC‑GGGGGGTGGGGGCGAccacaat                                                                                                                           >  2:1856123/1‑233 (MQ=255)
                                                                                            gACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGGGGGGATTGAGATGCAGGGATTCAGTCTCGCCCCGACATTTCACCGCGTGATTGC‑GGGGGGAGA‑GGCGCGTGGCG‑TGGGGCCGGCAG‑AAtattaa                                                                                                                      >  2:431660/1‑233 (MQ=255)
                                                                                             aCCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGGAGA‑GCCGC‑CTGGGGTGGGGGCCACGG‑AATTTTAAcc                                                                                                                    >  1:1386670/1‑238 (MQ=255)
                                                                                               cGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGATTGAGATGCAGGGGTTAAGGCTCGACCCGACATCTCACCGAGTGGTGGC‑GGGGGAGGAGGCCGG‑TGG‑GAATGGGCCGGCCG‑GAATTAAAc‑t                                                                                                                   >  2:2190405/1‑237 (MQ=255)
                                                                                                  cAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGGGGA‑GGCGC‑GGGGGGTGGGGcgcgc                                                                                                                                 >  2:467285/1‑218 (MQ=255)
                                                                                                               gaagaTCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGGGGTGATTGGGATGCAGGGATTAAGTCTCGACCCGACATCTCCCCGAGGGATGGCGGGGGGAGGA‑GCCGG‑CGGGGGGGGGGCCGACGC‑AATTTAAACGTGGTCCA‑Cattttttgt                                                                                                  >  2:2226584/1‑236 (MQ=255)
                                                                                                                       cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCCGGGATTAAGTCTCGTCCCGCCa                                                                                                                                                                                   >  1:2568575/1‑154 (MQ=255)
                                                                                                                       cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGt                                                                                                                                                                                               >  1:2145964/1‑142 (MQ=255)
                                                                                                                          gcgcGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGTGCTGGTGGCGCG‑ATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTTTCGACCCGACATCTCAccca                                                                                                                                                                          >  1:2124725/1‑157 (MQ=255)
                                                                                                                                                        tatCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGAt                                                                                                                                                                                                    >  1:59/1‑104 (MQ=255)
                                                                                                                                                          tCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGGAGA‑GCCGC‑CGGGGATGGGGCCGGCAG‑AATTTTAAC‑TGCTGCA‑CTTTTTTATTTAGGATCCTGGGCGCGGGGACAg                                                                           >  2:2317408/1‑217 (MQ=255)
                                                                                                                                                           cACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGAt                                                                                                                                                                                                    >  1:1172234/1‑101 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGC‑GGGGGGAGA‑GGCGC‑CGGGGGGGGGGGCCCCAGAAATTTAAAC‑‑GGTGCC‑CCTTTTTTTTAAGCATCCCTTGGGCGTGGGTAGCCCCGAGGCG‑CTGTTTAAACACCGCTgg                                              >  2:1266766/1‑239 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACAt                                                                                                                                                                                  >  1:2331096/1‑114 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGGTTGGGATGCAGGGATTAAGGCTCGACCCGACCTCTCACCGAGTGATGGC‑GGGGGAAAA‑GCCGC‑GGGGGATGGGGCCGGCCG‑GATTTTAAC‑TTCTGGC‑CTTTTTTTTTACGCAACCCGAGAGCGTGTACAGCCGCGAGCGA‑GTGTTTAAcccttctcct                                              >  2:1626369/1‑231 (MQ=255)
                                                                                                                                                                aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGTGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGGAAA‑GGCGC‑CGGGGGTGGGGGCGACCG‑CAATTAAAA‑TGCTGCC‑CTTTTTTTTGATGCATACTGAGCGCGGGTTACGACGCGAGCAG‑CGGTTTAACCACGCCTgg                                              >  2:916081/1‑239 (MQ=255)
                                                                                                                                                                     ttttCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGGGGAGG‑GCCGC‑CGGGGGTGGGGGCGGCAG‑AAATTTAAC‑TTCTGCACCTTTTTTTTTAC‑CACTCCGAGCGCGGGTATAGCCGCCGGCAGCCGGTTAAAcccaccttgggc                                            >  2:715744/1‑227 (MQ=255)
                                                                                                                                                                                   gCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGGGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GCCGC‑TGGAGATGGGGCCGACAG‑AATTTAAAC‑TGCTGCA‑CTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAG‑CTGTTAAACCACGTCTGGGGAACTAACGTGTAtgtg                            <  1:1347201/238‑1 (MQ=255)
                                                                                                                                                                                            gcgAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCa                                                                                                                                                                              >  1:317838/1‑90 (MQ=255)
                                                                                                                                                                                            gcgAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTGGACCCGACATCTCACCGAGTGATGGC‑GGGGGAAGA‑GCCGC‑TGgaga                                                                                                                                           >  1:1660714/1‑122 (MQ=255)
                                                                                                                                                                                                             tGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGt                                                                                                                                                                                               >  1:1365911/1‑56 (MQ=255)
                                                                                                                                                                                                                         cGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGtc                                                                                                                                                                                              >  2:522796/1‑45 (MQ=255)
                                                                                                                                                                                                                          gCCAATGGCGGTGGAAGAGGTGa                                                                                                                                                                                                                   <  2:2659992/23‑1 (MQ=255)
                                                                                                                                                                                                                             aaTGGCGGTGGAAGAGGTGATTg                                                                                                                                                                                                                <  2:156671/23‑1 (MQ=255)
                                                                                                                                                                                                                             aaTGGCGGTGGAAGAGGTGATTg                                                                                                                                                                                                                <  2:137516/23‑1 (MQ=255)
                                                                                                                                                                                                                              aTGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GCCGC‑TGGAGATGGGGCCGACAG‑AATTTAAAC‑TGCTGCA‑CTTTTTTTTGACGCATCCTGAgcgc                                                                                   >  1:2300824/1‑141 (MQ=255)
                                                                                                                                                                                                                               tGGCGGTGGAAGAGGTGATTgag                                                                                                                                                                                                              <  2:1172234/23‑1 (MQ=255)
                                                                                                                                                                                                                               tGGCGGTGGAAGAGGTGATTgag                                                                                                                                                                                                              <  2:59/23‑1 (MQ=255)
                                                                                                                                                                                                                                    gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                         <  2:2145964/23‑1 (MQ=255)
                                                                                                                                                                                                                                    gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                         <  2:1365911/23‑1 (MQ=255)
                                                                                                                                                                                                                                             gTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GCCGC‑TGGGGATGGGGCCGCCAG‑AATTTAAAC‑TGCTGCA‑CTTTTTTATGACGCATCCTGAGCGCGGGTACAGCCCCGAGCAG‑CTGTTTAACCACGTCTGGGGGACTAACGTGTATGTGGAAAACCGCACGGTCGATGTCCACAt  >  2:2058698/1‑206 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                   
CTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGC‑GGGCGAAGA‑GCCGC‑TGGAGATGGGGCCGACAG‑AATTTAAAC‑TGCTGCA‑CTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAG‑CTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACAT  >  NC_000913/417295‑417737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑