Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,062,307 C→G 8.2% E678Q (GAG→CAG)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,062,3070CG 8.2% 190.2 / 6.5 61E678Q (GAG→CAG) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base G (3/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CCGCCGCCGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAAC  >  NC_002947/4062166‑4062451
                                                                                                                                             |                                                                                                                                                
ccgccgccgATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTcggccgg                                                                                                                                            >  8:111474/1‑148 (MQ=255)
    cgccgATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACgg                                                                                                                                       <  5:354613/149‑1 (MQ=255)
        gATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTg                                                                                                                                   <  7:111474/149‑1 (MQ=255)
                   catcatCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAgg                                                                                                                         >  6:61091/1‑148 (MQ=255)
                    atcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCcagca                                                                                                                           >  8:93183/1‑145 (MQ=255)
                    atcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCcagca                                                                                                                           <  7:93183/145‑1 (MQ=255)
                     tcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAg                                                                                                                      >  8:261890/1‑149 (MQ=255)
                     tcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAg                                                                                                                      >  8:257023/1‑149 (MQ=255)
                     tcatcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAg                                                                                                                      >  8:125356/1‑149 (MQ=255)
                       atcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCAGCCCGTGCCGGccga                                                                                                                              <  7:318300/139‑3 (MQ=255)
                       atcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCAGCCCGTGCCGGccga                                                                                                                              >  8:318300/1‑137 (MQ=255)
                               cGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAAc                                                                                                             >  7:124589/1‑148 (MQ=255)
                                 tCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAg                                                                                                          <  8:124589/149‑1 (MQ=255)
                                   gAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGTCCCAGCAGGTAGCCCTGCAACGAGTc                                                                                                        <  5:61091/149‑1 (MQ=255)
                                        ggAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagt                                                                                                                         <  2:220579/127‑2 (MQ=39)
                                        ggAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagt                                                                                                                         >  1:220579/1‑126 (MQ=39)
                                                aTGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACgg                                                                                           <  7:261890/149‑1 (MQ=255)
                                                aTGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACgg                                                                                           <  7:257023/149‑1 (MQ=255)
                                                    cTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTgg                                                                                                                                  <  8:298463/106‑1 (MQ=255)
                                                    cTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGt                                                                                         >  2:115073/1‑147 (MQ=255)
                                                     ttCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagttc                                                                                                                       >  8:218958/1‑113 (MQ=38)
                                                     ttCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagttc                                                                                                                       <  7:218958/116‑4 (MQ=38)
                                                     ttCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTc                                                                                        >  2:314923/1‑147 (MQ=255)
                                                     ttCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTc                                                                                        <  1:314923/147‑1 (MQ=255)
                                                        atttccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGaa                                                                                   <  1:237545/145‑1 (MQ=255)
                                                        aTACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGaa                                                                                   <  7:125356/149‑1 (MQ=255)
                                                            ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACgg                                                                                           <  7:191255/137‑1 (MQ=255)
                                                            ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACgg                                                                                           >  8:191255/1‑137 (MQ=255)
                                                            ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAATCACATCCCAGACGGGTCAGGAAGTcc                                                                               >  6:158292/1‑149 (MQ=255)
                                                             cGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCct                                                                              >  5:331769/1‑149 (MQ=255)
                                                              gCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCctg                                                                             >  5:60276/1‑149 (MQ=255)
                                                               ccTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGaa                                                                                   >  8:105053/1‑142 (MQ=255)
                                                               ccTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGaa                                                                                   <  7:105053/142‑1 (MQ=255)
                                                                 tttACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCctgctg                                                                          >  2:279305/1‑149 (MQ=255)
                                                                           gcgATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGtt                                                                <  6:331769/149‑1 (MQ=255)
                                                                               tGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCaa                                                            <  1:115073/149‑1 (MQ=255)
                                                                                  aGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCaa                                                            <  3:101558/146‑1 (MQ=255)
                                                                                  aGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCaa                                                            >  4:101558/1‑146 (MQ=255)
                                                                                  aGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAccc                                                         >  6:228717/1‑149 (MQ=255)
                                                                                        gCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                   <  5:228717/149‑1 (MQ=255)
                                                                                                   tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTcggccg                                                                                                                                             >  1:73111/1‑48 (MQ=255)
                                                                                                   tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTcggccg                                                                                                                                             <  2:73111/48‑1 (MQ=255)
                                                                                                   tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGATGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACg                                         >  2:239330/1‑148 (MQ=255)
                                                                                                      gTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAgg                                     >  2:229918/1‑149 (MQ=255)
                                                                                                         ccGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                   <  6:227803/132‑1 (MQ=255)
                                                                                                         ccGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                   >  5:227803/1‑132 (MQ=255)
                                                                                                              cccATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGc                             >  5:72635/1‑149 (MQ=255)
                                                                                                                 aTGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCCCATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTg                          >  2:30236/1‑149 (MQ=255)
                                                                                                                 aTGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGACCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTg                          >  3:329296/1‑149 (MQ=255)
                                                                                                                      cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTg                          >  1:102522/1‑144 (MQ=255)
                                                                                                                      cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTg                          <  2:102522/144‑1 (MQ=255)
                                                                                                                               gCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgacga            >  8:230396/1‑149 (MQ=255)
                                                                                                                                 tCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTcc                                                                               >  7:113329/1‑80 (MQ=255)
                                                                                                                                 tCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTcc                                                                               <  8:113329/80‑1 (MQ=255)
                                                                                                                                 tCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGAtt          >  1:208821/1‑149 (MQ=255)
                                                                                                                                 tCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGAtt          >  8:11230/1‑149 (MQ=255)
                                                                                                                                 tCATGAACTGCTCGGCCGGCACGGGCTGG‑CCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTg         >  7:82331/1‑149 (MQ=255)
                                                                                                                                  cATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGt                   >  5:219003/1‑139 (MQ=255)
                                                                                                                                  cATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGt                   <  6:219003/139‑1 (MQ=255)
                                                                                                                                  cATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGt                   >  5:213831/1‑139 (MQ=255)
                                                                                                                                  cATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGt                   <  6:213831/139‑1 (MQ=255)
                                                                                                                                        cTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa   <  2:208821/149‑1 (MQ=255)
                                                                                                                                         tGCTGGGCCAGCCCGTGCCGGCCGAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa   >  3:142862/1‑148 (MQ=255)
                                                                                                                                          gCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAAc  >  2:90683/1‑148 (MQ=255)
                                                                                                                                             cggccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgacga            <  2:152043/135‑1 (MQ=255)
                                                                                                                                             cggccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgacga            >  1:152043/1‑135 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CCGCCGCCGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAAC  >  NC_002947/4062166‑4062451

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