Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,809,188 T→A 7.6% L80Q (CTG→CAG)  eno → enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,809,1880TA 7.6% 199.3 / 5.3 66L80Q (CTG→CAG) enoenolase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base A (2/3);  total (32/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 = 7.40e-01

ATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACC  >  NC_002947/1809043‑1809325
                                                                                                                                                 |                                                                                                                                         
atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAggt                                                                                                                                        >  6:219311/1‑148 (MQ=255)
   ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGc                                                                                                                                       >  5:290985/1‑147 (MQ=255)
    gCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCacg                                                                                                                                    <  5:219311/149‑3 (MQ=255)
        cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATc                                                                                                                                <  1:167548/149‑1 (MQ=255)
        cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATc                                                                                                                                >  3:81773/1‑149 (MQ=255)
        cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATc                                                                                                                                <  1:271889/149‑1 (MQ=255)
          cTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                              <  4:81773/149‑3 (MQ=255)
           ttGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                             >  7:187082/1‑149 (MQ=255)
             gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                             <  3:34123/147‑1 (MQ=255)
                       cGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAgg                                                                                                                 >  7:172819/1‑149 (MQ=255)
                         gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgctg                                                                                                                                          >  4:119683/1‑122 (MQ=255)
                          tGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGctc                                                                                                              <  8:187082/149‑1 (MQ=255)
                            cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                            <  6:290985/149‑1 (MQ=255)
                            cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                            >  4:174037/1‑149 (MQ=255)
                            cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGAAGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCCTGACCTGCTGCTGGGCAAGGATCCTTCCGAACAGAAGGCTCTg                                                                                                            >  4:248849/1‑149 (MQ=255)
                               cTACCGTCTCGCGCGAAGCGCTGGAGCTGCGTTATGGCGAAAAGAGCCTTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCTATCCGTGACCTGCTGCTTGGCAAGGATCCTTCCGATCAGAATGctcttg                                                                                                           <  4:154196/147‑3 (MQ=255)
                                 aCCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcg                                                                                                       >  2:96100/1‑149 (MQ=255)
                                 aCCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcg                                                                                                       >  3:79111/1‑149 (MQ=255)
                                  ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGAc                                                                                                         <  3:174037/146‑1 (MQ=255)
                                  ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcgc                                                                                                      <  4:79111/149‑1 (MQ=255)
                                         gcgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGa                                                                                               >  5:135565/1‑149 (MQ=255)
                                          cgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGTCCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGAt                                                                                              <  1:202149/149‑1 (MQ=255)
                                           gcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCa                                                                                                  <  6:355967/144‑1 (MQ=255)
                                           gcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCa                                                                                                  >  5:355967/1‑144 (MQ=255)
                                            cgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATc                                                                                             >  3:310027/1‑148 (MQ=255)
                                               aGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         >  1:71712/1‑149 (MQ=255)
                                                 cgcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTg                                                                                       <  5:148704/149‑1 (MQ=255)
                                                      gAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCAGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACgg                                                                                  >  6:66752/1‑149 (MQ=255)
                                                             gTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCgag                                                                           <  4:310027/149‑1 (MQ=255)
                                                                     gACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                                   <  7:150927/149‑1 (MQ=255)
                                                                       cAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGcc                                                                  >  1:9761/1‑148 (MQ=255)
                                                                        aaGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGcc                                                                  >  5:123002/1‑147 (MQ=255)
                                                                          gagCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCg                                                                             >  4:49355/1‑134 (MQ=255)
                                                                          gagCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCg                                                                             <  3:49355/134‑1 (MQ=255)
                                                                           agCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAg                                                                     <  1:145779/141‑1 (MQ=255)
                                                                           agCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAg                                                                     >  2:145779/1‑141 (MQ=255)
                                                                             ccGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         >  6:80331/1‑119 (MQ=255)
                                                                             ccGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         <  5:80331/119‑1 (MQ=255)
                                                                                       ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         <  5:349688/141‑1 (MQ=255)
                                                                                       ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         >  6:349688/1‑141 (MQ=255)
                                                                                       ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACg                                                  <  2:9761/148‑1 (MQ=255)
                                                                                            gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAg                                                               <  4:220213/130‑1 (MQ=255)
                                                                                            gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAg                                                               >  3:220213/1‑130 (MQ=255)
                                                                                             ggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGTCGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCt                                           <  6:123002/149‑1 (MQ=255)
                                                                                                 tGCTGAAGGCCGTCGCCAACATTAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCt                                       >  7:283096/1‑149 (MQ=255)
                                                                                                    tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                                   >  4:142235/1‑118 (MQ=255)
                                                                                                    tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  1:23058/1‑149 (MQ=255)
                                                                                                    tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  7:84797/1‑149 (MQ=255)
                                                                                                    tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  1:254770/1‑149 (MQ=255)
                                                                                                      aaGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    <  5:66752/147‑1 (MQ=255)
                                                                                                          ccGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCt                               <  2:193804/148‑1 (MQ=255)
                                                                                                          ccGTCGCCAACATCAACGGCCCGATCCGAGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCt                               >  1:193804/1‑148 (MQ=255)
                                                                                                           cGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTgg                             <  3:240827/149‑1 (MQ=255)
                                                                                                            gTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGc                            <  8:84797/149‑1 (MQ=255)
                                                                                                            gTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGc                            <  2:254770/149‑1 (MQ=255)
                                                                                                                  aaCATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTgccgc                      >  5:229198/1‑149 (MQ=255)
                                                                                                                    cATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTg                                          >  6:144776/1‑127 (MQ=255)
                                                                                                                    cATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTg                                          <  5:144776/127‑1 (MQ=255)
                                                                                                                     aTCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCa                    >  1:326293/1‑148 (MQ=255)
                                                                                                                        aaCGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAg                 <  8:283096/148‑1 (MQ=255)
                                                                                                                                 aTCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCcca       <  6:229198/149‑1 (MQ=255)
                                                                                                                                   ccGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc        >  1:83451/1‑146 (MQ=255)
                                                                                                                                   ccGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc        <  2:83451/146‑1 (MQ=255)
                                                                                                                                    cGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccagg     <  3:248849/148‑1 (MQ=255)
                                                                                                                                      tGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggacc  <  6:135565/149‑1 (MQ=255)
                                                                                                                                      tGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggac   >  4:15695/1‑148 (MQ=255)
                                                                                                                                       gACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         <  2:231942/93‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
ATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACC  >  NC_002947/1809043‑1809325

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: