Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,110 G→A 10.1% P1402S (CCG→TCG)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,1100GA10.1% 192.9 / 8.8 69P1402S (CCG→TCG) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (3/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AACAGCGCACCGCCCTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTGG  >  NC_002947/4782962‑4783252
                                                                                                                                                    |                                                                                                                                              
aaCAGCGCACCGCCCTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCgg                                                                                                                                                <  2:347296/149‑1 (MQ=255)
             ccTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCa                                                                                                                                     >  1:142367/1‑147 (MQ=255)
             ccTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCa                                                                                                                                     >  7:241428/1‑147 (MQ=255)
              cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  8:241428/149‑3 (MQ=255)
              cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  1:118524/149‑3 (MQ=255)
              cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  4:80979/149‑3 (MQ=255)
               tGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAActg                                                                                                                                 >  1:76230/1‑149 (MQ=255)
                 cTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAActgct                                                                                                                               <  1:194087/149‑1 (MQ=255)
                  tCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaa                                                                                                                                   <  5:140702/144‑2 (MQ=255)
                  tCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaa                                                                                                                                   >  6:140702/1‑143 (MQ=255)
                     cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        <  2:266956/136‑1 (MQ=255)
                     cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCtt                                                                                                                           <  6:41248/149‑1 (MQ=255)
                     cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCAGGGACCTCgg                                                                                                                                        >  1:266956/1‑136 (MQ=255)
                        aaGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  1:212354/149‑1 (MQ=255)
                         aGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCg                                                                                                                                         <  5:103037/131‑1 (MQ=255)
                         aGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCg                                                                                                                                         >  6:103037/1‑131 (MQ=255)
                             gCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             <  3:345982/139‑1 (MQ=255)
                             gCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             >  4:345982/1‑139 (MQ=255)
                                aTCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        >  7:208948/1‑125 (MQ=255)
                                aTCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        <  8:208948/125‑1 (MQ=255)
                                    aaCTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCt                                                                                                            >  2:345981/1‑149 (MQ=255)
                                      cTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTgg                                                                                                          >  5:145964/1‑149 (MQ=255)
                                      cTGGCCGAGGTACGTGAAGGGCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTACTGCTTGATGTGTTTGATGCTgg                                                                                                          <  7:140861/149‑1 (MQ=255)
                                         gCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGgccg                                                                                                       >  2:65979/1‑149 (MQ=255)
                                                     gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  8:105518/1‑147 (MQ=255)
                                                     gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  5:37179/1‑147 (MQ=255)
                                                     gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  7:210147/1‑147 (MQ=255)
                                                      aaGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAg                                                                                          >  6:187203/1‑149 (MQ=255)
                                                       aGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGc                                                                                         <  6:29809/149‑1 (MQ=255)
                                                          tCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgt                                                                                      <  5:341558/149‑1 (MQ=255)
                                                          tCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgt                                                                                      <  5:262369/149‑1 (MQ=255)
                                                           cACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        >  6:192213/1‑114 (MQ=255)
                                                           cACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  5:192213/114‑1 (MQ=255)
                                                            aCCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgtgt                                                                                    <  6:37179/149‑1 (MQ=255)
                                                              ccGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGc                                                                                                 >  8:305359/1‑134 (MQ=255)
                                                              ccGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGc                                                                                                 <  7:305359/134‑1 (MQ=255)
                                                                 cgGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGtt                                                                               <  2:76230/149‑1 (MQ=255)
                                                                    cTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCa                                                                            >  5:77301/1‑149 (MQ=255)
                                                                        ggCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCgg                                                                        <  5:37213/149‑1 (MQ=255)
                                                                              gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCaa                                                                                                                                    >  5:178743/1‑83 (MQ=255)
                                                                              gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCaa                                                                                                                                    <  6:178743/83‑1 (MQ=255)
                                                                              gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGaa                                                                     >  4:194323/1‑146 (MQ=255)
                                                                                  ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCAcc                                                                          <  1:276291/137‑1 (MQ=255)
                                                                                  ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCAcc                                                                          >  2:276291/1‑137 (MQ=255)
                                                                                  ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              >  5:42899/1‑149 (MQ=255)
                                                                                  ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              >  3:97154/1‑149 (MQ=255)
                                                                                  ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              <  3:194323/149‑1 (MQ=255)
                                                                                   cAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTgg                                                             <  1:217424/149‑1 (MQ=255)
                                                                                              cGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             >  7:79984/1‑74 (MQ=255)
                                                                                              cGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             <  8:79984/74‑1 (MQ=255)
                                                                                                  cgcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACgg                                              <  8:210147/149‑1 (MQ=255)
                                                                                                   gcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        >  2:355164/1‑74 (MQ=255)
                                                                                                   gcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  1:355164/74‑1 (MQ=255)
                                                                                                     ggCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGc                                           <  1:345981/149‑1 (MQ=255)
                                                                                                           gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCt                                       >  5:214918/1‑147 (MQ=255)
                                                                                                           gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCa                                     <  6:77301/149‑1 (MQ=255)
                                                                                                           gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCa                                     <  1:65979/149‑1 (MQ=255)
                                                                                                                ggTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCCCCGGGAACAGAGTGGTGAACCAGCCAACGGTGCGGCGCgggtcg                                >  4:330408/1‑149 (MQ=255)
                                                                                                                    aCGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCaggtcgagg                             <  5:187203/148‑1 (MQ=255)
                                                                                                                        agcaCGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGgtcgtc                        <  6:214918/149‑1 (MQ=255)
                                                                                                                           aCGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGaa                     <  1:338595/149‑1 (MQ=255)
                                                                                                                             gCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGaa                     >  8:5241/1‑147 (MQ=255)
                                                                                                                             gCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACa                   <  4:97154/149‑1 (MQ=255)
                                                                                                                                       aTACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTcgcg         <  6:42899/149‑1 (MQ=255)
                                                                                                                                         aCCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGc       >  3:86826/1‑149 (MQ=255)
                                                                                                                                           cTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCAACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCg     >  4:129491/1‑149 (MQ=255)
                                                                                                                                            tttGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCaggtcgag                              <  4:285736/123‑1 (MQ=255)
                                                                                                                                              tGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  >  6:210715/1‑149 (MQ=255)
                                                                                                                                              tGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  >  6:215840/1‑149 (MQ=255)
                                                                                                                                               ggCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  <  7:105518/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AACAGCGCACCGCCCTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTGG  >  NC_002947/4782962‑4783252

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