Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,957,897 G→C 8.1% P219P (CCG→CCC PP_2589 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,957,8970GC 8.1% 193.4 / 7.5 62P219P (CCG→CCCPP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (3/2);  total (32/30)
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.83e-01

TGGCCGCCTGGAAGCTGGCCCCGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGG  >  NC_002947/2957752‑2958043
                                                                                                                                                 |                                                                                                                                                  
tGGCCGCCTGGAAGCTGGCCCCGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGc                                                                                                                                                 <  2:62429/149‑1 (MQ=255)
          gAAGCTGGCCCCGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGATGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAg                                                                                                                                        >  2:180757/1‑148 (MQ=255)
              cTGGCCCCGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGc                                                                                                                                   >  8:36710/1‑149 (MQ=255)
                   cccGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCggg                                                                                                                               >  5:239516/1‑148 (MQ=255)
                   cccGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCggg                                                                                                                               >  7:215931/1‑148 (MQ=255)
                   cccGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCggg                                                                                                                               >  2:167741/1‑148 (MQ=255)
                   cccGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCCGTCCTGCTCGCCCAGGCCCGGc                                                                                                                              <  8:288694/149‑1 (MQ=38)
                    ccGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCggg                                                                                                                               >  7:288493/1‑147 (MQ=255)
                             gCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgg                                                                                                                     >  1:50870‑M2/1‑146 (MQ=255)
                                 ccGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCt                                                                                                                 <  3:137218/148‑1 (MQ=255)
                                 ccGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCt                                                                                                                 <  4:77136/148‑1 (MQ=255)
                                 ccGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTg                                                                                                                >  3:263038/1‑149 (MQ=255)
                                    gTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGc                                                                                                                                                 <  5:340617/113‑1 (MQ=255)
                                    gTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGc                                                                                                                                                 >  6:340617/1‑113 (MQ=255)
                                       accgCGTGGTCCTCAAGCCCGCCGAGCACTCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCGGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCAccc                                                                                                           <  8:288493/145‑1 (MQ=255)
                                               gTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGa                                                                                                  <  7:36710/149‑1 (MQ=255)
                                                 aCTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACt                                                                                                >  8:310394/1‑149 (MQ=255)
                                                   tCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGc                                                                                              >  4:264403/1‑149 (MQ=255)
                                                     aaGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGc                                                                                              <  3:270305/147‑1 (MQ=255)
                                                        cccGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCGGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTg                                                                                           >  5:317981/1‑147 (MQ=255)
                                                        cccGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTg                                                                                           <  6:317981/147‑1 (MQ=255)
                                                        cccGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGt                                                                                         <  3:264403/149‑1 (MQ=255)
                                                        cccGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGt                                                                                         <  5:247532/149‑1 (MQ=255)
                                                              gAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCAc                                                                                   <  1:223121/149‑1 (MQ=255)
                                                              gAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCAc                                                                                   <  7:352656/149‑1 (MQ=255)
                                                                  aaTCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGgcg                                                                                                                         >  7:126563/1‑107 (MQ=255)
                                                                  aaTCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGgcg                                                                                                                         <  8:126563/107‑1 (MQ=255)
                                                                  aaTCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTCGAAGCCGGTGTGCCGGAAGGCGGGATGAACGTGCTGCCGGGCCTGGGCGAGCAGGCCGCGCAGGCGCTGGGCCTGCACCCGGATGAGGACTGCCTGGTATTCACCgg                                                                                >  5:206304/1‑148 (MQ=255)
                                                                      aCCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGATGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTcc                                                                            <  1:180757/148‑1 (MQ=255)
                                                                      aCCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTc                                                                                     >  7:133327/1‑139 (MQ=255)
                                                                      aCCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTc                                                                                     <  8:133327/139‑1 (MQ=255)
                                                                              cGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGt                                                                    <  4:263038/148‑1 (MQ=255)
                                                                              cGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGt                                                                    <  7:310394/148‑1 (MQ=255)
                                                                                 cccTGCGCCTGGCCGAGCTGGGGCTGGAAGGCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGGCGGGGCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCTGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGc                                                                <  6:206304/149‑1 (MQ=255)
                                                                                  ccTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgggccgaaaacggtgattgctcggcgggcttgagtaccacgctgttaccggcagc                                                                <  1:236240‑M2/148‑56 (MQ=255)
                                                                                  ccTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCAc                                                                          <  6:217948/138‑1 (MQ=255)
                                                                                  ccTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCAc                                                                          >  5:217948/1‑138 (MQ=255)
                                                                                   cTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCaa                                                              <  1:167741/149‑1 (MQ=255)
                                                                                        ccTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgggccgaaaatggtgattgctcggcgggcttgagtaccacgctgttaccggcagccagcgcc                                                         <  3:246299‑M2/149‑63 (MQ=255)
                                                                                             cccAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATg                                                    <  6:239516/149‑1 (MQ=255)
                                                                                             cccAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATg                                                    >  3:59789/1‑149 (MQ=255)
                                                                                              ccAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATg                                                    >  4:87169/1‑148 (MQ=255)
                                                                                              ccAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATg                                                    >  1:209170/1‑148 (MQ=255)
                                                                                                 gCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgggccgaaaacggtgattgctcggcgggcttgagtaccacgctgttaccggcagccagcgccggggccagc                                                >  2:236240‑M2/1‑78 (MQ=255)
                                                                                                       gcTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgg                                                                                                                     >  8:305142‑M2/1‑72 (MQ=255)
                                                                                                       gcTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCgg                                          >  4:257745/1‑149 (MQ=255)
                                                                                                           gAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGgcg                                        >  7:263875/1‑147 (MQ=255)
                                                                                                                        cGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAg                         <  2:209170/149‑1 (MQ=255)
                                                                                                                         ggAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTggg                                                                 <  1:359822/108‑1 (MQ=255)
                                                                                                                         ggAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTggg                                                                 >  2:359822/1‑108 (MQ=255)
                                                                                                                            aGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTgggccgaaaatggtgattgctcggcgggcttgagtaccacgctgttaccggcagccagcgccggggccagcttccaggcggccatatcgagcgggaag                     >  4:246299‑M2/1‑51 (MQ=255)
                                                                                                                                 tGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTCCTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTgg                >  4:227302/1‑149 (MQ=255)
                                                                                                                                   cTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCt              >  1:72897/1‑149 (MQ=255)
                                                                                                                                      aaCGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTg             >  1:104951/1‑147 (MQ=255)
                                                                                                                                          gcctgCCCGGCCTGCTCGCCCAGGCCCGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTgc       >  2:271620/4‑149 (MQ=255)
                                                                                                                                           cctgCCCGGCCTGCTCGCCCAGGCCCGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTgcg      >  4:27236/3‑149 (MQ=255)
                                                                                                                                             tgCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGgc                                         <  3:367172/112‑1 (MQ=255)
                                                                                                                                             tgCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACACGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTAATTCATGCAGTACTCGgc                                         >  4:367172/1‑112 (MQ=255)
                                                                                                                                              gCCCGGCCTGCTCGCCCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTc                            >  3:115121/1‑124 (MQ=255)
                                                                                                                                              gCCCGGCCTGCTCGCCCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTc                            <  4:115121/124‑1 (MQ=255)
                                                                                                                                               ccGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGcggcgg  <  8:263875/149‑1 (MQ=255)
                                                                                                                                               ccGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGcggcgg  <  3:87169/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGCCGCCTGGAAGCTGGCCCCGGCGCTGGCTGCCGGTAACAGCGTGGTACTCAAGCCCGCCGAGCAATCACCGTTTTCGGCCCTGCGCCTGGCCCAGCTGGCGCTGGAAGCCGGTGTGCCGGAAGGCGTGCTGAACGTGGTGCCGGGCCTGGGCGAGCAGGCCGGGCAGGCGCTGGGCCTGCACCCGGATGTGGACTGCCTGGTATTCACCGGCTCCACCCAGGTGGGCAAGTACTTCATGCAGTACTCGGCGCAGTCCAACCTCAAGCAGGTGTGGCTGGAGTGCGGCGG  >  NC_002947/2957752‑2958043

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: