Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,454 C→G 22.7% intergenic (+59/‑119) 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,4540CG22.7% 131.0 / 38.4 70intergenic (+59/‑119)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (33/20);  new base G (7/9);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

GCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAGGCCC  >  NC_002947/5346308‑5346603
                                                                                                                                                     |                                                                                                                                                      
gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                                                     <  2:5304752/150‑1 (MQ=255)
       gTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                              <  1:3448399/150‑1 (MQ=255)
         cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                              <  1:112124/150‑6 (MQ=25)
         cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                              <  1:3492351/150‑6 (MQ=25)
         cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGt                                                                                                                                            <  2:4575096/150‑1 (MQ=255)
          ggTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCAG‑‑CCTACCAGGCAACGCGTTACCTCGGAGGAGCGGGtt                                                                                                                                           >  2:4183504/1‑150 (MQ=255)
           gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                              <  1:2171408/146‑1 (MQ=255)
           gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                              >  2:2171408/1‑146 (MQ=255)
            tAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTa                                                                                                                                         <  2:1979654/150‑1 (MQ=25)
               ccTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAccc                                                                                                                                      <  2:2217872/150‑1 (MQ=21)
                cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                                        >  1:2761463/1‑131 (MQ=255)
                cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                                        <  2:2761463/131‑1 (MQ=255)
                 tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                                                    <  1:848547/150‑1 (MQ=255)
                 tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                                                    <  1:4587956/150‑1 (MQ=255)
                    aaGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                 <  1:5144359/150‑1 (MQ=255)
                      gACAGCGCCTAAGACGGTTCCACCCACGGCCTGATCAACTACTTCCGCGGCCGTCAGCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGTTAACGCGTTGCCTGGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                 <  2:3375274/150‑1 (MQ=21)
                      gACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACGTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                               <  2:83669/150‑1 (MQ=255)
                       aCAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCACTGGGCGGGGGAAACCCC‑‑CCACCCAGCCACCCCGTTACCACGAGGGGGGGGGTTTGCCCCCGacgg                                                                                                                               >  1:1891544/1‑147 (MQ=19)
                        cAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGTGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                             <  1:1677516/150‑1 (MQ=255)
                               gAAGCCGGTTCCACCCATGGCCTGATCAACTAGGGCGGGGGCCGTCACCGGGGTTGATCCGTGTTTGTGCCGGCCTCTGCGGGGGTAAACCCG‑‑CCTACCAGGCAAGGGGTTACCTGGTTGGTGGGGGTTTCCCCGCGAAGAGGCCGGccc                                                                                                                     <  2:4832859/150‑1 (MQ=255)
                               gAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCGCGTCGCGGGGGAACCCG‑‑CCGACCCGGCAACGCGTTACCTCGTGGGGGCGGGTTTTCCCGCGGAGGGGGCGGccc                                                                                                                     >  1:4823307/1‑150 (MQ=255)
                                aaGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACCCGTTACCAAGTAGGAGCCGGTTTACCCCCGAAGAGGCCGGCCCa                                                                                                                    >  1:3572136/1‑150 (MQ=39)
                                    cGGTTCCACCCATGGCATGATCAACTACTTCCGCGGCCGTCACCGGGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGAAACGCGTTACCTCGTAGGAGCGGTTTTACCCGCGAAGTGGCCGGCCCAGGtt                                                                                                                <  2:4979188/150‑1 (MQ=38)
                                    cGGTTCCACCCACGGCCTTATCAACTACTTCCGCGGCCGTCACGGCGGTTGATCCGAGCTTGTGCCGGCCTCTGCGCGGGTAAACCCG‑‑CCTACCCGGCAACGCGTTCCCTCGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGtt                                                                                                                <  2:5444794/150‑1 (MQ=38)
                                     ggTTCCACCCAAGGCCTTATCAACTACTTCCGCGGCCGTCACCGCGTTTGATCCGAGCTTGTGCCGGCATCTTCGCGGTTAAACCCG‑‑CCTACCAGGAAACGCTTTACCTCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTg                                                                                                               <  2:2244991/150‑1 (MQ=38)
                                             ccAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGGAAAACCG‑‑CCTACCAGGCCACCCCGTACCACGGAGGAGCGGGGTTTACCGCGAAGAGGCCCGCCCACGTTGAAAAcccc                                                                                                       >  1:886098/1‑148 (MQ=255)
                                                     tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                              >  1:4195297/1‑104 (MQ=255)
                                                     tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                              <  2:4195297/104‑1 (MQ=255)
                                                     gtaTCAACTCCTCCCGCGGCCGTCACCGGGGTTTATCCGCGCTTGGCCCGGCCACGTCGCGGTAAAACCCG‑‑CCAACCCGGCAACGGGTTACCTCGTAGGCGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCCCCCACCt                                                                                               <  2:3572136/148‑1 (MQ=255)
                                                      gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                              <  2:351423/105‑6 (MQ=17)
                                                      gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                              >  1:351423/1‑100 (MQ=17)
                                                      gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAAACCC‑‑CCGACCAGGCAACGCGGTTACACGTAGGCGCGGGGTTACCCCCGGAGAGGCCGGCCCAGGGTGAACACCGCGCCCAGcga                                                                                              >  1:2221786/1‑148 (MQ=255)
                                                          aaCTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACt                                                                                          >  2:2233920/1‑150 (MQ=255)
                                                                         cGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACAGTCTTATCGAATAGc                                                                           >  2:3623106/1‑150 (MQ=255)
                                                                         cGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCGTACCGGGTTAACCCC‑‑CCAACCAGGCCACCCCGTACCCCGTAGGGGCGGGGTTTCCCCCCAAGAGGCCCGGCCCGGGTTAACCCCGCTCCCACCGACACAGTCTTAAAGAATAGc                                                                           >  2:2959779/1‑150 (MQ=255)
                                                                                 cgGTTGATCCGAGCTTGTGCCGGCCTCGTCGCGGGTAAAACCG‑‑CCTACCAGGCAACGCGTTACCTCGGAGGACCGGGGTTACCCCCGAAGAGGCCGGCCCAGGGTGAACACCGCACCCAGCTACACAGTCTTATCGAAGAGCCGGTGCAg                                                                   >  2:4475626/1‑150 (MQ=255)
                                                                                        tCCGAGCTTGTGCCGGCCTCGTCGCGGGGAAAACCC‑‑CCGAACACGCCACCCGGTACCACCGAGGGGCGGGTTTAACCCCGCACAGGCCCGCCCCCGGTGTACACCCCACCCCCCAACACAGGCTTATCGAAAAGCCGGTGCAGTACaacc                                                            >  2:696842/1‑147 (MQ=255)
                                                                                               ttGTGCCGGCCTCTTAGCGGGTAAACCCC‑‑CCAACCAGGCAACGCCTTAACTCGTAGGAGCCGGTTTACCCCCGAAAAGGCCCGCCCAGGTTGAACACCGCTCCCAGCTACACTGTATTAACGAAAAACCGTTGCAGTCCATCGCCCCGTa                                                     >  2:1291994/1‑150 (MQ=255)
                                                                                                  ggCCGGCGTCTTCGGGGGTAACCCCG‑‑CCTACCTGGCACGGGGTTCCCTCGTTGGCGGGGGTTTACCGGGGACGGGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACaa                                                  <  2:3780820/149‑1 (MQ=255)
                                                                                                       gCCTCTGCGGGGTAAAACCCG‑‑CATACCTGGAAACGCTTTACCTCTTTGTTGGGTTTTTACCCGCTAAGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCa                                             <  2:1424298/150‑1 (MQ=255)
                                                                                                        ccTCTTCGGGGGTAAACCCG‑‑CCTACCAGGCAAGGCGTTACGTGGTTGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAg                                            <  2:3255662/150‑1 (MQ=39)
                                                                                                          tctTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAggc                                          >  2:3175721/1‑149 (MQ=255)
                                                                                                            ttCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGacc                                        >  2:2637770/1‑150 (MQ=255)
                                                                                                              cgcgGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGaccac                                      >  2:478230/1‑150 (MQ=255)
                                                                                                                  ggttAACCCCG‑‑CCTACCCGGCACCGCGTTACGTCGTCGGCGGGGTTTTCCCCGCTATGCGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg                                   <  1:2233920/147‑1 (MQ=255)
                                                                                                                  gggTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg                                   >  1:1368915/1‑150 (MQ=255)
                                                                                                                   ggTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGc                                  >  2:3348225/1‑150 (MQ=255)
                                                                                                                    gTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                                 >  1:4332888/1‑150 (MQ=21)
                                                                                                                    gTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                                 >  2:1012599/1‑150 (MQ=255)
                                                                                                                    ctAAACCCG‑‑CCACCCCGGAAACGGTTTACCTCGTTGGAGCGGGTTTACCGGGTATGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                                 <  2:4592796/149‑1 (MQ=255)
                                                                                                                        aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt                             >  1:4016776/1‑150 (MQ=255)
                                                                                                                        aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt                             >  1:1074599/1‑150 (MQ=255)
                                                                                                                         cccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGGGCAGGACCACCCGCCATGtt                            >  1:2888797/7‑150 (MQ=25)
                                                                                                                          ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           >  2:5310536/6‑150 (MQ=25)
                                                                                                                          ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           >  2:2437345/6‑150 (MQ=25)
                                                                                                                          ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           >  1:5385930/6‑150 (MQ=25)
                                                                                                                          ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           >  1:3693510/1‑150 (MQ=255)
                                                                                                                          ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           >  2:162846/1‑150 (MQ=255)
                                                                                                                          ccG‑‑CCAACCCGGCACCGCGTTCCCGCTTTGGCGCGGGTTCACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                           <  2:158979/150‑1 (MQ=255)
                                                                                                                            g‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGa                         >  2:4991957/1‑150 (MQ=255)
                                                                                                                                 tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTAc                                                    >  1:2460264/1‑118 (MQ=25)
                                                                                                                                 tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTAc                                                    <  2:2460264/118‑1 (MQ=25)
                                                                                                                                  aCCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATc                     >  2:411124/1‑150 (MQ=255)
                                                                                                                                    cAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCg                   >  1:4805984/1‑150 (MQ=255)
                                                                                                                                      ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa                 >  2:768766/1‑150 (MQ=255)
                                                                                                                                      ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa                 >  2:768793/1‑150 (MQ=255)
                                                                                                                                      ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa                 >  2:768797/1‑150 (MQ=255)
                                                                                                                                          aCGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATa             >  2:4027394/1‑150 (MQ=255)
                                                                                                                                           cgcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAc            >  1:3458223/1‑150 (MQ=255)
                                                                                                                                            gcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAcc           >  1:3702099/1‑150 (MQ=255)
                                                                                                                                             cgTTACCTCGTCGGTGGGGGTTTACCCGCGACGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAccc          <  2:4823307/150‑1 (MQ=255)
                                                                                                                                               ttACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAccccc        >  2:2128269/1‑150 (MQ=255)
                                                                                                                                                     cGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAGGccc  >  2:2318966/1‑150 (MQ=255)
                                                                                                                                                     cGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAGGccc  >  1:4307285/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCCAGGCCC  >  NC_002947/5346308‑5346603

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: