Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,298,032 G→C 7.9% E83D (GAG→GAC) ‡ pvdQ → acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,298,0320GC 7.9% 195.3 / 5.8 63E83D (GAG→GAC) ‡pvdQacyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGCTGGCGG  >  NC_002947/3297888‑3298180
                                                                                                                                                |                                                                                                                                                    
cAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGAGTCGGATATGCATAAGCCAAGAACAACCAATGAATGAACGCCAATAAATAGAAAACCATGAGCGGTAAgcgc                                                                                                                                                   >  8:240004/1‑148 (MQ=255)
     cccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACgg                                                                                                                                             <  3:258980/149‑2 (MQ=255)
     cccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACg                                                                                                                                              <  1:70236/148‑1 (MQ=255)
        aGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAt                                                                                                                                          <  8:215549/149‑1 (MQ=255)
        aGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACggtca                                                                                                                                          <  6:355930/149‑5 (MQ=255)
              gCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGgc                                                                                                                                    >  6:121436/1‑149 (MQ=255)
               cTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGgcg                                                                                                                                   >  3:332453/1‑149 (MQ=255)
                     ccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCaatg                                                                                                                             >  3:219651/1‑147 (MQ=255)
                     ccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAgg                                                                                                                             >  5:58569/1‑149 (MQ=255)
                      cGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAggg                                                                                                                            <  5:289582/149‑1 (MQ=255)
                          aCGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCgg                                                                                                                                     <  5:271764/136‑1 (MQ=255)
                          aCGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCgg                                                                                                                                     >  6:271764/1‑136 (MQ=255)
                           cGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAg                                                                                                                        <  3:345821/148‑1 (MQ=255)
                             gCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAcc                                                                                                                     >  7:202092/1‑149 (MQ=255)
                             gCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAcc                                                                                                                     >  5:4401/1‑149 (MQ=255)
                              cGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAccc                                                                                                                    <  6:157041/149‑1 (MQ=255)
                                tACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTg                                                                                                                  <  2:229868/149‑1 (MQ=255)
                                  ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTgg                                                                                                                 <  6:4401/148‑1 (MQ=255)
                                  ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTgg                                                                                                                 <  8:202092/148‑1 (MQ=255)
                                  ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                                <  7:360001/149‑1 (MQ=255)
                                  ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                                <  5:121436/149‑1 (MQ=255)
                                     cATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTTAGCGGTGAGCGCTCACGCTATTTCGGCTCCAAGGGCCAGACCCTGGagca                                                                                                             <  7:240004/149‑1 (MQ=255)
                                           gcggCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGa                                                                                                       <  5:258613/147‑1 (MQ=255)
                                            tGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAAGCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAATTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGa                                                                                                       >  8:350885/1‑148 (MQ=255)
                                            gggCCCAGGATGTGCGTGGCCTGGGCTATGGCGTCGGCTATTCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGGTGACCGTGAGCGGTGTGCGCTCAGGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGAc                                                                                                      <  6:126711/148‑1 (MQ=255)
                                               ccAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAttt                                                                                                                                        >  7:243576/1‑112 (MQ=255)
                                               ccAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAttt                                                                                                                                        <  8:243576/112‑1 (MQ=255)
                                               ccAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAAc                                                                                                   >  1:153159/1‑149 (MQ=255)
                                                cAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAAcc                                                                                                  >  5:279529/1‑149 (MQ=255)
                                                 aaGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCt                                                                                                 <  6:58569/149‑1 (MQ=255)
                                                  aGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAgcgc                                                                                                                                                   >  6:74385/1‑98 (MQ=255)
                                                  aGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAgcgc                                                                                                                                                   <  5:74385/98‑1 (MQ=255)
                                                    gATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGc                                                                                              <  4:332453/149‑1 (MQ=255)
                                                        aGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGc                                                                                          <  6:279529/149‑1 (MQ=255)
                                                               ccTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGt                                                                                   >  3:285621/1‑149 (MQ=255)
                                                                  gggCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttct                                                                                >  8:22034/1‑149 (MQ=255)
                                                                        tGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                            <  4:219651/147‑1 (MQ=255)
                                                                          gCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTgg                                                                        >  1:62343/1‑149 (MQ=255)
                                                                          gCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAATTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTgg                                                                        <  7:350885/149‑1 (MQ=255)
                                                                           cGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                                 <  1:200265/139‑1 (MQ=255)
                                                                           cGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                                 >  2:200265/1‑139 (MQ=255)
                                                                                gCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCaa                                                                   >  6:51773/1‑148 (MQ=255)
                                                                                  tATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                            >  6:212996/1‑137 (MQ=255)
                                                                                  tATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                            <  5:212996/137‑1 (MQ=255)
                                                                                       cTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCa                                                                                                                         <  4:347228/87‑1 (MQ=255)
                                                                                       cTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCa                                                                                                                         >  3:347228/1‑87 (MQ=255)
                                                                                       cTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAg                                                           >  3:57041/1‑149 (MQ=255)
                                                                                            cccAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgcggcg                                                      <  7:22034/149‑1 (MQ=255)
                                                                                              cAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGt                                                    >  1:263956/1‑149 (MQ=255)
                                                                                                  aCAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGAt                                                <  2:62343/149‑1 (MQ=255)
                                                                                                      ccTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCt                                            >  3:160861/1‑149 (MQ=255)
                                                                                                          tGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTc                                          <  4:285621/147‑1 (MQ=255)
                                                                                                           gCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTg                                       >  8:99521/1‑149 (MQ=255)
                                                                                                            ccTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTg                                       <  7:164833/148‑1 (MQ=255)
                                                                                                              tGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAg                                    >  8:244924/1‑149 (MQ=255)
                                                                                                                  cGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCa                                                                                            <  8:89097/89‑1 (MQ=255)
                                                                                                                  cGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCa                                                                                            >  7:89097/1‑89 (MQ=255)
                                                                                                                  cGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGt                                                    <  7:192104/129‑1 (MQ=255)
                                                                                                                  cGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGt                                                    >  8:192104/1‑129 (MQ=255)
                                                                                                                       aTGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGGTCTTCACCTGGCTCAATAGACCAGCGGCGGTCGATGCATTCCTGCATGCGCAAccc                           >  8:329887/1‑149 (MQ=255)
                                                                                                                        tGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCcg                          >  8:193406/1‑149 (MQ=255)
                                                                                                                        tGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCcg                          >  1:263015/1‑149 (MQ=255)
                                                                                                                               cTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGt                   >  2:207524/1‑149 (MQ=255)
                                                                                                                                tGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTg                  >  3:306244/1‑149 (MQ=255)
                                                                                                                                   gcGTGAGCGCTCACCGCTCACGGTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg               >  6:344584/2‑149 (MQ=255)
                                                                                                                                   gcGTGAGCGCTCACCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCAGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg               >  1:95794/2‑149 (MQ=255)
                                                                                                                                     gtgagcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAgg              >  2:306571/1‑148 (MQ=255)
                                                                                                                                     gtgagcgctcaccgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAgg              >  7:236120/1‑148 (MQ=255)
                                                                                                                                        agcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCa            >  1:136801/1‑147 (MQ=255)
                                                                                                                                                gcgcTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGCTGGCgg  >  2:125754/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGCTGGCGG  >  NC_002947/3297888‑3298180

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: ‑