Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,142,565 C→G 6.8% A233G (GCC→GGC)  PP_3646 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,142,5650CG 6.8% 235.7 / 6.1 73A233G (GCC→GGC) PP_3646aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base G (2/3);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

ACCATCCTCGAGCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTC  >  NC_002947/4142420‑4142712
                                                                                                                                                  |                                                                                                                                                      
aCCATCCTCGAGCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCGGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCAc                                                                                                                                                     >  1:335405/1‑149 (MQ=255)
       tCGAGCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCAccgcccgc                                                                                                                                              >  8:57416/1‑149 (MQ=255)
          aGCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGGCGGCAccgccggtgaa                                                                                                                                           <  5:75297/149‑7 (MQ=255)
           gCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGGCGGCAccgccggtgaag                                                                                                                                          <  6:222566/149‑8 (MQ=255)
               gCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTgg                                                                                                                                       >  4:159485/1‑148 (MQ=255)
               gCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGt                                                                                                                                      <  8:336873/149‑1 (MQ=255)
                cccGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTa                                                                                                                                     >  8:304615/1‑149 (MQ=255)
                    aCTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAATGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGt                                                                                                                                      <  2:351414/144‑1 (MQ=255)
                    aCTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAATGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGt                                                                                                                                      >  1:351414/1‑144 (MQ=255)
                      tGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACgcag                                                                                                                                >  6:306172/1‑148 (MQ=255)
                           cTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑gcgc                                                                                                                       >  3:110536/1‑149 (MQ=255)
                                     gCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGGACTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCagca‑‑‑gc                                                                                                                         >  1:235990/1‑138 (MQ=255)
                                     gCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCagca‑‑‑gc                                                                                                                         <  2:235990/137‑1 (MQ=255)
                                     gCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTc                                                                                                             <  4:105151/149‑1 (MQ=255)
                                     gCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTc                                                                                                             <  3:4744/149‑1 (MQ=255)
                                          ccggccggTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCaa                                                                                                        <  4:350592/149‑1 (MQ=255)
                                           cggccggTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAg                                                                                                       >  1:157089/1‑149 (MQ=255)
                                                         aaCGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGGCGGCAccgccggt                                                                                                                                              >  1:211228/1‑96 (MQ=255)
                                                         aaCGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGGCGGCAccgccggt                                                                                                                                              <  2:211228/99‑4 (MQ=255)
                                                         aaCGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACt                                                                                         >  6:248995/1‑149 (MQ=255)
                                                         aaCGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACt                                                                                         >  5:203624/1‑149 (MQ=255)
                                                             tGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGggc                                                                                     <  1:59155/149‑1 (MQ=255)
                                                                 cACCGA‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGggcggc                                                                                  >  8:138515/1‑148 (MQ=255)
                                                                  aCCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGggc                                                                                     >  8:218471/1‑144 (MQ=255)
                                                                  aCCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGggc                                                                                     <  7:218471/144‑1 (MQ=255)
                                                                   ccGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAg                                                                               <  7:27586/149‑1 (MQ=255)
                                                                    cGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGt                                                                              >  1:245553/1‑149 (MQ=255)
                                                                    cGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGt                                                                              >  7:281969/1‑149 (MQ=255)
                                                                            ggCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGaa                                                                       >  5:368995/1‑149 (MQ=255)
                                                                            ggCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGaa                                                                       >  8:322381/1‑149 (MQ=255)
                                                                              ccccAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAACTGGGGCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCg                                                                         >  1:303574/1‑148 (MQ=255)
                                                                              ccccAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAAca                                                                     <  7:57416/149‑1 (MQ=255)
                                                                               cccAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAAcat                                                                    >  3:316352/1‑149 (MQ=255)
                                                                                ccAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAAcat                                                                    >  3:9764/1‑148 (MQ=255)
                                                                                 cAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTg                                                                                              <  7:204121/121‑1 (MQ=255)
                                                                                 cAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTg                                                                                              >  8:204121/1‑121 (MQ=255)
                                                                                  aGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAAcatcat                                                                 <  5:306172/149‑1 (MQ=255)
                                                                                   gCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACatcatc                                                                <  7:138515/149‑1 (MQ=255)
                                                                                         gcgcggcgcTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTcgcc                                                          <  2:335405/149‑1 (MQ=255)
                                                                                                gcTGATCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTcgccgacgccg                                                   >  4:196272/1‑149 (MQ=255)
                                                                                                 cTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTcgccgacgccga                                                  <  4:316352/149‑1 (MQ=255)
                                                                                                    aCCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCgccgacgccgac                                                 <  4:110536/147‑1 (MQ=255)
                                                                                                       cGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGa                                            <  2:157089/149‑1 (MQ=255)
                                                                                                       cGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGa                                            >  6:251635/1‑149 (MQ=255)
                                                                                                         ccACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGa                                            >  1:298452/1‑147 (MQ=255)
                                                                                                         ccACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACa                                          <  7:304615/149‑1 (MQ=255)
                                                                                                          cACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAg                                         >  6:56134/1‑149 (MQ=255)
                                                                                                           aCCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAgc                                        >  2:208476/1‑149 (MQ=255)
                                                                                                               gCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCa                                    >  5:263737/1‑149 (MQ=255)
                                                                                                                cTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCAt                                   <  7:322381/149‑1 (MQ=255)
                                                                                                                   gTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATc                                  <  1:208476/147‑1 (MQ=255)
                                                                                                                      cgcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATg                              <  6:263737/148‑1 (MQ=255)
                                                                                                                       gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgc                            <  6:203624/149‑1 (MQ=255)
                                                                                                                       gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgc                            >  8:136557/1‑149 (MQ=255)
                                                                                                                       gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgc                            >  6:367631/1‑149 (MQ=255)
                                                                                                                         aaGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                               <  5:339274/144‑1 (MQ=255)
                                                                                                                         aaGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                               >  6:339274/1‑144 (MQ=255)
                                                                                                                            aTCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGt                       >  5:104484/1‑149 (MQ=255)
                                                                                                                               gCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGc                    >  3:87912/1‑149 (MQ=255)
                                                                                                                               gCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGc                    <  5:248995/149‑1 (MQ=255)
                                                                                                                                ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGcc                   >  6:358654/1‑149 (MQ=255)
                                                                                                                                ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGcc                   <  2:298452/149‑1 (MQ=255)
                                                                                                                                ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGcc                   <  4:9764/149‑1 (MQ=255)
                                                                                                                                ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAACCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGcc                   >  6:363859/1‑149 (MQ=255)
                                                                                                                                 cTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCg                  >  1:328252/1‑149 (MQ=255)
                                                                                                                                    ggcgggcggTGGCGGCACCGCCGGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCa               >  5:24857/5‑149 (MQ=255)
                                                                                                                                      ccGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc             <  5:251635/149‑1 (MQ=255)
                                                                                                                                      ccGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc             <  2:245553/149‑1 (MQ=255)
                                                                                                                                       cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc             >  4:9701/1‑148 (MQ=255)
                                                                                                                                       cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc             >  3:98870/1‑148 (MQ=255)
                                                                                                                                       cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCt            <  2:303574/149‑1 (MQ=255)
                                                                                                                                         gcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTAt          <  3:159485/149‑1 (MQ=255)
                                                                                                                                          cggTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATg         <  5:367631/149‑1 (MQ=255)
                                                                                                                                                cgccACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGt   >  1:335656/1‑149 (MQ=255)
                                                                                                                                                cgccACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGt   >  8:8488/1‑149 (MQ=255)
                                                                                                                                                 gccACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATg                              <  4:119894/121‑1 (MQ=255)
                                                                                                                                                 gccACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATg                              >  3:119894/1‑121 (MQ=255)
                                                                                                                                                 gccACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTc  >  3:329843/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
ACCATCCTCGAGCTGGCCCGACTGGCCCTTGAGGCCGGCTTCCCGGCCGGTGTGGTCAACGTGGTCACCGG‑CTACGGCCCCAGCACCGGCGCGGCGCTGACCCGCCACCCGCTGGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCA‑‑‑GCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTC  >  NC_002947/4142420‑4142712

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