Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,431:2 +T 12.5% 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,4312.T12.5% 240.8 / 31.7 64intergenic (+36/‑142)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base . (26/30);  new base T (2/6);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

GAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGA  >  NC_002947/5346285‑5346580
                                                                                                                                                     |                                                                                                                                                      
gagctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCt                                                                                                                                                     <  2:347238/150‑1 (MQ=255)
   ctgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTAcc                                                                                                                                                  <  2:2456688/150‑1 (MQ=255)
    tgggcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCa                                                                                                                                                 <  1:1332168/150‑1 (MQ=255)
       gcaaggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGc                                                                                                                                              <  2:2751633/150‑1 (MQ=255)
          aggGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAc                                                                                                                                           <  2:3219351/150‑1 (MQ=255)
           gggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAc                                                                                                                                           >  2:3324894/1‑149 (MQ=255)
           gggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAc                                                                                                                                           <  1:3324894/149‑1 (MQ=255)
            ggCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAcgc                                                                                                                                         <  1:1463155/150‑1 (MQ=255)
                  tACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTAcc                                                                                                                                   <  1:144039/150‑1 (MQ=255)
                       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                                                                                                                       >  2:2171408/1‑146 (MQ=255)
                                  gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                       <  1:2171408/146‑1 (MQ=255)
                                   tAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTa                                                                                                                  <  2:1979654/150‑1 (MQ=25)
                                      ccTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAccc                                                                                                               <  2:2217872/150‑1 (MQ=21)
                                       cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                 <  2:2761463/131‑1 (MQ=255)
                                       cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                 >  1:2761463/1‑131 (MQ=255)
                                        tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                             <  1:4587956/150‑1 (MQ=255)
                                        tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                             <  1:848547/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                                                                                                                       <  2:4195297/104‑1 (MQ=255)
                                                                            tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                       >  1:4195297/1‑104 (MQ=255)
                                                                            gtaTCAACTCCTCCCGCGGCCGTCACCGGGGTTTATCCGCGCTTGGCCCGGCCACGTCGCGGTAAAACCCG‑‑CCAACCCGGCAACGGGTTACCTCGTAGGCGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCCCCCACCt                                                                        <  2:3572136/148‑1 (MQ=255)
                                                                             gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                       >  1:351423/1‑100 (MQ=17)
                                                                             gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                       <  2:351423/105‑6 (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)
                                                                                                                                                 ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    >  2:162846/1‑150 (MQ=255)
                                                                                                                                                 ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    >  1:3693510/1‑150 (MQ=255)
                                                                                                                                                 ccG‑‑CCAACCCGGCACCGCGTTCCCGCTTTGGCGCGGGTTCACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc    <  2:158979/150‑1 (MQ=255)
                                                                                                                                                   g‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGa  >  2:4991957/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGA  >  NC_002947/5346285‑5346580

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