Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,790,888 G→A 8.3% N984N (AAC→AAT)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,790,8880GA 8.3% 172.4 / 4.2 60N984N (AAC→AAT) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TTGTCGTCCAAGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGCCACCAGGTAG  >  NC_002947/4790742‑4791035
                                                                                                                                                  |                                                                                                                                                   
ttGTCGTCCAAGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCAtt                                                                                                                                                   >  7:126312/1‑149 (MQ=255)
     gtcCAAGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCa                                                                                                                                               <  2:164914/148‑1 (MQ=255)
      tcCAAGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCATTgactgc                                                                                                                                             >  1:1056/1‑144 (MQ=255)
          aGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAAt                                                                                                                                         <  7:46915/149‑1 (MQ=255)
              gACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCa                                                                                                                                     <  8:105564/149‑1 (MQ=255)
              gACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCa                                                                                                                                     <  8:220036/149‑1 (MQ=255)
                ccTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCATTGACTGCCAACGGCaag                                                                                                                                   <  2:1056/149‑3 (MQ=255)
                ccTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCATTGACTGCCAACGGCaag                                                                                                                                   <  8:126312/149‑3 (MQ=255)
                 cTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGc                                                                                                                                  >  5:309991/1‑149 (MQ=255)
                 cTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCATTGACTGCCAACGGCaagg                                                                                                                                  <  5:211025/149‑4 (MQ=255)
                      gACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTgg                                                                                                                                                 <  4:92106/129‑1 (MQ=255)
                      gACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTgg                                                                                                                                                 >  3:92106/1‑129 (MQ=255)
                           cccAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCa                                                                                                                                               >  5:41880/1‑126 (MQ=255)
                           cccAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCa                                                                                                                                               <  6:41880/126‑1 (MQ=255)
                           cccAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCt                                                                                                                                 <  7:16904/140‑1 (MQ=255)
                           cccAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCt                                                                                                                                 >  8:16904/1‑140 (MQ=255)
                                        gCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCt                                                                                                           >  7:302504/1‑149 (MQ=255)
                                         ccAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCt                                                                                                           <  5:219882/148‑1 (MQ=255)
                                         ccAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTg                                                                                                          >  4:266037/1‑149 (MQ=255)
                                           aaaGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTgg                                                                                                         >  3:178750/1‑148 (MQ=255)
                                                ttGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGTCTGTCGAGGAACAACACGCGGGCTGGCACCAt                                                                                                   <  1:213740/149‑1 (MQ=255)
                                                    ccAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCCGTCGAGGAACAACACGCGGGCTGGCACCATGTgc                                                                                               >  6:192327/1‑149 (MQ=255)
                                                           gATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGc                                                                                         <  6:309991/148‑1 (MQ=255)
                                                             tCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACg                                                                                      >  5:242165/1‑149 (MQ=255)
                                                               tGCAGTGGCCTCTGAGCCGGTGTGTGGCTGGCCACCGGGTCTTGTTCGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACCCGGGCTGGCACCAAGTGCGCCAGCACGtc                                                                                    >  2:121064/1‑148 (MQ=255)
                                                                 cAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCAc                                                                                       <  4:63583/144‑1 (MQ=255)
                                                                 cAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGt                                                                                     <  3:185621/146‑1 (MQ=255)
                                                                 cAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTg                                                                                    <  6:242165/147‑1 (MQ=255)
                                                                 cAGTGGCGTCTGAGCCGGTGTGTAGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGc                                                                                   >  4:6650/1‑148 (MQ=255)
                                                                      gCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACtt                                                                              <  4:177242/148‑1 (MQ=255)
                                                                            gAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCt                                                                       >  4:191237/1‑149 (MQ=255)
                                                                               ccGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTtgc                                                                    >  7:183485/1‑149 (MQ=255)
                                                                                   ccggtgGCCAGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCa                                                                >  7:3891/4‑149 (MQ=255)
                                                                                     tgtgGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACg                                                              <  3:266037/149‑1 (MQ=255)
                                                                                       tggctggcCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCtt                                                                      >  6:294412/1‑139 (MQ=255)
                                                                                       tggctggcCACCGGGTGTTGTTGGTCCGGCACCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCtt                                                                      <  5:294412/139‑1 (MQ=255)
                                                                                               cACCGGGTGTTGTTGGTCCGGCAGCGGCAGTCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCtgtg                                                    >  1:25018/1‑149 (MQ=255)
                                                                                                 ccGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCCGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTgcc                                                  <  5:192327/149‑1 (MQ=255)
                                                                                                          gttgGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCa                                         >  8:165986/1‑149 (MQ=255)
                                                                                                           ttgGTCCGGCAGCGGCAATCTGGTCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGATTGCTGCACGGCCTCCTGTGCCGCATCGGCAt                                        >  8:119129/1‑149 (MQ=255)
                                                                                                           ttgGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCAt                                        <  8:302504/149‑1 (MQ=255)
                                                                                                           ttgGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCAt                                        >  4:33613/1‑149 (MQ=255)
                                                                                                              gTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCAc                                     <  6:148355/149‑1 (MQ=255)
                                                                                                               tCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACg                                    >  6:92461/1‑149 (MQ=255)
                                                                                                                 cGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCt                                  >  4:47071/1‑149 (MQ=255)
                                                                                                                      gCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCgg                             <  7:165986/149‑1 (MQ=255)
                                                                                                                         gCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCGCGGCATCCCGCTCGCGGGc                           >  5:64647/1‑148 (MQ=255)
                                                                                                                            aTCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGAcag                        >  3:47505/1‑148 (MQ=255)
                                                                                                                               cggCAAGGTTGACCTGGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCagcagc                    >  7:234055/2‑149 (MQ=255)
                                                                                                                                 gCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCcg                  >  1:188254/1‑149 (MQ=255)
                                                                                                                                   cAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCcgcg                <  3:47071/149‑1 (MQ=255)
                                                                                                                                     ggTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCgcggc              <  3:6650/149‑1 (MQ=255)
                                                                                                                                         aaCCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCGTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGcca          >  8:165018/1‑149 (MQ=255)
                                                                                                                                         aaCCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGcca          >  2:47083/1‑149 (MQ=255)
                                                                                                                                            cTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCcgcg                <  1:295886/140‑1 (MQ=255)
                                                                                                                                            cTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCcgcg                >  2:295886/1‑140 (MQ=255)
                                                                                                                                            cTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGccacca       <  8:183485/149‑1 (MQ=255)
                                                                                                                                                ccGTTGGCAGTCAATGGCAGGCTGTCGAGGAACATCACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGCCACCAGGTa   <  1:121064/149‑1 (MQ=255)
                                                                                                                                                 cGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGCCACCAGGTAg  <  3:191237/149‑1 (MQ=255)
                                                                                                                                                 cGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGCCACCAGGTAg  <  5:92461/149‑1 (MQ=255)
                                                                                                                                                  gTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTgc                                                   >  3:125097/1‑99 (MQ=255)
                                                                                                                                                  gTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTgc                                                   <  4:125097/99‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGTCGTCCAAGCCGACCTGCGGACAGCCCAGCACTTCCTGCCAAAGGTTGGCCAGGTCGATCTGCAGTGGCGTCTGAGCCGGTGTGTGGCTGGCCACCGGGTGTTGTTGGTCCGGCAGCGGCAATCTGGCCAGGTCAACCTTGCCGTTGGCAGTCAATGGCAGGCTGTCGAGGAACAACACGCGGGCTGGCACCATGTGCGCCTGCACGTGCAACTTCAGCGCTTGCTGCACGGCCTCCTGTGCCGCCTCGGCATCACGCTCGCGGGCCAGCAGCCGCGGCGCCACCAGGTAG  >  NC_002947/4790742‑4791035

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