Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,547,7720TG40.6% 105.9 / 19.9 65V367G (GTA→GGA) murPN‑acetylmuramic acid permease, EIIBC component, PTS system
Reads supporting (aligned to +/- strand):  ref base T (11/27);  new base G (26/0);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.65e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCC  >  NC_000913/2547623‑2547922
                                                                                                                                                     |                                                                                                                                                       
gggCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTa                                                                                                                                                        <  1:715898/151‑1 (MQ=255)
     ttATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTAcg‑c                                                                                                                                                    <  2:900205/149‑1 (MQ=255)
       aTTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑Cgggg                                                                                                                                                <  2:876995/151‑1 (MQ=255)
           cTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCAACGCACAGTGGGCCCCCCCGTCCGGGACG‑CGGGGCg                                                                                                                                              >  1:425885/1‑149 (MQ=255)
             gTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGatta                                                                                                                                          <  2:286221/151‑1 (MQ=255)
               ttACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCACAGGGCGCTACCCAGGCAGGGACG‑Cggggggggtatt                                                                                                                                        >  1:867171/1‑143 (MQ=255)
                   cTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCCACGCAGGCCGGGACG‑CGGGGGGGTTATTtcc                                                                                                                                     >  1:555593/1‑150 (MQ=255)
                      gcgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGGCACCCCACCATGCGCCACGCAGGCACGGACG‑CGGGGGGGTTTTTTCCCGcc                                                                                                                                 >  1:467977/1‑151 (MQ=255)
                        gTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCctg                                                                                                                               <  1:872910/151‑1 (MQ=255)
                         ttAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCAACCCCCAATGCGCCACGCCGGCAGGGAAG‑CGGGGGGGTTATTTCCCGcccg                                                                                                                               >  1:484971/1‑148 (MQ=255)
                           aaTGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCctgctg                                                                                                                            <  2:328210/151‑1 (MQ=255)
                           aaTGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCctgctg                                                                                                                            <  1:342531/151‑1 (MQ=255)
                                  aGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCCACCGCCCCGCGCGCCACGCCGGCCGGGACG‑CCGGGCGGTTATTTCCCGCCTGCTGGGGGTTg                                                                                                                     >  1:918590/1‑151 (MQ=255)
                                     cAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTg                                                                                                                  <  1:934110/151‑1 (MQ=255)
                                     cAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTg                                                                                                                  <  2:470555/151‑1 (MQ=255)
                                     cAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGGGGGCTCCCCAGTCAGGGAAG‑GGGGGGGGTTTTTTCCCGCCCGCTgg                                                                                                                           >  1:530012/1‑142 (MQ=255)
                                      aGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCGCACAGGGCGCCAAGCCGGCAGGGACG‑CGGGGGGGGTTTTTCCCGCCCGCCGGGCGTTgg                                                                                                                    >  1:97640/1‑148 (MQ=255)
                                             cAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAATGCGCTACGCAGTCAGGGACG‑CGGGGCGATTTTTTCCGGCCCGCTGGGCGTTGGTGAAcc                                                                                                              >  2:807871/1‑147 (MQ=255)
                                                  gCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCAACCGCACAGGGCGCTACGCCGTCCGGGACG‑CGGGGGGGTTATTTCCGGCCTGCTGGGGGGTGGTGAAACGCCGAttt                                                                                                      >  1:837296/1‑150 (MQ=255)
                                                   cTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAAACGCACAGTGCGCTACGCAGGCAGGTACG‑CGGGGCGGGTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGAttt                                                                                                      >  2:593397/1‑149 (MQ=255)
                                                    ttATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCTCCCAGAGCGCTACGCCGTCAGGGGCGCCGGGGGGATATTTCCCGGCC‑GCCGGGCGGTGGGGGACCGCTGGttt                                                                                                      >  2:656545/1‑148 (MQ=255)
                                                    ttATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCACAGGGCGCTACGCCGGCAGGGACG‑CGGGGGGGTTATTTCCGGCCTGCTGGGCGGTGGTGAACCGCTGATTTACg                                                                                                   >  2:485546/1‑151 (MQ=255)
                                                          cccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCAAACGCAACGTGCGCGACGCAGGCCGGCGCG‑ACGGGGGATTATTTCCCGCCTGGTCGGCGGTgggggaaccctggatt                                                                                                      >  1:474174/1‑128 (MQ=255)
                                                               ccTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCCCCGCACAGTGCGCCACGCCGGCCGGGACG‑GCGGGCGGTTTTTTCCGGCCCGCTGGGGGTTGGGGGACCGCCGAATTTTGGTGTGAccc                                                                                          >  1:446876/1‑149 (MQ=255)
                                                               ccTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCt                                                                                         <  2:536289/150‑1 (MQ=255)
                                                                cTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCAACCGCCCAGTGCGCCACGCAGTCAGGGACG‑CGGGGCGGTTTTTTCCCGCCTGCTGGGCGGTGGGGGACCGCCGGTTTTCGGTGTCAcccc                                                                                         >  1:838672/1‑148 (MQ=255)
                                                                   tCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGGACG‑CGGGGCGGTTTTTTCCGGCCTGCTGGGCGGTGGTGGACCGCTGGTTTTCGGGGTCACCCCGcccc                                                                                    >  2:953678/1‑151 (MQ=255)
                                                                    cAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCg                                                                                   <  1:775113/151‑1 (MQ=255)
                                                                      aTGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCAACCCCACAGTGCGCTACGCAGTCAGGGACG‑CGGGGGGGTTATTCCCGGCCTGCTGGGCGTTGGTGGACCCCCGATTTACGGGGTCACCCTGcccc                                                                                    >  2:553531/1‑148 (MQ=255)
                                                                           aGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCTGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAg                                                                            <  2:552138/151‑1 (MQ=255)
                                                                           aGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCACAATGCGCCACCCAGTCAGGGACG‑CGGGGCGGTTATTTCCGGCCTGCTGGGGGGTGGTGGAACGCTGGTTTAAGGGGGCACCCCGCCCCCCCTGaa                                                                             >  2:494711/1‑150 (MQ=255)
                                                                            ggCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCCCAGAGCGCTACGCAGTCAGGGGCG‑CGGGGGGGTTTTTTCCGGCCCGCTGGGCGTTGGGGGACCGCCGGTTTTCGGGGTCACCCTGCCCCCCCTGAAg                                                                            >  2:338517/1‑150 (MQ=255)
                                                                               gcgGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCCGTCAGGGACG‑CGGGGCGGTTATTTCCCGCCTGCTGGGCGGTGGTGGAACGCTGATTTTCGGTGTCACCCTGCCCCGCCTGAAGcc                                                                          >  2:228980/1‑149 (MQ=255)
                                                                                  ggCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGGACG‑CGGGGGGATTATTCCCCGCCTGCTGGGCGGTGGTGAAACGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGgtt                                                                      >  2:499956/1‑150 (MQ=255)
                                                                                  ggCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGGACG‑CGGGGCGATTATTCCCGGCCCGCTGGGCGGTGGTGAACCGCTGATTTTCGGGGTCACCCTGCCCCCCATGAAACCGttt                                                                      >  2:941519/1‑150 (MQ=255)
                                                                                     cAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCGCACAGGGCGCTACGCCGTCAGGGAAG‑CGGGGCGGTTATTTCCCGCCCGCTGGGGGTTGGGGGAACGCTGATTTACGGTGTCCCCCCGCCCCGCCACAAGCCGGTtttt                                                                   >  1:733018/1‑150 (MQ=255)
                                                                                      aGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCACAGTGCGCTACGCGGTCCGGGAAG‑CGGGGCGGTTATTCCCCGCCTGCTTGGAGGTGGTGGACCGCTGATTTAACGTGTCACCCTGCCCCGCCTAAAGCCGgttgtt                                                                   >  1:234526/1‑149 (MQ=255)
                                                                                             cgcggcgcTGGCACTCTACTGGCGGGCGCCACCCCACACCCCACTACGCAGGCAGGGAAG‑CGGGGGGATTATTCCCCGGCTGCTGGGCGGTGGGGAACCGCTGGGTTACGGTCTGCCCCTGCCCCCCATGAAAGCGTTTGtt                                                                   >  1:142477/1‑142 (MQ=255)
                                                                                                ggcgcTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCCGTCAGGGACG‑CGGGGCGATTATTTCCCGCCTGCCGGGCGTTGGGGAACCGCTGATTTTCGGTGTCACCCTGCCCCGCCTGAAGCCGTTTGTTACCGCCTGttt                                                        >  1:671472/1‑150 (MQ=255)
                                                                                                  cgcTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGGCAGGGACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGGAACGCTGATTTACGGTGTCACCCTGCCCCGCATTAAGCCGTTTTTTAACGCCTGTTTTgg                                                     >  2:562409/1‑151 (MQ=255)
                                                                                                      ggCCCTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAggcgg                                                  <  2:704635/150‑1 (MQ=255)
                                                                                                      ggCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGCTTAggcg                                                   <  2:724117/149‑1 (MQ=255)
                                                                                                          ctctACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGgcggcgcg                                              <  2:241102/150‑1 (MQ=255)
                                                                                                                 ggcgggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCggcgggcg                                       <  2:867181/150‑1 (MQ=255)
                                                                                                                   cgggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGCTTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGgcgggcgg                                      >  1:364088/1‑149 (MQ=255)
                                                                                                                    gggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGgttt                                   <  2:162141/151‑1 (MQ=255)
                                                                                                                      gcgcAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGgtttgt                                 <  2:749093/151‑1 (MQ=255)
                                                                                                                       cgcAACCGCACAGTGCGCTACGCAGTCAGGTGCG‑CGGGGCGATTATTCCCGGCCTGCTGGGGGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCATGTTTAGGCGGCGCGGCGGGCGgtttgtt                                <  2:901257/151‑1 (MQ=255)
                                                                                                                       cgcAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGgtttgtt                                <  2:908030/151‑1 (MQ=255)
                                                                                                                         cAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTa                              <  1:97318/151‑1 (MQ=255)
                                                                                                                          aaCCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTAt                             >  2:587589/1‑151 (MQ=255)
                                                                                                                           aCCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATc                            <  2:708674/151‑1 (MQ=255)
                                                                                                                               cacaGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGcc                        >  1:91923/1‑151 (MQ=255)
                                                                                                                                 caGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTg                      <  2:50622/151‑1 (MQ=255)
                                                                                                                                  aGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGa                     <  2:1437/151‑1 (MQ=255)
                                                                                                                                    tGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGAt                    >  2:945319/1‑150 (MQ=255)
                                                                                                                                       gcTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATa                   <  2:61244/148‑1 (MQ=255)
                                                                                                                                       gcTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGcc                <  1:606440/151‑1 (MQ=255)
                                                                                                                                        cTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCt               >  2:472390/1‑151 (MQ=255)
                                                                                                                                         tACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCtg              >  2:808033/1‑151 (MQ=255)
                                                                                                                                                tCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTgggg        >  1:479434/1‑150 (MQ=255)
                                                                                                                                                 cAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTgggg        >  1:91941/1‑149 (MQ=255)
                                                                                                                                                 cAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGt       >  2:494188/1‑150 (MQ=255)
                                                                                                                                                 cAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTc      <  1:257582/151‑1 (MQ=255)
                                                                                                                                                    gTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGc   >  2:669747/1‑151 (MQ=255)
                                                                                                                                                     tACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGcc  >  2:872056/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACG‑CGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCC  >  NC_000913/2547623‑2547922

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