Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,717,746 A→T 7.4% M18K (ATG→AAG) ‡ paaH ← 3‑hydroxyacyl‑CoA dehydrogenase PaaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,717,7460AT 7.4% 240.1 / 7.5 81M18K (ATG→AAG) ‡paaH3‑hydroxyacyl‑CoA dehydrogenase PaaC
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (4/2);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

AGCGGCCAGTAACTTTCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCG  >  NC_002947/3717601‑3717894
                                                                                                                                                 |                                                                                                                                                      
aGCGGCCAGTAACTTTCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATg                                                                                                                                                      <  3:213374/148‑1 (MQ=255)
      cAGTAACTTTCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCATGGGCGCCgg                                                                                                                                               <  2:323756/149‑1 (MQ=255)
          aaCTTTCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCAcc                                                                                                                                           >  4:43693/1‑149 (MQ=255)
             tttCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGa                                                                                                                                         >  3:334970/1‑148 (MQ=255)
             tttCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGAt                                                                                                                                       >  8:350119/1‑149 (MQ=255)
             tttCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCAACCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGAt                                                                                                                                       >  2:328795/1‑149 (MQ=255)
             tttCCTTTGTCCACCAGCCGGGCCAGTTGGAGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGAt                                                                                                                                       >  6:92116/1‑149 (MQ=255)
                 cTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCAcc                                                                                                                                   <  1:328795/149‑1 (MQ=255)
                  tttGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGgcgc                                                                                                                                                  <  1:226559/134‑1 (MQ=255)
                  tttGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGgcgc                                                                                                                                                  >  2:226559/1‑134 (MQ=255)
                     gTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCAc                                                                                                                                            <  1:339527/137‑1 (MQ=255)
                     gTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCAc                                                                                                                                            >  2:339527/1‑137 (MQ=255)
                      tCCACGAGCCGTGCCAGTTGGCGATCGATTCCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTTGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCAc                                                                                                                              <  6:47271/149‑1 (MQ=255)
                      tCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCAc                                                                                                                              <  7:350119/149‑1 (MQ=255)
                       ccaccaGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCAcc                                                                                                                             <  3:10833/149‑1 (MQ=255)
                        caccaGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCt                                                                                                                            >  7:303272/1‑149 (MQ=255)
                         acccGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCGGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGACCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTg                                                                                                                           <  6:286156/149‑1 (MQ=255)
                         accaGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTg                                                                                                                           >  3:184891/1‑149 (MQ=255)
                          ccaGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTg                                                                                                                           >  1:84204/1‑148 (MQ=255)
                              ccGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACAt                                                                                                                      >  6:311453/1‑149 (MQ=255)
                                   ccAGTTGGGGATGGTTACCCGTGAGCGCCTGGGCGGCAGCACCCGGGGGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACAttgctt                                                                                                                 <  6:92978/149‑1 (MQ=255)
                                     aGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACAttgcttgc                                                                                                               >  1:102555/1‑149 (MQ=255)
                                       ttGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGa                                                                                                             <  8:326109/149‑1 (MQ=255)
                                       ttGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGa                                                                                                             <  8:262343/149‑1 (MQ=255)
                                         ggCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGt                                                                                                           >  8:155351/1‑149 (MQ=255)
                                            gatcgatACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACAGTGCTTGCGAGTGCg                                                                                                        >  4:81990/1‑149 (MQ=255)
                                             atcgatACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCgg                                                                                                       >  8:19326/1‑149 (MQ=255)
                                                 atACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCa                                                                                                    <  7:279578/148‑1 (MQ=255)
                                                 atACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCa                                                                                                    >  8:279578/1‑148 (MQ=255)
                                                     ccGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCAtttt                                                                                                >  4:172196/1‑148 (MQ=255)
                                                     ccGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCAtttt                                                                                                >  3:48357/1‑148 (MQ=255)
                                                      cGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCCCCGGCACATTGCTTGCGAGTGCGGCCCtttt                                                                                                >  6:261023/1‑147 (MQ=255)
                                                          aCCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAgcgc                                                                                         >  6:154631/1‑149 (MQ=255)
                                                          aCCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAgcg                                                                                          <  3:43693/148‑1 (MQ=255)
                                                             gCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGcccc                                                                                      <  2:84204/149‑1 (MQ=255)
                                                                  ggcggcAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGaa                                                                                 <  1:301395/149‑1 (MQ=255)
                                                                           aCCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGgcgt                                                                        >  5:243734/1‑149 (MQ=255)
                                                                           aCCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCATGGGCGCCggta                                                                                                                                             >  6:317441/1‑80 (MQ=255)
                                                                           aCCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCATGGGCGCCggta                                                                                                                                             <  5:317441/82‑3 (MQ=255)
                                                                            cccGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGgcgtg                                                                       >  4:123825/1‑149 (MQ=255)
                                                                             ccGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGgcgtgc                                                                      >  2:126634/1‑149 (MQ=255)
                                                                              cGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTg                                                                                   <  4:196920/135‑1 (MQ=255)
                                                                              cGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTg                                                                                   >  3:196920/1‑135 (MQ=255)
                                                                                 gCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGttt                                                                  <  2:219461/149‑1 (MQ=255)
                                                                                    gTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGtt                                                               <  4:184891/149‑1 (MQ=255)
                                                                                     ttATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGtt                                                               >  8:294886/1‑148 (MQ=255)
                                                                                     ttATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGtt                                                               >  2:42036/1‑148 (MQ=255)
                                                                                     ttATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTc                                                              <  5:311453/149‑1 (MQ=255)
                                                                                     ttATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTc                                                              <  7:19326/149‑1 (MQ=255)
                                                                                       aTCGTAGAGCTTGACGTGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCa                                                             <  8:66306/148‑1 (MQ=255)
                                                                                        tCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATg                                                           <  5:290650/149‑1 (MQ=255)
                                                                                           tAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAg                                                        <  3:123825/149‑1 (MQ=255)
                                                                                            agagCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGttt                                                                  <  7:350180/138‑1 (MQ=255)
                                                                                            agagCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGttt                                                                  >  8:350180/1‑138 (MQ=255)
                                                                                             gagCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGc                                                      <  1:126634/149‑1 (MQ=255)
                                                                                                  tGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGa                                                                                  <  8:143081/116‑1 (MQ=255)
                                                                                                  tGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGa                                                                                  >  7:143081/1‑116 (MQ=255)
                                                                                                    aCGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGAGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACg                                               >  5:307956/1‑149 (MQ=255)
                                                                                                      ggggTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGcc                                             >  1:135770/1‑149 (MQ=255)
                                                                                                         gTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTc                                          >  1:353238/1‑149 (MQ=255)
                                                                                                          tGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCAGATTGCGTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTca                                         <  5:261023/149‑2 (MQ=255)
                                                                                                          tGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGCGAACGCCTTca                                         <  4:334970/149‑2 (MQ=255)
                                                                                                             ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGt                                                                <  6:223014/123‑1 (MQ=255)
                                                                                                             ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGt                                                                >  5:223014/1‑123 (MQ=255)
                                                                                                             ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                      <  5:154631/149‑1 (MQ=255)
                                                                                                                      gCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCg                             >  6:113492/1‑149 (MQ=255)
                                                                                                                        cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                      <  4:164748/138‑1 (MQ=255)
                                                                                                                        cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                      >  3:164748/1‑138 (MQ=255)
                                                                                                                         gcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTg                          >  6:176825/1‑149 (MQ=255)
                                                                                                                            aCCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcgc                       <  2:353238/149‑1 (MQ=255)
                                                                                                                             ccTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGCGCGGCCATTTTCAGCGcccc                                                                                      <  7:310927/85‑1 (MQ=255)
                                                                                                                             ccTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGCGCGGCCATTTTCAGCGcccc                                                                                      >  8:310927/1‑85 (MQ=255)
                                                                                                                                 ggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGc                  >  5:244908/1‑149 (MQ=255)
                                                                                                                                      tgccGGCGCCATGGGCGCCGGCACCGAGTCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAgg              >  1:204115/3‑149 (MQ=255)
                                                                                                                                        ccGGCGCCATGGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGgcgc           >  6:376216/1‑149 (MQ=255)
                                                                                                                                        ccGGCGCCATGGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGgcgc           >  3:67097/1‑149 (MQ=255)
                                                                                                                                         cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcgc                       >  4:158323/1‑136 (MQ=255)
                                                                                                                                         cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcgc                       <  3:158323/136‑1 (MQ=255)
                                                                                                                                           gcgcCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACACCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTg         >  3:196513/1‑148 (MQ=255)
                                                                                                                                              cccATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAgg     >  3:58446/1‑149 (MQ=255)
                                                                                                                                                cATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCTAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGt    >  1:101764/1‑149 (MQ=255)
                                                                                                                                                 aTGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTcg  >  3:12518/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
AGCGGCCAGTAACTTTCCTTTGTCCACCAGCCGGGCCAGTTGGCGATCGATACCCGTCACCGCCTGGGCGGCAGCACCCGGGCGGTTATCGTAGAGCTTGACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCG  >  NC_002947/3717601‑3717894

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: