Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,528 C→G 6.8% V66L (GTG→CTG)  pde ← 3',5'‑cyclic‑nucleotide phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,590,5280CG 6.8% 180.8 / 2.6 59V66L (GTG→CTG) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base G (2/2);  total (28/31)
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.01e-01

CCCGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATC  >  NC_002947/5590381‑5590675
                                                                                                                                                    |                                                                                                                                                   
cccGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACg                                                                                                                                                    >  1:217241/1‑149 (MQ=255)
    ccACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTgg                                                                                                                                                >  3:43281/1‑149 (MQ=255)
     cACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGc                                                                                                                                               <  7:60202/149‑1 (MQ=255)
          gcCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCACAGGATGCCAGCGTGGCATcc                                                                                                                                           >  1:266964/1‑148 (MQ=255)
           ccGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCAGCGTGGCATCCta                                                                                                                                         <  1:14888/149‑2 (MQ=255)
            cGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTgg                                                                                                                                        <  7:143028/149‑1 (MQ=255)
              atatCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCACAGGATGCCAGCGTGGCATCCtatga                                                                                                                                      <  2:266964/149‑5 (MQ=255)
                 tCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACa                                                                                                                                    <  3:280032/148‑1 (MQ=255)
                   aGTCACCGCCTGCACCAGTTCAGCCGCCACCTGAGCCATCACCTCAGCCTCGT‑CATGGTTGCCCGGCAGACAGCCAGAGGGTACGGCAAACGATGCGGTCAGCTCCCGAAACGCTTCATAGGAAGCCACGCTGGCATCCTGGGACAgct                                                                                                                                 >  5:80592/1‑147 (MQ=255)
                     tCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTGCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTc                                                                                                                                >  1:89325/1‑149 (MQ=255)
                       aCCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGcc                                                                                                                             >  2:329916/1‑149 (MQ=255)
                        ccGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGccc                                                                                                                            <  4:43281/149‑1 (MQ=255)
                          gCCTGCACCAGTTCAGGCGCCACCTGCGCCATCAGCTGAGCCTCTT‑CATGGTTGGCCGGCAGCCAGCGAGTGGGTACGGCAAACTGTGCGGTCAGCTCCCGAAAGGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGt                                                                                                                          <  6:80592/149‑1 (MQ=255)
                          gCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGt                                                                                                                          <  2:217241/149‑1 (MQ=255)
                           ccTGCACCAGTTCCGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTg                                                                                                                         <  8:116274/149‑1 (MQ=255)
                              gCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGca                                                                                                                       >  2:166039/1‑148 (MQ=255)
                                aCCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcag                                                                                                                    >  5:100252/1‑149 (MQ=255)
                                    gTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACagcagc                                                                                                                <  1:329916/149‑1 (MQ=255)
                                          gcgcCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAg                                                                                                             >  2:206412/1‑146 (MQ=255)
                                          gcgcCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAg                                                                                                             <  1:206412/146‑1 (MQ=255)
                                             ccACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCTTCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGAc                                                                                                        >  5:281009/1‑148 (MQ=255)
                                               aCCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGAc                                                                                                        >  5:72816/1‑146 (MQ=255)
                                               aCCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGAc                                                                                                        <  6:72816/146‑1 (MQ=255)
                                                    cgcCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCgg                                                                                                  <  2:89325/147‑1 (MQ=255)
                                                       cATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCTTCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCAgctgct                                                                                             <  6:281009/149‑1 (MQ=255)
                                                        aTCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGc                                                                                                    <  3:112924/141‑1 (MQ=255)
                                                        aTCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGc                                                                                                    >  4:112924/1‑141 (MQ=255)
                                                          cACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTcgc                                                                                          >  1:178139/1‑149 (MQ=255)
                                                          cACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTAGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTcgg                                                                                          >  1:308482/1‑148 (MQ=255)
                                                           aCCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTcgcg                                                                                         >  6:333196/1‑149 (MQ=255)
                                                             cTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGcgcc                                                                                       >  4:329921/1‑149 (MQ=255)
                                                                accctcct‑CATCGTTGCCCGGCAGCCAGCGAGTCGCAACGGCAAACGGTGCGCTCAGCGCCAGAAACGCTTCAAAGGATCCCACGCTGGCATCCAGGGACAGGTCGCCCGTCCACAGCAGCAAGACGACAAGCCCCTGCACGCGCCGCa                                                                                    >  4:295442/3‑149 (MQ=255)
                                                                 gCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCa                                                                                    >  5:178446/1‑148 (MQ=255)
                                                                  ccTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCgcgc                                                                                        <  7:105521/143‑1 (MQ=255)
                                                                  ccTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCgcgc                                                                                        >  8:105521/1‑143 (MQ=255)
                                                                       t‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGc                                                                             <  1:155807/149‑1 (MQ=255)
                                                                       t‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGc                                                                             <  5:272009/149‑1 (MQ=255)
                                                                              ttGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAacac                                                                       >  2:324119/1‑149 (MQ=255)
                                                                                  ccGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTgg                                                                   <  1:324119/149‑1 (MQ=255)
                                                                                                gTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGGCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCg                                                     <  7:144610/149‑1 (MQ=255)
                                                                                                 tGGGTACGGCAATCGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACATCAGCAAGTCGACAAGCGGCTGCTCGCGCCGGACCTGGGCCAACACGTGGCGCAGGCTGTCTCgg                                                    <  3:295442/149‑1 (MQ=255)
                                                                                                        ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGa                                             <  2:178139/149‑1 (MQ=255)
                                                                                                        ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGa                                             <  6:100252/149‑1 (MQ=255)
                                                                                                        ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTAGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGa                                             <  2:308482/149‑1 (MQ=255)
                                                                                                         gCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTgg                                                                               <  8:277056/114‑1 (MQ=255)
                                                                                                         gCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTgg                                                                               >  7:277056/1‑114 (MQ=255)
                                                                                                         gCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAg                                            >  6:85921/1‑149 (MQ=255)
                                                                                                          cAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAgg                                           <  2:275482/149‑1 (MQ=255)
                                                                                                          cAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAgg                                           >  1:76232/1‑149 (MQ=255)
                                                                                                                   cGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCAtg                                  <  5:85921/149‑1 (MQ=255)
                                                                                                                      tCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGCTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGTGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCAtgctg                               >  3:287258/1‑149 (MQ=255)
                                                                                                                             ccGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCg                        >  3:206328/1‑149 (MQ=255)
                                                                                                                             ccGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCg                        >  2:27711/1‑149 (MQ=255)
                                                                                                                                    gCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGc                 <  6:178446/149‑1 (MQ=255)
                                                                                                                                    gCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGc                 <  3:123753/149‑1 (MQ=255)
                                                                                                                                        cATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGaa              <  5:333196/148‑1 (MQ=255)
                                                                                                                                          caGGATGCCAGCGTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGAGTCGGCGAACAg           >  1:305073/2‑149 (MQ=255)
                                                                                                                                              aTGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTgg       <  3:329921/149‑1 (MQ=255)
                                                                                                                                                   aCGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGCCTGCTCGCGCCGCACCTGGGCCAAAACGTGGCGCAGGCTGTCTCGTCTCTTGAGGACCAGCATGCTTCCAGACGGGTCGTCGAACAGGTCGGCATc  >  6:174129/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGT‑CATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATC  >  NC_002947/5590381‑5590675

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