Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,120,236 T→C 6.1% F218L (TTC→CTC)  adhB → alcohol dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,120,2360TC 6.1% 190.9 / 2.3 65F218L (TTC→CTC) adhBalcohol dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTCCTGAAGACCG  >  NC_002947/4120088‑4120383
                                                                                                                                                    |                                                                                                                                                   
cTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCTGCAGTTGCCATACGCCACGTggcctggcct                                                                                                                                                     >  6:333345/1‑149 (MQ=255)
 tGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCACCGGACACTGCGGCAGATGACATACGCCACGTggcctggcct                                                                                                                                                     >  8:264024/1‑148 (MQ=255)
    cccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCaa                                                                                                                                                 >  6:211472/1‑149 (MQ=255)
    cccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCaa                                                                                                                                                 <  6:78229/149‑1 (MQ=255)
     ccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAAt                                                                                                                                                >  8:41337/1‑149 (MQ=255)
     ccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAAt                                                                                                                                                <  3:77949/149‑1 (MQ=255)
     ccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAAt                                                                                                                                                <  1:320154/149‑1 (MQ=255)
     ccATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACACCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAAt                                                                                                                                                >  8:352250/1‑149 (MQ=255)
     ccATCGCCTTGTGGAACGGCCTGGTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAAt                                                                                                                                                <  5:333345/149‑1 (MQ=255)
           ccTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCTGGACACTGCGGCAGTTGCCATACGCAACGTGGCCTGGCCTTCAATGAAAAg                                                                                                                                          <  1:308886/149‑1 (MQ=255)
             ttGTGGAACGGCCTGTTCGCCGCCACTACGCCTTGCACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGgc                                                                                                                                        >  3:221183/1‑149 (MQ=255)
             ttGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGaaaa                                                                                                                                           <  4:197603/146‑1 (MQ=255)
             ttGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGaaaa                                                                                                                                           >  3:197603/1‑146 (MQ=255)
               gtgGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAATGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGgcgc                                                                                                                                      >  2:202774/1‑149 (MQ=255)
                    aCGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGAt                                                                                                                                 <  7:37291/149‑1 (MQ=255)
                       gCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGa                                                                                                                               >  1:88958/1‑148 (MQ=255)
                            ttCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGc                                                                                                                           <  8:27181/147‑1 (MQ=255)
                            ttCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCgg                                                                                                                         >  1:227827/1‑149 (MQ=255)
                                 cgccACTACGCCTTCCACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGc                                                                                                                    <  4:221183/149‑1 (MQ=255)
                                  gccACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGcc                                                                                                                   <  1:115096/149‑1 (MQ=255)
                                         cGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTg                                                                                                            <  7:264024/149‑1 (MQ=255)
                                            cTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTcc                                                                                                         >  7:4596/1‑149 (MQ=255)
                                            cTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTcc                                                                                                         >  8:282578/1‑149 (MQ=255)
                                               acacCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCg                                                                                                        <  5:211472/147‑1 (MQ=255)
                                                cacCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGgcg                                                                                                     <  6:277006/149‑1 (MQ=255)
                                                 acCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAATGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGgcgc                                                                                                    >  8:223390/1‑149 (MQ=255)
                                                 acCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGgcg                                                                                                     >  5:197065/1‑148 (MQ=255)
                                                 acCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGgcg                                                                                                     >  2:200258/1‑148 (MQ=255)
                                                  ccGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCa                                                                                                   >  4:79670/1‑149 (MQ=255)
                                                       aaGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCCTCAATGAAAAGGCGCTGGATGAc                                                                                                                              <  1:397237/117‑1 (MQ=255)
                                                       aaGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCCTCAATGAAAAGGCGCTGGATGAc                                                                                                                              <  1:403578/117‑1 (MQ=255)
                                                       aaGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCCTCAATGAAAAGGCGCTGGATGAc                                                                                                                              >  2:397237/1‑117 (MQ=255)
                                                       aaGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCCTCAATGAAAAGGCGCTGGATGAc                                                                                                                              >  2:403578/1‑117 (MQ=255)
                                                          gCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGa                                                                                           <  7:352250/149‑1 (MQ=255)
                                                              ggCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCGCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGc                                                                                       <  1:384237/149‑1 (MQ=255)
                                                                cAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTg                                                                                     <  1:331782/149‑1 (MQ=255)
                                                                cAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTg                                                                                     >  2:331782/1‑149 (MQ=255)
                                                                        acaGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCg                                                                                                                          >  7:103560/1‑104 (MQ=255)
                                                                        acaGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCAc                                                                             >  8:104230/1‑149 (MQ=255)
                                                                          aGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCAcc                                                                            >  1:2449/1‑148 (MQ=255)
                                                                           gTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCa                                                                          <  1:202774/149‑1 (MQ=255)
                                                                            tGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAg                                                                         >  7:191321/1‑149 (MQ=255)
                                                                              gAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGcc                                                                       >  7:133523/1‑149 (MQ=255)
                                                                                 ccGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTAc                                                                    <  6:197065/149‑1 (MQ=255)
                                                                                                cGTCCAGGGCCCCGGACACTGCGGCAGTTGGCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGa                                                             >  1:32532/1‑141 (MQ=255)
                                                                                                cGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGa                                                             <  2:32532/141‑1 (MQ=255)
                                                                                                  tCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACAcc                                                   >  8:111696/1‑149 (MQ=255)
                                                                                                        gCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTg                                             <  2:2449/149‑1 (MQ=255)
                                                                                                         ccccGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTgc                                            >  7:273784/1‑148 (MQ=255)
                                                                                                          cccGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTggg                                           >  6:132607/1‑149 (MQ=255)
                                                                                                                 acTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATG‑AAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTg                                     >  4:118783/1‑147 (MQ=255)
                                                                                                                  cTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGa                                   <  8:133523/149‑1 (MQ=255)
                                                                                                                   tGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCTCCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAg                                  >  6:340063/1‑149 (MQ=255)
                                                                                                                     cGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTgg                                    <  5:137106/145‑1 (MQ=255)
                                                                                                                     cGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTgg                                    >  6:137106/1‑145 (MQ=255)
                                                                                                                      ggCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGc                                 >  3:305170/1‑147 (MQ=255)
                                                                                                                        cAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAgg                             <  8:273784/149‑1 (MQ=255)
                                                                                                                        cAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAgg                             <  2:88958/149‑1 (MQ=255)
                                                                                                                                aTACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCa                           >  6:197773/1‑143 (MQ=255)
                                                                                                                                 tACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGAtag                    >  6:124790/1‑147 (MQ=255)
                                                                                                                                  aCGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCTCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATcgc                   >  5:141126/1‑149 (MQ=255)
                                                                                                                                   cGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATcgcg                  <  8:191321/149‑1 (MQ=255)
                                                                                                                                    gCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGa                               >  5:15016/1‑135 (MQ=255)
                                                                                                                                    gCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGa                               <  6:15016/135‑1 (MQ=255)
                                                                                                                                     ccACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCa                <  7:111696/149‑1 (MQ=255)
                                                                                                                                     ccACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCa                <  3:79670/149‑1 (MQ=255)
                                                                                                                                         gTGGCCTGGCCTTCAATG‑AAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTcc           <  3:118783/149‑1 (MQ=255)
                                                                                                                                                 gcctTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTCCTGAAGAc    >  5:322398/1‑149 (MQ=255)
                                                                                                                                                   ctTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTCCTGAAGACCg  >  4:328599/1‑149 (MQ=255)
                                                                                                                                                   ctTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTCCTGAAGACCg  >  3:76714/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATATCGTCCAGGGCCCCGGACACTGCGGCAGTTGCCATACGCCACGTGGCCTGGCCTTCAATGAAAAGGCGCTGGATGACAGCGGCAAGCCGTTCCTGTCCGGCGCACTGCTCGACGGCTGGTATGCACCCAGCCTACGCGCCGATCACAACACCGGGCTGGGGCGCTGGAGCGAGGCAGAGATCGCGCAGTTCCTGAAGACCG  >  NC_002947/4120088‑4120383

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