Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,467,800 A→T 7.8% S15S (TCA→TCT)  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,467,8000AT 7.8% 181.3 / 5.5 64S15S (TCA→TCT) PP_3959voltage‑gated chloride channel family protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base T (2/3);  total (30/34)
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.44e-01

AAGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCTCTTCT  >  NC_002947/4467654‑4467946
                                                                                                                                                   |                                                                                                                                                  
aaGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACg                                                                                                                                                  <  3:135142/149‑1 (MQ=255)
         gCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCgag                                                                                                                                          >  5:338560/1‑148 (MQ=255)
          ccACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGACgccat                                                                                                                                        >  2:15239/1‑145 (MQ=255)
           cACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGTCCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAAt                                                                                                                                       <  4:152494/149‑1 (MQ=255)
            aCTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTCACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAAt                                                                                                                                       <  1:261666/149‑1 (MQ=255)
            aCTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGaacg                                                                                                                                      >  2:208966/1‑147 (MQ=255)
              tCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGc                                                                                                                                    >  8:56121/1‑149 (MQ=255)
                ccGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGc                                                                                                                                    <  1:208966/147‑1 (MQ=255)
                ccGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                  <  6:166871/149‑2 (MQ=255)
                       ggCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGAc                                                                                                                           <  6:207032/149‑1 (MQ=255)
                        gCTACAACCTGTGTCGGGTAGGCCTGTGCTGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACt                                                                                                                          <  7:56121/149‑1 (MQ=255)
                        gCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACt                                                                                                                          <  7:33902/149‑1 (MQ=255)
                                cTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCt                                                                                                                  <  7:182576/149‑1 (MQ=255)
                                   tgtCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTg                                                                                                                 <  6:76826/147‑1 (MQ=255)
                                       gggTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCgg                                                                                                           <  3:82128/149‑1 (MQ=255)
                                        ggTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCggg                                                                                                          <  2:101969/149‑1 (MQ=255)
                                              ccTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCt                                                                                                                  <  4:204979/135‑1 (MQ=255)
                                              ccTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCt                                                                                                                  >  3:204979/1‑135 (MQ=255)
                                                  tgCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCggg                                                                                                          >  3:323027/1‑139 (MQ=255)
                                                  tgCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCggg                                                                                                          <  4:323027/139‑1 (MQ=255)
                                                  tgCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCtgg                                                                                                >  6:274549/1‑149 (MQ=255)
                                                     gggAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGTGTTCTGGATCGGGGCATTGCtggtg                                                                                              <  3:358987/148‑1 (MQ=255)
                                                     gggAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGTGTTCTGGATCGGGGCATTGCtggtg                                                                                              >  4:358987/1‑148 (MQ=255)
                                                       gAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGAGAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTgggg                                                                                           >  1:254141/1‑149 (MQ=255)
                                                       gAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCCTACACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTCGATCGGGGCATTGCtggtggcc                                                                                           >  2:305615/1‑147 (MQ=255)
                                                        aaaaCACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGc                                                                                          <  7:7825/149‑1 (MQ=255)
                                                               gCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTgc                                                                                   >  7:163192/1‑149 (MQ=255)
                                                                      tCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCt                                                                            >  3:24727/1‑149 (MQ=255)
                                                                      tCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCt                                                                            >  3:341128/1‑149 (MQ=255)
                                                                       cACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                  <  6:299177/94‑2 (MQ=37)
                                                                       cACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                  >  5:299177/1‑93 (MQ=37)
                                                                        aCGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCttt                                                                          >  8:67724/1‑149 (MQ=255)
                                                                        aCGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCttt                                                                          <  8:163192/149‑1 (MQ=255)
                                                                          ggCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCg                                                                        >  2:181815/1‑149 (MQ=255)
                                                                                          cATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGAcgc                                                        <  2:174754/149‑1 (MQ=255)
                                                                                                  gACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTc                                                 <  5:147855/148‑1 (MQ=255)
                                                                                                  gACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTc                                                 >  6:147855/1‑148 (MQ=255)
                                                                                                  gACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCg                                                <  6:338560/149‑1 (MQ=255)
                                                                                                   aCAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGc                                               >  1:11456/1‑149 (MQ=255)
                                                                                                     aaTGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGc                                                                             <  4:24929/117‑1 (MQ=255)
                                                                                                     aaTGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGc                                                                             >  3:24929/1‑117 (MQ=255)
                                                                                                     aaTGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCtt                                             >  7:159397/1‑149 (MQ=255)
                                                                                                       tGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCg                                           >  2:95471/1‑149 (MQ=255)
                                                                                                          ccGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGGGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGtt                                         >  7:91551/1‑148 (MQ=255)
                                                                                                            gACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCaa                                      <  4:24727/149‑1 (MQ=255)
                                                                                                              cccT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgc                                    <  2:11456/149‑1 (MQ=255)
                                                                                                              cccT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgc                                    >  8:265554/1‑149 (MQ=255)
                                                                                                                cT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGc                                                                                                                                    >  5:17629/1‑51 (MQ=255)
                                                                                                                cT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGc                                                                                                                                    <  6:17629/51‑1 (MQ=255)
                                                                                                                cT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgcgc                                  <  1:181815/149‑1 (MQ=255)
                                                                                                                    cccGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTg                               <  4:341128/149‑1 (MQ=255)
                                                                                                                       gACCTGCACCATCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTACTGGATCGGGGAATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTg                            >  4:76357/1‑149 (MQ=255)
                                                                                                                            gCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATg                       >  7:32644/1‑149 (MQ=255)
                                                                                                                              aCCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATggg                     <  8:91551/149‑1 (MQ=255)
                                                                                                                                cTTCACTCCGCCGCCGTTCACGCTGGCGAGTATGGCGTCAGCGTCTGGCGTTCTGGATCGGGGCTTTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGc                                                                  <  6:247363/102‑1 (MQ=255)
                                                                                                                                 ttCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCAc                  >  5:49938/1‑149 (MQ=255)
                                                                                                                                    aCTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCt               >  3:164800/1‑149 (MQ=255)
                                                                                                                                       cccccGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTg                            <  2:295963/133‑1 (MQ=255)
                                                                                                                                       cccccGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTg                            >  1:295963/1‑133 (MQ=255)
                                                                                                                                           cgccgTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAAt                        >  7:325289/1‑133 (MQ=255)
                                                                                                                                           cgccgTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAAt                        <  8:325289/133‑1 (MQ=255)
                                                                                                                                             ccgTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACctc      <  5:274549/149‑1 (MQ=255)
                                                                                                                                               gTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCtctt    <  8:32644/149‑1 (MQ=255)
                                                                                                                                                  cACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCtcttct  <  4:153235/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCT‑ACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCTCTTCT  >  NC_002947/4467654‑4467946

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