Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,431:1 +C 9.2% intergenic (+36/‑142) pgi‑II → / → acsA‑II glucose 6‑phosphate isomerase/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,346,4311.C 9.2% 309.6 / 26.1 76intergenic (+36/‑142)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base . (29/40);  new base C (2/5);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

AGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGA  >  NC_002947/5346283‑5346580
                                                                                                                                                     |                                                                                                                                                       
aggagctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑c                                                                                                                                                       <  2:3963969/150‑1 (MQ=255)
aggagctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑c                                                                                                                                                       <  2:199695/150‑1 (MQ=255)
 ggagctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAAcccc‑‑cc                                                                                                                                                      >  2:567902/1‑150 (MQ=255)
    gctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTGGCGGGGTAACCCG‑‑CCGAc                                                                                                                                                   >  1:4351650/1‑150 (MQ=255)
      tgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCa                                                                                                                                                 <  2:140468/150‑1 (MQ=255)
      tgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCa                                                                                                                                                 >  2:4671484/1‑150 (MQ=255)
       gggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAg                                                                                                                                                <  1:3032307/150‑1 (MQ=255)
              ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         <  1:1106822/150‑1 (MQ=255)
              ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         >  2:2880784/1‑150 (MQ=255)
              ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         <  2:4541831/150‑1 (MQ=255)
              ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         <  2:4563408/150‑1 (MQ=255)
              ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         <  2:5577/150‑1 (MQ=255)
                 gTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGtt                                                                                                                                      <  1:1904654/150‑1 (MQ=255)
                  tttACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTa                                                                                                                                     <  1:2909749/150‑1 (MQ=255)
                      ccAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                 <  2:4214809/150‑1 (MQ=255)
                       cAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCCGGCCACGCGGTACCTCg                                                                                                                                >  2:4276549/1‑150 (MQ=255)
                        aGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑c                                                                                                                                                       <  2:1069151/126‑1 (MQ=255)
                        aGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑c                                                                                                                                                       >  1:1069151/1‑126 (MQ=255)
                         gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                              <  2:1125664/150‑1 (MQ=255)
                         gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                              <  2:4999689/150‑1 (MQ=255)
                         gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                              <  2:405364/150‑1 (MQ=255)
                             ttGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAg                                                                                                                          <  2:3116332/150‑1 (MQ=255)
                              tGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggc                                                                                                                           <  1:2127880/150‑2 (MQ=37)
                              tGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggc                                                                                                                           <  1:2127632/150‑2 (MQ=37)
                              tGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACGAGGTAACGCGTTGCCTGGTAggcgg                                                                                                                         <  2:1605492/150‑4 (MQ=37)
                               ggTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCg                                                                                                                       <  2:437113/150‑1 (MQ=255)
                                gTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgg                                                                                                                         <  1:3764174/150‑4 (MQ=25)
                                gTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgg                                                                                                                         <  1:3763999/150‑4 (MQ=25)
                                 tCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcggg                                                                                                                        <  2:3131236/150‑5 (MQ=25)
                                  cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGt                                                                                                                    <  1:581910/150‑1 (MQ=255)
                                    gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGttt                                                                                                                  <  2:2094102/150‑1 (MQ=255)
                                      aGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTAc                                                                                                                <  1:1637955/150‑1 (MQ=255)
                                      aGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTAc                                                                                                                <  1:2474137/150‑1 (MQ=255)
                                       gCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTAcc                                                                                                               <  2:3291005/150‑1 (MQ=255)
                                         cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCGCGTAGCGGGGGAACCCC‑‑CCGACCACGCAACCCGGTACCACGTGGGCGCGGGGTTAcccc                                                                                                             >  1:3167070/1‑149 (MQ=25)
                                            gAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTAAACCcgctc                                                                                                          <  2:3055497/150‑3 (MQ=255)
                                            gAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACGTCCGCGGCCGCCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTAAACCcgctc                                                                                                          <  2:3055651/150‑3 (MQ=255)
                                             aaGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                         <  1:2119736/150‑1 (MQ=255)
                                             aaGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                         <  1:4595649/150‑1 (MQ=255)
                                             aaGACAGCGCCGAAGACGGGTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                         <  2:1273835/150‑1 (MQ=255)
                                               gACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                       <  1:4242483/150‑1 (MQ=255)
                                                   gcgcCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCAGTCACCGCGGTTTATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACATCGTAGGAGCGGGTTTACCCGCGaatattcc                                                                                                   <  2:4827326/150‑7 (MQ=39)
                                                    cgcCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                                                  <  2:1576078/150‑1 (MQ=255)
                                                         aaGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTTGCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCa                                                                                             >  2:4303647/1‑150 (MQ=255)
                                                                ttCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAAGCCC‑‑CCGACCAGGCCACGCGGTACCGCGGAGGGGCGGGGTTTACCGCCAAGATGCCGGCCCACGTTGaa                                                                                      >  2:4420/1‑150 (MQ=255)
                                                                ttCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAAACCC‑‑CCAACCAGGCAAACCCTTAACACGGAGGAGCGGGTTTACCCGCGAAGAGGCCCGCCCAGGTTTaa                                                                                      >  2:4248/1‑150 (MQ=255)
                                                                ttCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCGTCGCGGGTAAACACG‑‑CCTAACAGGCAAAGCGTTAACTCGTTGGAGCGGGTTTAACCGCGAAGAGGCCCGCCCAgggttaa                                                                                      >  2:4239/1‑145 (MQ=255)
                                                                  ccACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAaca                                                                                    <  2:4227515/150‑1 (MQ=255)
                                                                  ccACCCAAGGCCTGATCAACTACTTCCGCGGCCGCCCCCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAACCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTGAaca                                                                                    <  2:4227911/150‑1 (MQ=255)
                                                                      ccAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCGTCGCGGGTAAACCCC‑‑CCGCCCAGGCAACCCGGTACCGCGGAGGAGCGGGGTTAACCCCGCAGCGGCCGGCCCCGGGTTGACACCGc                                                                                >  2:2317604/1‑150 (MQ=255)
                                                                       cAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCC‑‑CCTACCACGCAACCCGTTACCACGTAGGTGCGGGGTTTACCCCGAAGAGGCCGGCCCAGGTTGGACACCGCt                                                                               >  1:3731383/1‑150 (MQ=255)
                                                                              tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCAACCAGCt                                                                        >  2:4152594/1‑150 (MQ=255)
                                                                                      ttcgcccGGGGCCGCCCCGGGGGTTGAGCCGTGTTTGTGCGGGCATTTTCGCTGTAAAACCCC‑‑CATCCCTGGCAACGCGTAACCTCTGTGGCGGGGGTTCACCCGCGAAGCGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTc                                                                <  2:504461/145‑1 (MQ=255)
                                                                                                     cACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGt                                                 <  2:3772042/150‑1 (MQ=255)
                                                                                                               gATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGGAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCAt                                       >  2:2599288/1‑150 (MQ=255)
                                                                                                                 tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCCACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCg                                     >  2:4233252/1‑150 (MQ=255)
                                                                                                                 tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCg                                     >  2:4233179/1‑150 (MQ=255)
                                                                                                                 tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCg                                     >  2:3782608/1‑150 (MQ=255)
                                                                                                                 tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCC‑‑CCGAACCCGCAACCCCTGAACTACTACGACCCCCATAACCCCCGAAGAGCCCGGCCCCAGTTGAAAACCGCGCCCACCCGCACAGGCTTATCGAAAAGACGATTCAGACCAACg                                     >  2:3277262/1‑150 (MQ=255)
                                                                                                                  cggtgCTTGTGCGGGCCTCTTCGCGGGAAAACCCG‑‑CCAACCCGGCAACGCTTTACCTCGTAGGAGGGGTTTTCCCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGc                                    <  2:3611686/146‑1 (MQ=255)
                                                                                                                        ttGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCCGGCAACGCTTTACCTCTTAGGAGCGGTTTTACCCGCTATTAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTa                              <  1:1157855/150‑1 (MQ=255)
                                                                                                                         tgggCCGGCCTCTTCGGGGGTAACCCCG‑‑CCTACCAGCCAACGCTTTCCCTCGTAGGCGCGGGTTAACCCGCTAAGCGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTAc                             <  1:3374106/147‑1 (MQ=255)
                                                                                                                             ggggCCTGTTCGGGGGTAACCCCG‑‑CTTACCCGGCACGGCGTTACGTGTTTGGTGGGGGTTTACCCGGTAAGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAaga                         <  2:4400029/148‑1 (MQ=255)
                                                                                                                             cgggCTTCTTGGGGGTTAACCCCG‑‑CCAACCTGGAACCGGGTTACCTGTTTGTTGGGTTTTTACCGGCTATGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAaga                         <  2:930474/148‑1 (MQ=255)
                                                                                                                             ccGGCCTCTTCGCGGTTAAACCCG‑‑CCTACCTGCCAACGCGTTACCTCGTTGGAGGGGTTTTACCCGCTATGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAaga                         <  2:732957/150‑1 (MQ=255)
                                                                                                                                gCCTCTTGGGGTTTAACCCCG‑‑CCTACCCGGAACCGGTTTACGTCTTTGGTGGGGTTTTACCGGCTATGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCa                      <  1:2880784/150‑1 (MQ=255)
                                                                                                                                gCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCACCGCTTTACCTCTTAGGTGCGGGTTTACCGGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCa                      <  2:315151/150‑1 (MQ=255)
                                                                                                                                   tctTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGa                   >  2:1494965/1‑150 (MQ=255)
                                                                                                                                   tctTCGCGGGTAAACCCG‑‑CCTACCAGGCAAAGCGTTAACTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATACAATAGCCGTTGCAGTCCAACGCCCCATACAAAAGCAGGa                   >  1:3986590/1‑150 (MQ=38)
                                                                                                                                     ttCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGacc                 >  2:3515672/1‑150 (MQ=21)
                                                                                                                                          gggTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg            >  2:825801/1‑150 (MQ=21)
                                                                                                                                          gggTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACAGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg            >  1:3089481/1‑150 (MQ=255)
                                                                                                                                            gTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc          >  2:4771660/1‑150 (MQ=255)
                                                                                                                                            gTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc          >  1:486420/1‑150 (MQ=255)
                                                                                                                                            gTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc          >  2:384665/1‑150 (MQ=255)
                                                                                                                                                aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt      >  1:4791057/1‑150 (MQ=255)
                                                                                                                                                aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt      >  1:2387875/1‑150 (MQ=255)
                                                                                                                                                  ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    >  1:443676/1‑150 (MQ=255)
                                                                                                                                                  ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    >  2:5031766/1‑150 (MQ=255)
                                                                                                                                                  ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    >  1:1138628/1‑150 (MQ=255)
                                                                                                                                                   cG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCg   >  1:3818789/1‑150 (MQ=255)
                                                                                                                                                    g‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAAGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGa  >  2:2943594/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGA  >  NC_002947/5346283‑5346580

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: