Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,450 A→G 27.4% intergenic (+55/‑123) 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,4500AG27.4% 144.9 / 61.5 93intergenic (+55/‑123)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (18/7);  total (54/39)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

TTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAG  >  NC_002947/5346301‑5346599
                                                                                                                                                       |                                                                                                                                                      
tttACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTa                                                                                                                                                        <  1:2909749/150‑1 (MQ=255)
    ccAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                                    <  2:4214809/150‑1 (MQ=255)
     cAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCCGGCCACGCGGTACCTCg                                                                                                                                                   >  2:4276549/1‑150 (MQ=255)
       gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                                                 <  2:4999689/150‑1 (MQ=255)
       gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                                                 <  2:1125664/150‑1 (MQ=255)
       gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                                                 <  2:405364/150‑1 (MQ=255)
           ttGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAg                                                                                                                                             <  2:3116332/150‑1 (MQ=255)
            tGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggc                                                                                                                                              <  1:2127632/150‑2 (MQ=37)
            tGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggc                                                                                                                                              <  1:2127880/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:2474137/150‑1 (MQ=255)
                    aGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTAc                                                                                                                                   <  1:1637955/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:4595649/150‑1 (MQ=255)
                           aaGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                            <  1:2119736/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:3782608/1‑150 (MQ=255)
                                                                                               tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCg                                                        >  2:4233179/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)
                                                                                                                               cccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGtt                        >  2:5021532/7‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:3764487/6‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:583479/6‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:2715018/6‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:3903221/6‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:3973247/6‑150 (MQ=25)
                                                                                                                                ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:3764560/6‑150 (MQ=25)
                                                                                                                                ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  1:1138628/1‑150 (MQ=255)
                                                                                                                                ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  2:5031766/1‑150 (MQ=255)
                                                                                                                                ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  1:443676/1‑150 (MQ=255)
                                                                                                                                 cgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       <  2:2949720/145‑1 (MQ=25)
                                                                                                                                 cgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                       >  1:2949720/5‑149 (MQ=25)
                                                                                                                                 cG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCg                      >  1:3818789/1‑150 (MQ=255)
                                                                                                                                  g‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAAGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGa                     >  2:2943594/1‑150 (MQ=255)
                                                                                                                                       tACGAGGTAACGCGTTGCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGatat                  >  1:4566149/1‑150 (MQ=255)
                                                                                                                                       tACGAGGTAACGCGTTGCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGatat                  >  1:3579071/1‑150 (MQ=255)
                                                                                                                                       tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATcc                >  1:3537429/1‑150 (MQ=37)
                                                                                                                                       tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATcc                >  1:2180630/1‑150 (MQ=37)
                                                                                                                                       tACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGatat                  >  1:4732860/1‑150 (MQ=255)
                                                                                                                                       aaCCCGGCACCGCGTTACCGCGTCGGTGCGGTTTTACCCGGGAAGCGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGatat                  <  2:91507/149‑1 (MQ=255)
                                                                                                                                          gaGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa             >  1:4150372/2‑150 (MQ=38)
                                                                                                                                           aGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTATTATCGAATAGCCGTTGCAGTCCAGCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCga              >  1:3881040/1‑149 (MQ=255)
                                                                                                                                            ggtaaCGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCaaa           >  1:1216573/4‑150 (MQ=38)
                                                                                                                                            ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCACCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa             >  2:2858888/1‑150 (MQ=255)
                                                                                                                                             gCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCaa            >  1:127332/1‑150 (MQ=255)
                                                                                                                                             gCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCaa            >  2:2258029/1‑150 (MQ=255)
                                                                                                                                              cAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCaaa           >  1:3078560/1‑150 (MQ=255)
                                                                                                                                                aCGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAcc       >  2:403706/1‑150 (MQ=255)
                                                                                                                                                 cgcgTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAccc      >  1:3094830/1‑150 (MQ=255)
                                                                                                                                                 cgcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTACAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAc        >  1:2128024/1‑150 (MQ=255)
                                                                                                                                                 cgcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTACAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAc        >  1:4698807/1‑150 (MQ=255)
                                                                                                                                                  gcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAcc       >  2:4731287/1‑150 (MQ=255)
                                                                                                                                                       aCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAg  >  1:155905/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAG  >  NC_002947/5346301‑5346599

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: