Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,782,601 C→G 9.4% P1571P (CCG→CCC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,782,6010CG 9.4% 181.3 / 6.5 63P1571P (CCG→CCC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (3/3);  total (29/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

GCACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAA  >  NC_002947/4782453‑4782740
                                                                                                                                                    |                                                                                                                                           
gcgccATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAAccc                                                                                                                                             >  1:190922/4‑149 (MQ=255)
gcgccATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAAccc                                                                                                                                             >  5:23563/4‑149 (MQ=255)
  aCCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCgg                                                                                                                                           <  3:21061/149‑1 (MQ=255)
     aTGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTc                                                                                                                                        <  1:11389/149‑1 (MQ=255)
        cgcTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                      <  6:23563/148‑1 (MQ=255)
        cgcTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                      <  2:190922/148‑1 (MQ=255)
         gcTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAg                                                                                                                                     >  5:113850/1‑148 (MQ=255)
          cTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGACCGAACCCGGTTCGAGc                                                                                                                                    >  7:281352/1‑148 (MQ=255)
               ccAGGTCACAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtg                                                                                                                              >  6:105500/1‑149 (MQ=255)
                    tCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                 >  6:252694/1‑141 (MQ=255)
                    tCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                 <  5:252694/141‑1 (MQ=255)
                        aGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcag                                                                                                                      >  1:209216/1‑148 (MQ=255)
                         gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcagca                                                                                                                    <  5:105500/149‑1 (MQ=255)
                         gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcag                                                                                                                      >  5:91490/1‑147 (MQ=255)
                                       gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                      <  3:165667/117‑1 (MQ=255)
                                       gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                      >  4:165667/1‑117 (MQ=255)
                                       gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATc                                                                                                      <  2:201030/149‑1 (MQ=255)
                                       gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATc                                                                                                      <  7:205375/149‑1 (MQ=255)
                                        cagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCg                                                                                                                                       <  8:252332/115‑1 (MQ=255)
                                        cagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCg                                                                                                                                       >  7:252332/1‑115 (MQ=255)
                                          gcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGc                                                                                                   <  6:113850/149‑1 (MQ=255)
                                              gTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGATGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGAc                                                                                                >  8:274723/1‑148 (MQ=255)
                                                        gcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGAt                                                                                     <  7:268605/149‑1 (MQ=255)
                                                              tccGGTTTCCATGCTTCGCTGAACGCGCCGTCATCTATGGCGTGGTCGACGATGTAGCAATACTGTATCAGATAGATGCCCGATCCCGTTTCGAGCAAGGTGTGCATCTGCAGCCCTTCCTGCATCTGCGACAACGGATAGATGTCtt                                                                                <  8:281352/147‑1 (MQ=255)
                                                              aCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTc                                                                               <  2:209216/149‑1 (MQ=255)
                                                                cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                             >  6:310297/1‑149 (MQ=255)
                                                                            ttCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCgg                                                                 >  8:123454/1‑149 (MQ=255)
                                                                               gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTAGATATCGCGCGCGGCCa                                                              >  2:191749/1‑149 (MQ=255)
                                                                                 tGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCAcc                                                            >  5:31441/1‑149 (MQ=255)
                                                                                            aCATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                 <  5:75034/149‑1 (MQ=255)
                                                                                                 aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                            <  5:310297/149‑1 (MQ=255)
                                                                                                 aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTAGATATCGCGCGCGGCCACCGGCAGGCGATCCAGt                                             <  1:191749/148‑1 (MQ=255)
                                                                                                  aTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTtg                                           >  6:41930/1‑149 (MQ=255)
                                                                                                       gTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGAAGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTgctgct                                      >  3:77177/1‑149 (MQ=255)
                                                                                                          gTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTgg                                    <  6:91490/148‑1 (MQ=255)
                                                                                                          gTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTgg                                    <  2:147017/148‑1 (MQ=255)
                                                                                                             gACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCTCGCGGCCACCGGCAGCCGATCCAGATGCTGCTGGCTCa                                >  6:147637/1‑149 (MQ=255)
                                                                                                               cgcgGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCgg                                                                 >  1:242727/1‑114 (MQ=255)
                                                                                                               cgcgGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCgg                                                                 <  2:242727/114‑1 (MQ=255)
                                                                                                               cgcgGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCgg                                                                 <  2:246537/114‑1 (MQ=255)
                                                                                                               cgcgGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCgg                                                                 >  1:246537/1‑114 (MQ=255)
                                                                                                                       cAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAg                      >  5:54421/1‑149 (MQ=255)
                                                                                                                        aaTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCa                                               >  5:305534/1‑123 (MQ=255)
                                                                                                                        aaTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCa                                               <  6:305534/123‑1 (MQ=255)
                                                                                                                        aaTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGt                     <  6:54421/149‑1 (MQ=255)
                                                                                                                         aTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTg                    <  5:41930/149‑1 (MQ=255)
                                                                                                                         aTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTGCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTg                    <  7:274723/149‑1 (MQ=255)
                                                                                                                           aCTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                   >  5:276372/1‑148 (MQ=255)
                                                                                                                             tGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTggg                  >  2:224650/1‑147 (MQ=255)
                                                                                                                               cATCAGATAGATGCCCGAACCCGGTTCGAGCTAGGTGTGCATCAGCAGGCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACGGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGaa                <  5:147637/147‑1 (MQ=255)
                                                                                                                                aTCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCGCAGCCCGGCCAGTGGGAAGTc             >  6:312194/1‑149 (MQ=255)
                                                                                                                                 tCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCgg                                                          <  7:236820/103‑1 (MQ=255)
                                                                                                                                 tCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCgg                                                          >  8:236820/1‑103 (MQ=255)
                                                                                                                                  cAGAAAGATCCCCGACCCCGGTTCGAGCAACGTGCGCATCAGCAGCCCCTCCTGCATCGGCGACAACGCATAGATGCCTTCCATATCGCGCGCGGCCACCGGCAGGCGATCCACTTGCCGCTGGCTCAGCCCGG‑CAGTGGGAAGTCGCt          <  4:77177/149‑1 (MQ=255)
                                                                                                                                   agatagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGc           <  7:195090/148‑1 (MQ=255)
                                                                                                                                     atagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTg        <  7:123454/149‑1 (MQ=255)
                                                                                                                                        gatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCGCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGc      <  5:275541/148‑1 (MQ=255)
                                                                                                                                        gatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCg     >  2:137952/1‑149 (MQ=255)
                                                                                                                                         ttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGaa                >  1:114628/2‑137 (MQ=255)
                                                                                                                                         ttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGaa                <  2:114628/136‑1 (MQ=255)
                                                                                                                                           gctcgcaCCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCCGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCACAGCCCGGCCAGTGGGAAGTCGCtt         >  5:103562/7‑142 (MQ=255)
                                                                                                                                           gctcGAACGGGGTTCTGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACATCGGATAGATGTCTTCGATATCGAGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGGCCGGCCAGTGGGAAGTCGCtt         <  6:103562/139‑1 (MQ=255)
                                                                                                                                           gCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTaa  >  5:25866/1‑149 (MQ=255)
                                                                                                                                            ctcGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAg               >  2:271683/3‑135 (MQ=255)
                                                                                                                                            ctcGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAg               <  1:271683/133‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GCACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAA  >  NC_002947/4782453‑4782740

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