Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,378,9680CT45.8% 32.7 / 40.4 73E98K (GAA→AAA) menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase; SEPHCHC synthase
Reads supporting (aligned to +/- strand):  ref base C (14/25);  new base T (33/0);  total (48/25)
Fisher's exact test for biased strand distribution p-value = 1.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAAT  >  NC_000913/2378819‑2379118
                                                                                                                                                         |                                                                                                                                                        
ccAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTtc                                                                                                                                                          <  1:1468695/150‑1 (MQ=255)
 cAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTtct‑c                                                                                                                                                       <  1:1186664/151‑1 (MQ=255)
     gTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGt                                                                                                                                                   >  1:754595/1‑151 (MQ=255)
          gggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑AATC‑A‑GTTTTTTT‑CCgggttac                                                                                                                                               >  2:111152/1‑145 (MQ=255)
          gggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAAGG‑GAAC‑A‑GTTTTTTT‑CCCGTTTAc                                                                                                                                               >  2:1183958/1‑150 (MQ=255)
            gATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑GATC‑A‑GTTTTTTT‑CCCGTTAAacc                                                                                                                                             >  1:410223/1‑150 (MQ=255)
            gATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAGAAATC‑G‑GTTTTT‑T‑CCGGTTAACCCg                                                                                                                                            >  2:372482/1‑151 (MQ=255)
            gATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAAAG‑AATT‑A‑GTTTTTTT‑CCGGGTAAACCg                                                                                                                                            >  1:582259/1‑151 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑GAAA‑C‑GTTTTTTT‑CCGGGTAACCCgg                                                                                                                                           >  1:1433029/1‑149 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAACCCGGc                                                                                                                                          >  1:571386/1‑150 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTATGCG‑AATA‑ACGGTTTTTT‑CCGTTTAAccgg                                                                                                                                            >  1:42328/1‑147 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATT‑A‑GTTTTTTT‑CCGGTTTACCCgg                                                                                                                                           >  2:1013538/1‑149 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGGTAAacccg                                                                                                                                           >  2:393017/1‑147 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAAAG‑AATC‑A‑GTTTTTTT‑CCGGTTAAccccg                                                                                                                                           >  2:1319197/1‑147 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAAACCgg                                                                                                                                           >  2:216482/1‑149 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGGTAAccccg                                                                                                                                           >  2:507646/1‑147 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAG‑AATC‑A‑GTTTTTTT‑CCCGTTAAccccg                                                                                                                                           >  2:1075008/1‑147 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGGCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAAAG‑AATCAA‑GTTTTTTT‑CGGGTGAACCCgc                                                                                                                                           >  1:33620/1‑149 (MQ=255)
               atCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGGGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTCGAG‑TATT‑C‑GTTTTTCTCCCGGTTAACCCgc                                                                                                                                           >  1:160056/1‑149 (MQ=255)
                    gggTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AAAT‑A‑GGTTTTTT‑CCCGGTAAACCGGCCTCAAt                                                                                                                                    >  2:465363/1‑151 (MQ=255)
                      gTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATc                                                                                                                                   <  2:1282373/150‑1 (MQ=255)
                      gTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATc                                                                                                                                   <  2:819025/150‑1 (MQ=255)
                      gTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATc                                                                                                                                   <  1:1446487/150‑1 (MQ=255)
                           gcgcGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑c                                                                                                                             <  2:751622/150‑1 (MQ=255)
                                gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑C‑GTTTTTTT‑TCCGTTTACCCGGCTTCAATCAGTggccggta                                                                                                                        >  1:896420/1‑144 (MQ=255)
                                gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGAt                                                                                                                        <  1:730729/150‑1 (MQ=255)
                                gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGAt                                                                                                                        <  2:1324031/150‑1 (MQ=255)
                                gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGAt                                                                                                                        <  1:83362/150‑1 (MQ=255)
                                gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGAt                                                                                                                        <  1:1198022/150‑1 (MQ=255)
                                 cAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATag                                                                                                                      <  2:1360566/151‑1 (MQ=255)
                                  aaTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑AATC‑A‑GGTTTTTT‑CCCGGTAAACCCGCTTCAATCAATG‑CCGGGTaga                                                                                                                     >  1:338406/1‑151 (MQ=255)
                                   aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAAACCGGCTTCAATCAGTG‑CCGGGTagag                                                                                                                    >  2:176942/1‑151 (MQ=255)
                                   aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑TCGGCTAACCCGGCTTCAATCAGTG‑CCGGATagag                                                                                                                    >  2:268669/1‑151 (MQ=255)
                                   aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATA‑A‑GTTTTTTT‑TCGGGTAACCCGGCGTCCATCAGAG‑CCGGGAagag                                                                                                                    >  1:306527/1‑151 (MQ=255)
                                   aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATA‑A‑GGTTTTTT‑CCCGTTTAACCGGGTTCAAACAGTG‑CCGGGTagag                                                                                                                    >  2:508559/1‑151 (MQ=255)
                                    tGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGGTAgaga                                                                                                                   >  2:1487529/1‑151 (MQ=255)
                                     gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑C‑GTTTTTTT‑CCCGTTAACCCCGCCTCAATCAGTG‑CCGGATAGAGAt                                                                                                                  >  2:396282/1‑151 (MQ=255)
                                        aTACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAACCCGGCCTCAATCAGTG‑CCGGATAGAGAtttt                                                                                                               >  2:489574/1‑150 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGGTAACCCCGCTTCAATCAGTG‑CCGGGTAGAGAttttc                                                                                                              >  2:1165612/1‑148 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGcc                                                                                                             <  1:123758/151‑1 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGcc                                                                                                             <  1:1205049/151‑1 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGcc                                                                                                             <  2:572059/151‑1 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAAAG‑AATC‑A‑GTTTTTTT‑CCGGTTAACCCCGCTTCCATCAGTG‑CCCGGTAGAGAatttt                                                                                                              >  1:1383161/1‑148 (MQ=255)
                                          aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAAAG‑AATT‑C‑GTTTTTCT‑CTCGGTAACCCCGCGTTTATGAGGG‑CCCGGTAAAGAttttt                                                                                                              >  1:1138414/1‑148 (MQ=255)
                                             gtgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGAG‑AATC‑C‑GGTTTTTT‑TCCGTTAAACCCGCCTCACTCCGTG‑CCGGGTAGAGAATTGTCAc                                                                                                           >  1:1409269/1‑150 (MQ=255)
                                                cGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAAG‑G‑AATC‑A‑GTTTTTCT‑TCGGTTTACCCCGGTTCAAACACGG‑TCGGGTAGGGATTTTCCAcccc                                                                                                        >  2:770157/1‑148 (MQ=255)
                                                   ggggTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATG‑C‑GATTTTTT‑CCCGTTAACCCGGCTTCAAACCGTG‑CCGGGTAGAGATTTGCCAcccccgcg                                                                                                    >  1:1048707/1‑149 (MQ=255)
                                                    gggTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGGTAGAGATTTGCCACCGCCGTg                                                                                                    >  2:1215695/1‑150 (MQ=255)
                                                     ggTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AAAC‑A‑GTTTTTCT‑CCCGGTAAACCCGCTTCAATCAAGG‑CCGGGTAGAGATTTTCCAcccccccgc                                                                                                   >  2:630642/1‑146 (MQ=255)
                                                          gagGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAG‑GATC‑A‑GTTTTTTT‑CCGGGTAAACCCGCCTCCATCCGTG‑GCCGATAAAGAATTTCCACCCCCGTGcccgggg                                                                                             >  1:867031/1‑146 (MQ=255)
                                                                  cATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTACCCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATc                                                                                     <  1:1486164/151‑1 (MQ=255)
                                                                  cATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATc                                                                                     <  2:116851/151‑1 (MQ=255)
                                                                                          tttCGCGCCGCGGTCAATTAGCCCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCC‑CCGGTTACCCCGGCTTCAATCAGTC‑CCGGATAGAGATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACt                                                             <  2:841756/150‑1 (MQ=255)
                                                                                                   gcAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGGTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTTCCACCCCCGTGCCGGGGGTCACAATCACCCCCCCCCGCCGCCTGCTGACTTTTGcccg                                                     >  2:905844/1‑148 (MQ=255)
                                                                                                      gTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑AATC‑A‑GTTTTTTT‑TCGGGTTACCCGGCCTCAATCAGTG‑CCGGGTAGAGATTTTCCACCGCCCTGCCGGAGGGCACAATCAACGCCACCCGCTGCTTTCTGACTTTTGCCAGccc                                                  >  1:1503122/1‑150 (MQ=255)
                                                                                                      gTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAG‑AATC‑A‑GTTTTTTT‑CCGGTTAACCCGGCTTCAATCAGAG‑CCGGATAGAGAATTTCCACCGCCGGGCCCGAGGGCACAATCACCCCCACCCGCTGCTTGCTGGCTTTCCCCCGccc                                                  >  1:318623/1‑150 (MQ=255)
                                                                                                      gTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCGGGTAACCCGGCTTCAATCCGTG‑CCGGGTAGAGAATTTCCACCGCCGTGCCGGAGGGCACAATCACCGCCACCGGCCGCTTGCTGGCTTTTGCCAGccc                                                  >  1:124470/1‑150 (MQ=255)
                                                                                                        cAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTTT‑CCCGGTTACCCGGGTTCAAACAGTG‑CCGGGTAAAGATTTTCCACCCCCCTGCCGGAGGTCACAATCACCGCCACCCGCCGCCTGCTGACTTTTGCCAGGGCCAg                                               >  1:811260/1‑151 (MQ=255)
                                                                                                        cAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAAAA‑AATC‑A‑GTTTTTTT‑CCGGGTAAACCCGCTTCAATCACTG‑CCCGGTAGAGATTTGCCACCCCCGTGCCCGAAGTCACAAACACCCCCCCCCGCGTCCTGCCGACCTTTGCCCCCCCCa                                                >  1:1126135/1‑150 (MQ=255)
                                                                                                           ttacgtgcggcggGCGCTCGGC‑GGTTAAGGA‑AATC‑A‑GTTTTTCT‑CCGGTTAACACGGCTTCAATC‑GTGTCTGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAACCACCGCCTCCGGCTGCTTGCTGACTTTCGCCAGCCCCAgcgc                                            <  1:857553/144‑1 (MQ=255)
                                                                                                                     cgggcgATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGccc                                   <  1:1096649/150‑1 (MQ=255)
                                                                                                                               gc‑ggTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCCATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGGCACAATCACCGCCACCGGCTGCTTGCTGGCTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc                         >  1:1025662/1‑150 (MQ=255)
                                                                                                                                   gTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCAt                     >  1:1167759/1‑151 (MQ=255)
                                                                                                                                   gTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCAt                     >  1:1269024/1‑151 (MQ=255)
                                                                                                                                         ag‑aATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGaa                 <  1:1372970/149‑1 (MQ=255)
                                                                                                                                          g‑aATC‑A‑GTTTTTCT‑CCGGTTACCCCGGCTTCAATCAGTG‑CCGGTTAGAGATTTGCCACCCCCGCGCCGGAGGACACAATCATCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCCCGCTCATCGAAATg              >  1:548482/1‑151 (MQ=255)
                                                                                                                                            aaTC‑A‑GTTTTCCT‑CCGGTTACCCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATg              <  1:328398/150‑1 (MQ=255)
                                                                                                                                                   gTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAg                                                                                              >  2:12312/1‑65 (MQ=255)
                                                                                                                                                   gTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGtgtg        >  1:145706/1‑151 (MQ=255)
                                                                                                                                                   gTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCAACCGCCGTGCCGGAg                                                                                              <  1:12298/65‑1 (MQ=255)
                                                                                                                                                    tttttCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGgtg          <  2:1183135/148‑1 (MQ=255)
                                                                                                                                                    tttttCC‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGgtg          <  2:1190029/148‑1 (MQ=255)
                                                                                                                                                    tttttCC‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGgtg          <  2:599778/148‑1 (MQ=255)
                                                                                                                                                    tttttCC‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTgtgg       <  1:176889/151‑1 (MQ=255)
                                                                                                                                                     ttttCC‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTgtggt      <  2:1503909/151‑1 (MQ=255)
                                                                                                                                                         ct‑cCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGtgaat  <  2:1383848/151‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                        
CCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAG‑AATC‑A‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTG‑CCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAAT  >  NC_000913/2378819‑2379118

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