Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,700,369 T→C 8.3% S283P (TCG→CCG)  phaA → poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,700,3690TC 8.3% 168.0 / 4.5 60S283P (TCG→CCG) phaApoly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (2/3);  total (27/33)
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.00e+00

AGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGA  >  NC_002947/5700222‑5700517
                                                                                                                                                   |                                                                                                                                                    
aGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTc                                                                                                                                                     <  8:145448/149‑1 (MQ=255)
aGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTc                                                                                                                                                     <  8:41983/149‑1 (MQ=255)
  cccGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGa                                                                                                                                                   <  4:145476/149‑1 (MQ=255)
  cccGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCgg                                                                                                                                                   >  6:127893/1‑148 (MQ=255)
  cccGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGAAGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCgg                                                                                                                                                   >  4:173851/1‑148 (MQ=255)
      gAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTAc                                                                                                                                               <  3:229347/149‑1 (MQ=255)
             gCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgcagc                                                                                                                                        <  8:57597/149‑1 (MQ=255)
             gCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagc                                                                                                                                        <  2:93001/149‑1 (MQ=255)
                  gcgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAgg                                                                                                                                    <  5:42667/148‑1 (MQ=255)
                  gcgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagcacc                                                                                                                                     <  5:127893/147‑3 (MQ=255)
                  gcgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGACGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagcacc                                                                                                                                     <  3:173851/147‑3 (MQ=255)
                    gcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGAcc                                                                                                                                 <  4:147121/149‑1 (MQ=255)
                        ttCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGaa                                                                                                                             >  1:99248/1‑149 (MQ=255)
                          ctgcctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAAc                                                                                                                            <  6:139912/148‑1 (MQ=255)
                            gcctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgcag                                                                                                                                         <  7:278898/133‑1 (MQ=255)
                            gcctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgcag                                                                                                                                         >  8:278898/1‑133 (MQ=255)
                             cctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGc                                                                                                                        >  7:93729/1‑149 (MQ=255)
                                  gcTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCAGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGgc                                                                                                                   >  7:153014/1‑149 (MQ=255)
                                      gcagcagcaGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGgcgc                                                                                                                 <  2:322218/147‑1 (MQ=255)
                                        agcagcagACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGc                                                                                                             >  3:243314/1‑149 (MQ=255)
                                        agcagcagACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTAACGGCAGCAAGGACCTGAAAATGCTCTGAGCCTGc                                                                                                             >  6:70645/1‑149 (MQ=255)
                                          cagcagACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTc                                                                                                           <  7:36256/149‑1 (MQ=255)
                                            gcagACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCAGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCg                                                                                                         <  8:153014/149‑1 (MQ=255)
                                                                      aCCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCt                                                                                                               >  7:50358/1‑117 (MQ=255)
                                                                      aCCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCt                                                                                                               <  8:50358/117‑1 (MQ=255)
                                                                               aGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTAt                                                                      >  1:200285/1‑149 (MQ=255)
                                                                                 gCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgc                                                                    >  2:92547/1‑149 (MQ=255)
                                                                                   ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAAc                                                                                                                            >  7:174299/1‑91 (MQ=255)
                                                                                   ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAAc                                                                                                                            <  8:174299/91‑1 (MQ=255)
                                                                                   ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgccg                                                                  >  3:64251/1‑149 (MQ=255)
                                                                                     aGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgccgc                                                                 >  7:13840/1‑148 (MQ=255)
                                                                                     aGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgccgc                                                                 >  1:132651/1‑148 (MQ=255)
                                                                                           aaTGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGc                                                          <  7:34153/149‑1 (MQ=255)
                                                                                            aTGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCg                                                         <  3:131890/149‑1 (MQ=255)
                                                                                             tGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCga                                                        <  2:63790/149‑1 (MQ=255)
                                                                                             tGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCga                                                        <  5:266806/149‑1 (MQ=255)
                                                                                             tGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCga                                                        <  2:200285/149‑1 (MQ=255)
                                                                                              ggggCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGc                                                          >  6:128197/1‑146 (MQ=255)
                                                                                              ggggCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGCCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGc                                                          >  2:177018/1‑146 (MQ=255)
                                                                                                ggCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtg                                                                                  >  2:65018/1‑120 (MQ=255)
                                                                                                ggCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtg                                                                                  <  1:65018/120‑1 (MQ=255)
                                                                                                   cTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACaa                                                  <  1:291274/149‑1 (MQ=255)
                                                                                                      tcctccTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGTCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGATAACaat                                                 <  5:70645/147‑2 (MQ=255)
                                                                                                      tCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGt                                               <  8:13840/149‑1 (MQ=255)
                                                                                                           cTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACg                                          >  2:111646/1‑149 (MQ=255)
                                                                                                           cTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACg                                          >  1:243593/1‑149 (MQ=255)
                                                                                                              cATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGccc                                       <  4:64251/149‑1 (MQ=255)
                                                                                                               aTCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCt                                      <  1:111646/149‑1 (MQ=255)
                                                                                                                 cGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGa                                    >  2:93989/1‑149 (MQ=255)
                                                                                                                  gATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGAc                                   <  2:99248/149‑1 (MQ=255)
                                                                                                                     gcgcTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCct                                >  4:91912/1‑149 (MQ=255)
                                                                                                                       gcTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCctgc                              >  4:18336/1‑149 (MQ=255)
                                                                                                                         tGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCctgctg                            >  7:70521/1‑149 (MQ=255)
                                                                                                                          gAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACTCCCTGACCCTGCTgg                           >  5:100425/1‑149 (MQ=255)
                                                                                                                                     tCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTg                                     <  7:169431/128‑1 (MQ=255)
                                                                                                                                     tCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTg                                     >  8:169431/1‑128 (MQ=255)
                                                                                                                                     tCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTg                <  4:243314/149‑1 (MQ=255)
                                                                                                                                          gcgGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGacac           <  3:91912/149‑1 (MQ=255)
                                                                                                                                            ggTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGAcacca         >  5:76449/1‑149 (MQ=255)
                                                                                                                                                cTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCAt     <  3:18336/149‑1 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCa                        >  4:81646/1‑127 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCa                        <  3:81646/127‑1 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGa  <  6:76449/149‑1 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGa  >  1:49721/1‑149 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGa  >  1:45065/1‑149 (MQ=255)
                                                                                                                                                   tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGa  <  8:70521/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGGA  >  NC_002947/5700222‑5700517

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