Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,469,026 G→C 5.8% A409A (GCG→GCC PP_3959 → voltage‑gated chloride channel family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,469,0260GC 5.8% 262.8 / 7.7 86A409A (GCG→GCCPP_3959voltage‑gated chloride channel family protein
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base C (2/3);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

TTGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGAGCAGC  >  NC_002947/4468878‑4469166
                                                                                                                                                    |                                                                                                                                            
ttGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTTGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcggcgg                                                                                                                                             >  1:68722/1‑150 (MQ=255)
  gACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGcc                                                                                                                                           >  2:1283090/1‑150 (MQ=255)
     ccTCAGGCAGGCGCACTACTTGGTATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtc                                                                                                                                        >  1:457287/1‑150 (MQ=255)
     ccTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcggcgg                                                                                                                                             <  1:191744/145‑1 (MQ=255)
     ccTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcggcgg                                                                                                                                             >  2:191744/1‑145 (MQ=255)
     ccTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtc                                                                                                                                        >  1:1871008/1‑150 (MQ=255)
      cTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtct                                                                                                                                       >  1:625826/1‑150 (MQ=255)
        caggcaggCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCCGCCGCCAtcagt                                                                                                                                     <  2:2272718/150‑4 (MQ=255)
                     tgctTGGCGGGGGCGAATACCTGGTCGGGGTGACACGGTCTCGACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTAGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGtctc                                                                                                                        >  2:1895731/3‑150 (MQ=255)
                     tactTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGtctc                                                                                                                        >  2:762580/1‑150 (MQ=255)
                     tactTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGtctc                                                                                                                        >  2:2202697/1‑150 (MQ=255)
                       ctTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCg                                                                                                                      >  2:1492988/1‑150 (MQ=255)
                            cATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGtctc                                                                                                                 >  2:1018312/1‑150 (MQ=255)
                             aTGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCg                                                                                                                >  1:2250335/1‑150 (MQ=255)
                             aTGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCg                                                                                                                >  1:1156321/1‑150 (MQ=255)
                                 gcgcATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGTGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCttt                                                                                                            <  2:1950927/150‑1 (MQ=255)
                                  cgcATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTa                                                                                                           <  2:757900/150‑1 (MQ=255)
                                    cATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAcc                                                                                                         <  2:457287/150‑1 (MQ=255)
                                    cATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAcc                                                                                                         <  1:2134475/150‑1 (MQ=255)
                                     aTACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCCGCCGCCAtc                                                                                                                                        >  1:1639101/1‑118 (MQ=255)
                                     aTACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCCGCCGCCAtc                                                                                                                                        <  2:1639101/118‑1 (MQ=255)
                                            gCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTg                                                                                                 <  2:1610404/150‑1 (MQ=255)
                                                 tgtgACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCaa                                                                                            >  1:2190391/1‑150 (MQ=255)
                                                 tgtgACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCaa                                                                                            <  1:256712/150‑1 (MQ=255)
                                                  gtgACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAg                                                                                           <  2:1871008/150‑1 (MQ=255)
                                                     acacGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATggg                                                                                                                              >  2:1037105/1‑112 (MQ=255)
                                                     acacGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATggg                                                                                                                              <  1:1037105/112‑1 (MQ=255)
                                                       acGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGAt                                                                                      >  2:251502/1‑150 (MQ=255)
                                                         gCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTa                                                                                    >  1:1813014/1‑150 (MQ=255)
                                                          cTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCAt                                                                                                       >  1:248890/1‑130 (MQ=255)
                                                          cTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCAt                                                                                                       <  2:248890/130‑1 (MQ=255)
                                                           tCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTAcc                                                                                  <  1:1283090/150‑1 (MQ=255)
                                                           tCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCCGCCGCCAtcagt                                                                                                                                     >  1:281330/1‑96 (MQ=255)
                                                           tCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCCGCCGCCAtcagt                                                                                                                                     <  2:281330/99‑4 (MQ=255)
                                                           gcGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTAcc                                                                                  >  2:1504999/2‑150 (MQ=255)
                                                              ccACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGAc                                                                               >  1:843583/1‑150 (MQ=255)
                                                               cACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtctc                                                                                                                                      >  1:126347/1‑94 (MQ=255)
                                                               cACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtctc                                                                                                                                      <  2:126347/94‑1 (MQ=255)
                                                               cACTGACCACCTCTGTCATCGGGCTTGAACTTACACACAGCCCGGACCTGGCCATTCCCATGTAGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGCTGTCTCCGGTCTCGCTGTACCGTGCCTTGGCCAAGCAGATTACCGACa                                                                              >  2:1097324/1‑150 (MQ=255)
                                                                aCTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACaa                                                                             >  1:1578639/1‑150 (MQ=255)
                                                                  tGACCGCCTCGGTCATGGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCTACTGCGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGc                                                                           <  1:987431/150‑1 (MQ=255)
                                                                  tGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGc                                                                           <  2:843583/150‑1 (MQ=255)
                                                                      cGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCa                                                                       >  1:231723/1‑150 (MQ=255)
                                                                          tCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCAtgct                                                                   >  1:2321552/1‑150 (MQ=255)
                                                                          tCGGTCATAGAGCTTGAACGGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCAtgct                                                                   <  2:1813014/150‑1 (MQ=255)
                                                                                 tCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTa                                                                                                           >  1:2091123/1‑103 (MQ=255)
                                                                                 tCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTa                                                                                                           <  2:2091123/103‑1 (MQ=255)
                                                                                  cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  2:2321552/150‑1 (MQ=255)
                                                                                  cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  2:231723/150‑1 (MQ=255)
                                                                                  cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  1:1018312/150‑1 (MQ=255)
                                                                                    tGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGTCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAc                                                                                                          <  2:814017/101‑1 (MQ=255)
                                                                                    tGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGTCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAc                                                                                                          >  1:814017/1‑101 (MQ=255)
                                                                                    tGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtctt                                                         <  2:2190391/150‑1 (MQ=255)
                                                                                     gCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCAt                                                                                                       >  1:514262/1‑103 (MQ=255)
                                                                                     gCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCAt                                                                                                       <  2:514262/103‑1 (MQ=255)
                                                                                       ttGAACTGACACACAGCCCGGACGTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCGATCTCTGGATGGGTGTCTCCGGTCTCGCGTTACCAACCCCTGGCCAAGCAGATTACCGACAAACTTCATGCTGCGGCGTgttctt                                                      >  2:264083/1‑150 (MQ=255)
                                                                                       ttGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGtctc                                                                                                                        >  1:2270728/1‑84 (MQ=255)
                                                                                       ttGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGtctc                                                                                                                        <  2:2270728/84‑1 (MQ=255)
                                                                                         gAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTcc                                                    <  2:68722/150‑1 (MQ=255)
                                                                                           aCTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCcac                                                  <  1:1504999/150‑1 (MQ=255)
                                                                                             cgACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAg                                                >  2:1771144/2‑150 (MQ=255)
                                                                                                   acaGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCTTTTTCTTCCACAGAACccg                                          >  1:1284700/1‑150 (MQ=255)
                                                                                                    caGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACccgc                                         <  1:251502/150‑1 (MQ=255)
                                                                                                      gCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACccgc                                         >  1:1418572/1‑148 (MQ=255)
                                                                                                      gCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACccgc                                         <  2:1418572/148‑1 (MQ=255)
                                                                                                              cTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCaa                               >  2:1314388/1‑150 (MQ=255)
                                                                                                                        ccATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTg                                                                                                 >  1:16742/1‑74 (MQ=255)
                                                                                                                        ccATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTg                                                                                                 <  2:16742/74‑1 (MQ=255)
                                                                                                                        ccATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTa                     >  2:1306889/1‑150 (MQ=255)
                                                                                                                           tGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCg                  >  2:432112/1‑150 (MQ=255)
                                                                                                                                gCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTgg                                                                                                >  1:1350534/1‑67 (MQ=255)
                                                                                                                                gCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTgg                                                                                                <  2:1350534/67‑1 (MQ=255)
                                                                                                                                 ccGCGACACTGATGGCGGCGGCCATCTCTGGATgg                                                                                                                               >  2:824112/1‑35 (MQ=255)
                                                                                                                                 ccGCGACACTGATGGCGGCGGCCATCTCTGGATgg                                                                                                                               <  1:824112/35‑1 (MQ=255)
                                                                                                                                 ccGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCa                                                                                                        >  2:92702/1‑58 (MQ=255)
                                                                                                                                 ccGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCa                                                                                                        <  1:92702/58‑1 (MQ=255)
                                                                                                                                    cgACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCtt         <  1:1771144/150‑1 (MQ=255)
                                                                                                                                    cgACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCtt         <  1:1492988/150‑1 (MQ=255)
                                                                                                                                    cgACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCtt         <  2:1578639/150‑1 (MQ=255)
                                                                                                                                     gACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTg        <  1:762580/150‑1 (MQ=255)
                                                                                                                                      acacTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGa       <  2:625826/150‑1 (MQ=255)
                                                                                                                                          tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCAtgc                                                                    >  1:1335899/1‑85 (MQ=255)
                                                                                                                                          tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCAtgc                                                                    <  2:1335899/85‑1 (MQ=255)
                                                                                                                                          tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTg        >  1:1415998/1‑145 (MQ=255)
                                                                                                                                          tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTg        <  2:1415998/145‑1 (MQ=255)
                                                                                                                                          tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGagcag   <  2:1156321/150‑1 (MQ=255)
                                                                                                                                           gATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGTTCTCGCTTGACCATGCCTTGGCAACGCAGAGTACCGACAAGCTTCATGCTGCCGCCTCTTCTTCCACAGACCCCGCCGCTCGCCAGCCTAGCACTACCGAAAAACCTTGaccagc  <  2:938770/150‑1 (MQ=255)
                                                                                                                                             tGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCATATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCtt         >  2:795751/1‑141 (MQ=255)
                                                                                                                                             tGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCtt         <  1:795751/141‑1 (MQ=255)
                                                                                                                                                    ggcCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTcc                                                    <  2:539414/91‑1 (MQ=255)
                                                                                                                                                    ggcCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTcc                                                    >  1:539414/1‑91 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TTGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGAGCAGC  >  NC_002947/4468878‑4469166

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