Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,822,941 G→C 7.5% A157A (GCG→GCC curA → NADPH‑dependent curcumin/dihydrocurcumin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,822,9410GC 7.5% 290.2 / 12.1 93A157A (GCG→GCCcurANADPH‑dependent curcumin/dihydrocurcumin reductase
Reads supporting (aligned to +/- strand):  ref base G (50/36);  new base C (3/4);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGACTT  >  NC_002947/2822802‑2823090
                                                                                                                                           |                                                                                                                                                     
gCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGc                                                                                                                                             >  2:1686015/1‑150 (MQ=255)
 cGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccaccag                                                                                                                                            <  2:1597050/150‑6 (MQ=255)
 cGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccaccag                                                                                                                                            <  2:1829387/150‑6 (MQ=255)
  gATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTggcggcggc                                                                                                                                                     >  1:204861/1‑140 (MQ=255)
  gATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTggcggcggc                                                                                                                                                     <  2:204861/140‑1 (MQ=255)
  gATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCg                                                                                                                                           >  2:809566/1‑150 (MQ=255)
  gATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccaccagg                                                                                                                                           <  2:1590700/150‑7 (MQ=255)
  gATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccaccagg                                                                                                                                           >  1:1590700/1‑144 (MQ=255)
         ccGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcg                                                                                                                                    >  1:792414/1‑150 (MQ=255)
         ccGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcg                                                                                                                                    <  2:935430/150‑1 (MQ=255)
          cgcgTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcgg                                                                                                                                   >  1:57716/1‑149 (MQ=255)
             gTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGc                                                                                                                                             <  2:1969101/137‑1 (MQ=255)
             gTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGc                                                                                                                                             >  1:1969101/1‑137 (MQ=255)
             gTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGAc                                                                                                                                >  1:1964465/1‑150 (MQ=255)
              ttGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACt                                                                                                                               >  2:1423905/1‑150 (MQ=255)
                gggTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACtgt                                                                                                                             >  1:2360380/1‑150 (MQ=255)
                 ggTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACtgtg                                                                                                                            >  1:79840/1‑150 (MQ=255)
                  gTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgg                                                                                                                           <  1:1699116/150‑1 (MQ=255)
                  gTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgg                                                                                                                           >  1:861523/1‑150 (MQ=255)
                  gTGACTCGACGGCGCACACTTCTTATGCCTTTGGCATCTTTGTCATGCAGGGCTTCGGAGCTTATAAGGGCTTGCTCGATCTTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgg                                                                                                                           <  2:684106/150‑1 (MQ=255)
                   tGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTggg                                                                                                                          >  2:557431/1‑150 (MQ=255)
                   tGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTggg                                                                                                                          <  2:513467/150‑1 (MQ=255)
                    gAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgggg                                                                                                                         <  2:2360380/150‑1 (MQ=255)
                     aGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGATTGCTCGATATTGGCCGGCCGCAGGCGGGTTAAACACTGGTGGTGGCGGCGGCCCCCGGGCCGGGTGGCGCGGCTGTGGGGc                                                                                                                        >  2:588278/1‑150 (MQ=255)
                      gTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGAATGTGGGGCa                                                                                                                       >  1:34686/1‑150 (MQ=255)
                      gTCGCCGGCGCATCCTTCTTATGCCCTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGc                                                                                                                                             >  1:1868337/1‑128 (MQ=255)
                      gTCGCCGGCGCATCCTTCTTATGCCCTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGc                                                                                                                                             <  2:1868337/128‑1 (MQ=255)
                         gCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAAACCTGGTGGTGGCGGCGGACACGGGGCCGGTCGGCGAGACTGTGGGGCAGAt                                                                                                                    >  1:877406/1‑150 (MQ=255)
                           cGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCg                                                                                                                  <  2:861523/150‑1 (MQ=255)
                             gcgcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGc                                                                                                                >  2:698393/1‑150 (MQ=255)
                               gcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgggg                                                                                                                         >  1:957443/1‑139 (MQ=255)
                               gcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTgggg                                                                                                                         <  2:957443/139‑1 (MQ=255)
                               gcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCaa                                                                                                              <  1:1423905/150‑1 (MQ=255)
                                    cttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGCCGCGACTGTGGGGCAGATCGGCAAGATCa                                                                                                         >  1:1108501/1‑150 (MQ=255)
                                     ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                        <  2:1657006/150‑1 (MQ=255)
                                     ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                        <  2:792414/150‑1 (MQ=255)
                                     ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                        <  1:809566/150‑1 (MQ=255)
                                     ttcttATGCCTTGGGCATCCTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                        >  2:1786587/1‑150 (MQ=255)
                                        ttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCg                                                                                                                                                >  2:1804/1‑107 (MQ=255)
                                        ttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCg                                                                                                                                                <  1:1804/107‑1 (MQ=255)
                                        ttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAggg                                                                                                     >  2:551251/1‑150 (MQ=255)
                                             ccTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGcc                                                                                                <  1:557431/150‑1 (MQ=255)
                                                tGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAtg                                                                                             <  2:1619517/150‑1 (MQ=255)
                                                   gCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTgg                                                                                          <  2:34686/150‑1 (MQ=255)
                                                   gCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTgg                                                                                          >  1:2033633/1‑150 (MQ=255)
                                                   gCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCAGCAAGATCAAGGGCTGCCATGTgg                                                                                          >  2:377224/1‑150 (MQ=255)
                                                     aTCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCg                                                                                                                                           >  1:527610/1‑99 (MQ=255)
                                                     aTCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCg                                                                                                                                           <  2:527610/99‑1 (MQ=255)
                                                        tttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGCGGTGGCCGCCGCCAccacca                                                                                                                                             >  2:1989162/1‑90 (MQ=255)
                                                        tttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGCGGTGGCCGCCGCCAccacca                                                                                                                                             <  1:1989162/94‑5 (MQ=255)
                                                         ttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGc                                                                                                    <  2:2141481/134‑1 (MQ=255)
                                                         ttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGc                                                                                                    >  1:2141481/1‑134 (MQ=255)
                                                              aTGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGcc                                                                               >  2:1008331/1‑150 (MQ=255)
                                                              aTGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGcc                                                                               >  2:575030/1‑150 (MQ=255)
                                                               tGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCg                                                                              >  2:1764817/1‑150 (MQ=255)
                                                                  cGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTg                                                                           >  1:582885/1‑150 (MQ=255)
                                                                  cGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTg                                                                           <  2:1964465/150‑1 (MQ=255)
                                                                  cGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGCGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTg                                                                           >  2:1615136/1‑150 (MQ=255)
                                                                         cACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGa                                                                    >  2:2202715/1‑150 (MQ=255)
                                                                          aCCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGaa                                                                   >  1:553307/1‑150 (MQ=255)
                                                                          aCCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGaa                                                                   >  2:2870/1‑150 (MQ=255)
                                                                           ccGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGaaa                                                                  >  2:931850/1‑150 (MQ=255)
                                                                            cGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCg                                                                                                                  <  1:1861784/101‑1 (MQ=255)
                                                                            cGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCg                                                                                                                  >  2:1861784/1‑101 (MQ=255)
                                                                            cGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGaaaa                                                                 >  2:2159478/1‑150 (MQ=255)
                                                                                 atatGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCCCCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGcc                                                            <  2:1108501/150‑1 (MQ=255)
                                                                                  tatGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCg                                                           <  2:123748/150‑1 (MQ=255)
                                                                                    tGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGCTAAAAGGCTGCCAGGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGc                                                         >  1:1365429/1‑150 (MQ=255)
                                                                                      ggCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCAt                                                       <  1:2202715/150‑1 (MQ=255)
                                                                                      ggCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACACTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGACATGTGGTCGGTGTGGCCGGTGGCGCGGAAAAATGCCGGCag                                                       >  1:1058733/1‑149 (MQ=255)
                                                                                        cTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGc                                                     >  2:1893905/1‑150 (MQ=255)
                                                                                            cTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCgg                                                                                                                 >  2:1849426/1‑86 (MQ=255)
                                                                                            cTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCgg                                                                                                                 <  1:1849426/86‑1 (MQ=255)
                                                                                            cTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTc                                                 >  1:556930/1‑150 (MQ=255)
                                                                                             tCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCg                                                >  1:1326979/1‑150 (MQ=255)
                                                                                               gATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAg                                              <  1:2870/150‑1 (MQ=255)
                                                                                                   ttGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGc                                          <  1:517677/150‑1 (MQ=255)
                                                                                                   ttGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGACCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGc                                          >  1:1341517/1‑150 (MQ=255)
                                                                                                      gccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAt                                                                                              <  2:1074929/95‑1 (MQ=255)
                                                                                                      gccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAt                                                                                              >  1:1074929/1‑95 (MQ=255)
                                                                                                      gccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCg                                                                                       >  2:1182519/1‑102 (MQ=255)
                                                                                                      gccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCg                                                                                       <  1:1182519/102‑1 (MQ=255)
                                                                                                       ccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCg                                                                              <  1:824076/110‑1 (MQ=255)
                                                                                                       ccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCg                                                                              >  2:824076/1‑110 (MQ=255)
                                                                                                               aGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGa                                               <  1:1758884/133‑1 (MQ=255)
                                                                                                               aGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCAGGTGGCGCGGAAAAATGCCGGCATGCGGTCGa                                               >  2:1758884/1‑133 (MQ=255)
                                                                                                                   gggTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTg                          >  1:137098/1‑150 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCCGGTGGGGGCGCTTCATGTCTTCCATGTGGTCGTTGTCGCCGGTGGTGCGGAAATATGTTTGCATTCGTTCGAGGTGCTGGGTTTTGATGCCTGTc                        <  1:551251/150‑1 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                        <  2:1823769/150‑1 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                        <  2:1823974/150‑1 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                        >  1:1823974/1‑150 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                        >  1:1823769/1‑150 (MQ=255)
                                                                                                                     gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                        <  2:1326979/150‑1 (MQ=255)
                                                                                                                         aaCCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGa                    <  1:1008331/150‑1 (MQ=255)
                                                                                                                              gcccggtggCCGCCGCCACCAGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCAcc               >  1:364457/5‑150 (MQ=255)
                                                                                                                                      gcggcggcCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCg       >  2:250868/1‑150 (MQ=255)
                                                                                                                                           ggcCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGACtt  <  2:582885/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGACTT  >  NC_002947/2822802‑2823090

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: