Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,701,229 G→C 9.7% intergenic (+27/‑36) phaA → / → phaB poly(3‑hydroxyalkanoate) polymerase/poly(3‑hydroxyalkanoate) depolymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,701,2290GC 9.7% 175.6 / 9.9 62intergenic (+27/‑36)phaA/phaBpoly(3‑hydroxyalkanoate) polymerase/poly(3‑hydroxyalkanoate) depolymerase
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (3/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

CACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCAT  >  NC_002947/5701083‑5701374
                                                                                                                                                  |                                                                                                                                                 
cACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGgc                                                                                                                                                  <  6:286854/148‑1 (MQ=255)
  cTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGCAGGCCGCCACGgcgct                                                                                                                                               <  6:178661/149‑2 (MQ=255)
  cTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGCAGGCCGCCACGgcgct                                                                                                                                               <  7:211228/149‑2 (MQ=255)
     gCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACg                                                                                                                                            >  1:296743/1‑149 (MQ=255)
     gCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCAc                                                                                                                                             <  2:113719/148‑1 (MQ=255)
       agagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCAGCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGgt                                                                                                                                          >  2:206560/1‑149 (MQ=255)
       agagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCACggc                                                                                                                                          <  6:113903/149‑2 (MQ=255)
        gagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGAGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGgtg                                                                                                                                         <  8:76820/149‑1 (MQ=255)
              gCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACt                                                                                                                                   >  5:381788/1‑149 (MQ=255)
                tGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACtt                                                                                                                                  >  8:226396/1‑148 (MQ=255)
                  ggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAc                                                                                                                               <  4:18226/149‑1 (MQ=255)
                     gAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTAc                                                                                                                                    <  1:179975/141‑1 (MQ=255)
                      aGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCaa                                                                                                                           <  7:270865/149‑1 (MQ=255)
                        cGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                         <  6:381788/149‑1 (MQ=255)
                        cGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                         <  8:154817/149‑1 (MQ=255)
                          tGCCGGCGCGCTGAAAACGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGt                                                                                                                       <  2:307010/149‑1 (MQ=255)
                            ccGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCATCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTc                                                                                                                      <  1:206560/148‑1 (MQ=255)
                             cGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCAc                                                                                                                    <  7:226396/149‑1 (MQ=255)
                             cGGCGCGCTGAAAAAGGCACCGACCCGCCAGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGAGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCAc                                                                                                                    >  2:73126/1‑149 (MQ=255)
                                   gcTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCAt                                                                                                              >  1:342338/1‑149 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                             <  2:100857/132‑1 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                             >  1:100857/1‑132 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGc                                                                                                            >  2:171493/1‑149 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCACg                                                                                                                                            >  1:239362/1‑117 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCACg                                                                                                                                            <  2:239362/117‑1 (MQ=255)
                                     tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATACCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGc                                                                                                            >  1:211493/1‑149 (MQ=255)
                                             gCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTAACCCAAGAGTCACGTGCATGCCGCAAccc                                                                                                    >  4:90791/1‑149 (MQ=255)
                                                aagACCCTCCTGGGCAAACGTACCTATGCCGCCGGCTAAGCCTCCCCAGGCACCAACGTTCACGATCGTTGAGTTACCGCGCCGTGGAGGCCTTAACGGCGCCACGGTGTTTACTTAAACCAAGAGTCACGTTCATGCAGCAAACCta                                                                                                  >  4:251344/3‑148 (MQ=255)
                                                   aCCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGc                                                                                                         >  5:256719/1‑138 (MQ=255)
                                                   aCCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGc                                                                                                         <  6:256719/138‑1 (MQ=255)
                                                       gCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                          <  2:296743/149‑1 (MQ=255)
                                                        ccTGGGCAACCGTACCTATGGCGCCGGCGAAGCCTCCGCAGGCACCTACGTTCTCGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTGACGTGCATGCCGCAACCCTATATTTTCAg                                                                                         <  8:188914/149‑1 (MQ=255)
                                                             gCAACCGTACCTATGCCGCCGGCGAAGCCTTCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCg                                                                                    >  4:196468/1‑149 (MQ=255)
                                                               aaCCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCAGCAACCCTATATTTTCAGGACCGTc                                                                                  >  3:10330/1‑149 (MQ=255)
                                                                 ccGTACCTATACCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGa                                                                                <  2:211493/149‑1 (MQ=255)
                                                                      ccTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTg                                                                            <  2:334552/148‑1 (MQ=255)
                                                                      ccTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTg                                                                            <  7:261574/148‑1 (MQ=255)
                                                                         aTGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGAc                                                                         >  8:195779/1‑148 (MQ=255)
                                                                            ccgccgGCGAAGCCTTCCCAGGCACCTACGTTCACGAGCGTTGAGGTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAAc                                                                      <  3:196468/148‑1 (MQ=255)
                                                                            ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAAc                                                                      <  1:171493/148‑1 (MQ=255)
                                                                            ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAAc                                                                      <  7:195779/148‑1 (MQ=255)
                                                                             cgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCg                                                                                    <  1:232499/133‑1 (MQ=255)
                                                                             cgccgGCGAAGCCTACCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCg                                                                                    >  2:232499/1‑133 (MQ=255)
                                                                                  gCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCa                                                               >  7:120333/1‑149 (MQ=255)
                                                                                    gAAGCCTCCCCAGGCAGCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGGCGCAACCCTATATTTTCAGGACCTTCGAGCTGGACATCCAGTACATc                                                             >  8:151086/1‑149 (MQ=255)
                                                                                         cTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCAc                                                        <  1:276527/149‑1 (MQ=255)
                                                                                           ccccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCg                                                      >  2:287226/1‑149 (MQ=255)
                                                                                            cccAGGCCCCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGAGTTTACTTAACCCAAGAGTCACGTGAATGACGCAACCCGATATTTTCAGGACAGTCGAGCTGGACAACCAGTCCATCCGCACCGc                                                     >  8:235032/1‑149 (MQ=255)
                                                                                            cccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGACACGGTGTTTACTTCACCCAAGAGTAACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGc                                                     >  7:53236/1‑149 (MQ=255)
                                                                                             ccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCt                                                    >  1:47715/1‑149 (MQ=255)
                                                                                             ccAGGCACCTACGATCACGAGCGTGGAGTTCCAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCt                                                    <  2:40993/149‑1 (MQ=255)
                                                                                                   aCCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGc                                              >  3:34539/1‑149 (MQ=255)
                                                                                                             aCGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAg                                     <  1:357843/148‑1 (MQ=255)
                                                                                                             aCGACCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGATTCACGTGCAAGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAg                                     <  7:151086/148‑1 (MQ=255)
                                                                                                               gAGCGTTGAGATACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTGACCCAAGATTCACGTGCATGCCGCAACCCTTTATTTTCATGACCTTCGAGCTGGACAACCAGTCCATCCTCCTCGCTGTTCGCCCCGGCAAGCCg                                  <  7:235032/149‑1 (MQ=255)
                                                                                                                  cGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCAc                               <  1:287226/149‑1 (MQ=255)
                                                                                                                    ttGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCt                                                                             >  5:360327/1‑101 (MQ=255)
                                                                                                                    ttGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCt                                                                             <  6:360327/101‑1 (MQ=255)
                                                                                                                    ttGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCt                             <  4:34539/149‑1 (MQ=255)
                                                                                                                     tGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCt                             >  6:270626/1‑148 (MQ=255)
                                                                                                                          gtaaacaccGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGCCAACCAGTCCCTCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGcc                       >  5:357152/8‑149 (MQ=255)
                                                                                                                             aacaccGTGGCGCCGTGCAGGCCGCCACGGCGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGtt                    >  8:58834/5‑149 (MQ=255)
                                                                                                                             aGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGtt                    >  8:144038/1‑149 (MQ=255)
                                                                                                                                   gTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGAt              >  1:281925/1‑149 (MQ=255)
                                                                                                                                               aCGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCAt  >  8:206716/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCAT  >  NC_002947/5701083‑5701374

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