Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,299,392 C→T 6.2% Q537* (CAG→TAG)  pvdQ → acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,299,3920CT 6.2% 193.0 / 2.0 64Q537* (CAG→TAG) pvdQacyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/2);  total (30/34)
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

GACCCGGCCTGGATGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGT  >  NC_002947/3299248‑3299538
                                                                                                                                                 |                                                                                                                                                  
gACCCGGCCTGGATGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAgcg                                                                                                                                                >  2:33163/1‑149 (MQ=255)
 aCCCGGCCTGGATGGCCAACCCGGCGCAGCCGCTGAGCCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAgcg                                                                                                                                                >  1:61006/1‑149 (MQ=255)
   ccGGCCTGGATGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCTACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCAt                                                                                                                                             <  3:33998/149‑1 (MQ=255)
      gCCTGGATGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGt                                                                                                                                          <  3:285645/149‑1 (MQ=255)
             tGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGAt                                                                                                                                   <  1:353074/149‑1 (MQ=255)
                     ccGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACTGATGAAGAGgt                                                                                                                           <  8:129874/149‑1 (MQ=255)
                     ccGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGgt                                                                                                                           >  1:210756/1‑149 (MQ=255)
                     ccGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGgt                                                                                                                           >  7:215562/1‑149 (MQ=255)
                          gcAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTAcc                                                                                                                      >  2:102993/1‑149 (MQ=255)
                             gCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTgg                                                                                                                   <  8:212578/149‑1 (MQ=255)
                              ccGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGc                                                                                                                  >  1:94836/1‑149 (MQ=255)
                                       cGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTAc                                                                                                                       <  4:284955/136‑1 (MQ=255)
                                       cGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTAc                                                                                                                       >  3:284955/1‑136 (MQ=255)
                                             cTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGc                                                                                                    <  1:102993/149‑1 (MQ=255)
                                                gCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGGCAGCCTGGTGTTGc                                                                                                    >  8:303774/1‑146 (MQ=255)
                                                    cTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACct                                                                                             <  6:206897/149‑1 (MQ=255)
                                                     tGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACctg                                                                                            >  3:156428/1‑149 (MQ=255)
                                                         cAGCCGCAATGACAAACCACTGGGCATGCGAGGCCGCTTTGCCTTGCAACGCCTGCAGGCCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACctgctg                                                                                         >  8:290462/1‑148 (MQ=255)
                                                          aGCCGCAATGACCAACCTCTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCtgctgc                                                                                        <  2:210852/148‑1 (MQ=255)
                                                            ccGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAAt                                                                                     <  1:9058/149‑1 (MQ=255)
                                                             cGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGgt                                                                                                                           <  6:251275/110‑1 (MQ=255)
                                                             cGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGgt                                                                                                                           >  5:251275/1‑110 (MQ=255)
                                                               cAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGt                                                                                  <  2:276226/149‑1 (MQ=255)
                                                               cAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGt                                                                                  <  2:376140/149‑1 (MQ=255)
                                                                       aaccaCGGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGgc                                                                          <  7:39794/149‑1 (MQ=255)
                                                                        accaCTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGgcg                                                                         >  7:295761/1‑149 (MQ=255)
                                                                          caCTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGgcgc                                                                        >  3:101885/1‑148 (MQ=255)
                                                                             tGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAgc                                                                    >  4:28946/1‑149 (MQ=255)
                                                                                  aTGCGTGGCCGCTTTGCCATGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGAATGGTCACCGATGAAGAGGTGAACCTGCCCAGCCTAGTGTTGCCAGACCTGCTCCAATGATTAAAGGGAGCCAgag                                                                   >  3:300460/1‑143 (MQ=255)
                                                                                     cGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGAtg                                                             <  2:94836/148‑1 (MQ=255)
                                                                                           cGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGc                                                      <  5:265140/149‑1 (MQ=255)
                                                                                            gCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGcc                                                     >  4:288283/1‑149 (MQ=255)
                                                                                                tGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCgt                                                   <  7:303774/147‑1 (MQ=255)
                                                                                                 gCCTTGCAACGCCTGCATGGCATAGCCCGGTTATGTGTAGTTGAATTGCAGCGCATGGTCACCGATGAATAGGTGTACCTGGCCAGCCTGGTGTTGGCAGACCTGGTGCAATGGTGCAAGGTCGCCAGCGCTGATGTGCAGGCCgtgt                                                 <  8:295761/148‑1 (MQ=255)
                                                                                                 gCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCgtgt                                                 <  2:61006/148‑1 (MQ=255)
                                                                                                 gCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCgtgt                                                 <  4:101885/148‑1 (MQ=255)
                                                                                                  ccTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGTAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGtgtg                                                >  4:99571/1‑148 (MQ=255)
                                                                                                  ccTTGCAACGCATGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGTAGCGCATGGTCACAGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGATGAAATGGTGCAAGGGCTCCAGCGCATCTGTGCATGCCGtgtg                                                >  4:99573/1‑148 (MQ=255)
                                                                                                   cTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGtgtg                                                >  3:364696/1‑147 (MQ=255)
                                                                                                    ttGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCAAGCGCAGATGTGCAGGCCGTGTgcag                                             >  7:261639/1‑149 (MQ=255)
                                                                                                      gCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTTGCCAGCCTGgt                                                                                                         >  4:311138/1‑87 (MQ=255)
                                                                                                      gCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGgt                                                                                                         <  3:311138/87‑1 (MQ=255)
                                                                                                       cAACGCCTGCAGGCCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTgcagca                                           <  7:290462/148‑1 (MQ=255)
                                                                                                           gCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTg                                      <  4:156428/149‑1 (MQ=255)
                                                                                                              tGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCAGACCTGCTGCAATGGTGCAAGGGTGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGgccg                                  >  2:242502/1‑149 (MQ=255)
                                                                                                               gCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGgccg                                  >  1:267845/1‑149 (MQ=255)
                                                                                                                  tggCAAAGGCCGGTTAGGTGTAGATGAATTGCTGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCt                               <  8:261639/148‑1 (MQ=255)
                                                                                                                  gggCAAAGCCCTGGTAGGTGACGGTGAATTTCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCt                               <  8:347101/149‑1 (MQ=255)
                                                                                                                   ggCATAGCGCGGTTAGGTGTAGATGAATTGCAGCTCATGGTCACCGATGAAGAGGTGTACCTGGCGAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTTTGAATCAGGCTGGCCGACt                               <  4:300460/148‑1 (MQ=255)
                                                                                                                   ggCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCagcagc                                         <  8:153564/138‑1 (MQ=255)
                                                                                                                   ggCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCagcagc                                         >  7:153564/1‑138 (MQ=255)
                                                                                                                    gCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTgg                             <  8:215562/149‑1 (MQ=255)
                                                                                                                     cAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTgg                             >  8:251479/1‑148 (MQ=255)
                                                                                                                          cccGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGc                       <  7:310057/149‑1 (MQ=255)
                                                                                                                                   gtgtAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGAc              >  1:220454/1‑149 (MQ=255)
                                                                                                                                        gATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCATGGCCGACCTCGa         >  7:191703/1‑149 (MQ=255)
                                                                                                                                        gATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGa         >  3:285387/1‑149 (MQ=255)
                                                                                                                                        gATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGa         >  7:149534/1‑149 (MQ=255)
                                                                                                                                          tGAATTGTAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCATACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTTGAACGGCAAGGCCGACCTCGACa       <  3:99573/149‑1 (MQ=255)
                                                                                                                                          tGAATTGTAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACa       <  3:99571/149‑1 (MQ=255)
                                                                                                                                          tGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACa       <  7:251479/149‑1 (MQ=255)
                                                                                                                                             aTTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTg    <  3:288283/149‑1 (MQ=255)
                                                                                                                                              ttGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTgg   >  4:262639/1‑149 (MQ=255)
                                                                                                                                                gCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGt  <  8:337985/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GACCCGGCCTGGATGGCCAACCCGGCGCAGCCGCTGACCGGCTACTCGCCACTGGTCAGCCGCAATGACCAACCACTGGGCATGCGTGGCCGCTTTGCCTTGCAACGCCTGCAGGGCAAAGCCCGGTTAGGTGTAGATGAATTGCAGCGCATGGTCACCGATGAAGAGGTGTACCTGGCCAGCCTGGTGTTGCCAGACCTGCTGCAATGGTGCAAGGGCGCCAGCGCAGATGTGCAGGCCGTGTGCAGCAGCCTGGCCGCCTGGAACGGCAAGGCCGACCTCGACAGTGGT  >  NC_002947/3299248‑3299538

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: