Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,469,023 G→C 5.8% A408A (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,0230GC 5.8% 261.5 / 7.7 86A408A (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.25e-01

GATGTTGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGAGCAGC  >  NC_002947/4468874‑4469166
                                                                                                                                                     |                                                                                                                                               
gATGTTGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcg                                                                                                                                                 <  2:1803659/150‑1 (MQ=255)
    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)
                                 aTGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcggc                                                                                                                                               >  2:1349881/1‑119 (MQ=255)
                                 aTGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATggcggc                                                                                                                                               <  1:1349881/119‑1 (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                                                                                                                                      <  2:126347/94‑1 (MQ=255)
                                                                   cACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCAtctc                                                                                                                                      >  1:126347/1‑94 (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                                                                                                           <  2:2091123/103‑1 (MQ=255)
                                                                                     tCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTa                                                                                                           >  1:2091123/1‑103 (MQ=255)
                                                                                      cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  1:1018312/150‑1 (MQ=255)
                                                                                      cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  2:231723/150‑1 (MQ=255)
                                                                                      cGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGtc                                                           <  2:2321552/150‑1 (MQ=255)
                                                                                        tGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGTCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAc                                                                                                          >  1:814017/1‑101 (MQ=255)
                                                                                        tGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGTCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTAc                                                                                                          <  2:814017/101‑1 (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                                         <  2:1418572/148‑1 (MQ=255)
                                                                                                          gCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACccgc                                         >  1:1418572/1‑148 (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                                                                                                <  2:1350534/67‑1 (MQ=255)
                                                                                                                                    gCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTgg                                                                                                >  1:1350534/1‑67 (MQ=255)
                                                                                                                                     ccGCGACACTGATGGCGGCGGCCATCTCTGGATgg                                                                                                                               <  1:824112/35‑1 (MQ=255)
                                                                                                                                     ccGCGACACTGATGGCGGCGGCCATCTCTGGATgg                                                                                                                               >  2:824112/1‑35 (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                                                                    <  2:1335899/85‑1 (MQ=255)
                                                                                                                                              tGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCAtgc                                                                    >  1:1335899/1‑85 (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)
                                                                                                                                                     |                                                                                                                                               
GATGTTGACCCTCAGGCAGGCGCACTACTTGGCATGGGCGCATACCTGGCCGGTGTGACACGCTCGCCACTGACCGCCTCGGTCATCGTGCTTGAACTGACACACAGCCCGGACCTGGCCATTCCCATGTTGGCCGCGACACTGATGGCGGCGGCCATCTCTGGATGGGTGTCTCCGGTCTCGCTTTACCATGCCTTGGCCAAGCAGATTACCGACAAGCTTCATGCTGCCGCGTCTTCTTCCACAGAACCCGCCGCTCGCCAACCGAGCACTACCGAAAAACCTTGAGCAGC  >  NC_002947/4468874‑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: