Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,062,321 T→C 8.3% Q673R (CAG→CGG)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,062,3210TC 8.3% 166.9 / 5.0 60Q673R (CAG→CGG) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base C (2/3);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGTCA  >  NC_002947/4062174‑4062467
                                                                                                                                                   |                                                                                                                                                  
gaggttgagCGCATCATCGCGGCCGTGGAGGCGGTAATGTCTGAGTTCATACCCGCCTTTACCTTGCGCGATCCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTg                                                                                                                                                   <  6:309428/143‑1 (MQ=255)
gATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTg                                                                                                                                                   <  7:19999/149‑1 (MQ=255)
gATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTg                                                                                                                                                   >  3:239870/1‑149 (MQ=255)
gATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTg                                                                                                                                                   >  3:204464/1‑149 (MQ=255)
 aTGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTgg                                                                                                                                                  <  2:183696/149‑1 (MQ=255)
      gagCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCca                                                                                                                                              <  3:227548/148‑1 (MQ=255)
       agCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagc                                                                                                                                            <  2:225974/149‑1 (MQ=255)
            atcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTACAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAgg                                                                                                                                         >  8:301422/1‑147 (MQ=255)
               atcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagttcat                                                                                                                                     <  1:282778/148‑6 (MQ=38)
                  gAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCt                                                                                                                                  <  7:178496/148‑1 (MQ=255)
                       cGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAAc                                                                                                                             >  8:130228/1‑148 (MQ=255)
                       cGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAAc                                                                                                                             >  7:21438/1‑148 (MQ=255)
                        gTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGa                                                                                                                           <  4:204464/149‑1 (MQ=255)
                        gTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGa                                                                                                                           <  4:239870/149‑1 (MQ=255)
                           gAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGt                                                                                                                         <  1:344172/148‑1 (MQ=255)
                                ggAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCAGCCCGTGCCGGCCgag                                                                                                                                             >  6:219033/1‑123 (MQ=255)
                                ggAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCAGCCCGTGCCGGCCgag                                                                                                                                             <  5:219033/123‑1 (MQ=255)
                                        aTGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACgg                                                                                                           >  7:36863/1‑149 (MQ=255)
                                                 tACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAg                                                                                                  <  1:271312/149‑1 (MQ=255)
                                                    ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTcc                                                                                               >  6:144204/1‑149 (MQ=255)
                                                    ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTcc                                                                                               >  6:10214/1‑149 (MQ=255)
                                                     cGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCct                                                                                              <  4:92703/149‑1 (MQ=255)
                                                       ccTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTACAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCctgc                                                                                            <  7:301422/149‑1 (MQ=255)
                                                              cTTGCGCGATGAAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCCCATCCCAGACGGGTCAGGAAGTCGTGCTGCTTGt                                                                                     >  8:343451/1‑149 (MQ=255)
                                                                   gcgATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGtt                                                                                <  8:36863/149‑1 (MQ=255)
                                                                          aGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAccc                                                                         >  2:175676/1‑149 (MQ=255)
                                                                            tCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAcccc                                                                        >  6:33491/1‑148 (MQ=255)
                                                                              tAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCACCCCCCTc                                                                     >  4:349833/1‑149 (MQ=255)
                                                                               aGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCg                                                                    >  4:281383/1‑149 (MQ=255)
                                                                                 cccGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAgg                                                                                                     <  6:227467/114‑1 (MQ=255)
                                                                                 cccGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAgg                                                                                                     >  5:227467/1‑114 (MQ=255)
                                                                                   cGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCaa                                                                >  5:145332/1‑149 (MQ=255)
                                                                                        cagTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATa                                                           >  1:277930/1‑149 (MQ=255)
                                                                                           tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACg                                                         >  4:103344/1‑148 (MQ=255)
                                                                                           tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGc                                                        <  3:103344/149‑1 (MQ=255)
                                                                                               ttCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGc                                                    <  5:33491/149‑1 (MQ=255)
                                                                                                 ccGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGcc                                                   >  5:33870/1‑148 (MQ=255)
                                                                                                    ggCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGGTTCGACCCCCTCGGCAACGATCCGCCGGCCCAAc                                               >  8:338603/1‑149 (MQ=255)
                                                                                                       ccATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCaa                                                                <  5:283778/129‑1 (MQ=255)
                                                                                                       ccATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCaa                                                                >  6:283778/1‑129 (MQ=255)
                                                                                                        cATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCt                                           >  5:175039/1‑149 (MQ=255)
                                                                                                            gCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGcc                                       <  8:21438/149‑1 (MQ=255)
                                                                                                                    cTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcga                               <  7:130228/149‑1 (MQ=255)
                                                                                                                     ttGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAAAGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgac                              >  4:75390/1‑149 (MQ=255)
                                                                                                                       gCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgacga                            >  3:96617/1‑149 (MQ=255)
                                                                                                                         tCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGAtt                          >  1:76396/1‑149 (MQ=255)
                                                                                                                           aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAc                                                                           >  4:119258/1‑98 (MQ=255)
                                                                                                                           aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAc                                                                           <  3:119258/98‑1 (MQ=255)
                                                                                                                           aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGc                        <  6:175039/149‑1 (MQ=255)
                                                                                                                           aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTCCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGAAGATTGc                        >  5:163001/1‑149 (MQ=255)
                                                                                                                                ctgctgGGCCAGCCCGTGCCGGCCGAGCAGGTAGCCCTGCAACGAGTCACATCCCAGCCGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa                   >  5:313091/1‑149 (MQ=255)
                                                                                                                                cTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa                   <  3:75390/149‑1 (MQ=255)
                                                                                                                                   cTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaaaaa                >  8:341367/1‑149 (MQ=255)
                                                                                                                                       gccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGc            <  5:144204/149‑1 (MQ=255)
                                                                                                                                        ccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCt           >  4:76406/1‑149 (MQ=255)
                                                                                                                                        ccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCt           <  2:277930/149‑1 (MQ=255)
                                                                                                                                            cACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCAt       >  2:71281/1‑149 (MQ=255)
                                                                                                                                            cACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCACTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTgact       >  8:282571/1‑146 (MQ=255)
                                                                                                                                               gggCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGt    <  6:33870/149‑1 (MQ=255)
                                                                                                                                                 gCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGTCa  >  3:110515/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGTCA  >  NC_002947/4062174‑4062467

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