Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,569,756 T→G 7.4% K182Q (AAG→CAG)  nnrD ← ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,569,7560TG 7.4% 188.5 / 7.7 68K182Q (AAG→CAG) nnrDADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (27/36);  new base G (2/3);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

AGCACGCCAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTCGCGCCAGCAAATTCAAGGCATCGGCATC  >  NC_002947/5569611‑5569900
                                                                                                                                                 |                                                                                                                                                
aGCACGCCAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGAACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTgcgc                                                                                                                                              >  2:1644306/1‑150 (MQ=255)
   aCGCCAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCAAGCTGgcgcgcc                                                                                                                                           <  1:1050068/150‑2 (MQ=255)
       cAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCAc                                                                                                                                       >  2:588293/1‑150 (MQ=255)
       cAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCAc                                                                                                                                       <  1:45609/150‑1 (MQ=255)
          tCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCat                                                                                                                                    >  1:946668/1‑149 (MQ=255)
            aGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCAACAGGGTGCCGGCACCTTTGAGTACTCAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAggg                                                                                                                                  <  2:136149/150‑1 (MQ=255)
               acaTCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCTCCAGCTTGCGCGCAGCACCAGGGCgg                                                                                                                               >  1:822481/1‑150 (MQ=255)
               acaTCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCgg                                                                                                                               <  1:159772/150‑1 (MQ=255)
                 aTCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCTCGCGCGCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTc                                                                                                                             <  2:1621409/150‑1 (MQ=255)
                 aTCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTc                                                                                                                             >  2:890027/1‑150 (MQ=255)
                    gCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGc                                                                                                                          <  1:386348/150‑1 (MQ=255)
                     cccAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGcc                                                                                                                         <  2:1304875/150‑1 (MQ=255)
                     cccAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGcc                                                                                                                         >  1:955792/1‑150 (MQ=255)
                        aGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGt                                                                                                                      <  2:696699/150‑1 (MQ=255)
                          gccggcgccggcCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCAACTTTGAGTACACAGACGCTGTTATAACGGCGTGTCAGCTTGCGCGAAGCACGAGGGAGGTCAGCCTGTAc                                                                                                                    >  1:15148/1‑150 (MQ=255)
                               cgccggcCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCt                                                                                                               >  1:271143/1‑150 (MQ=255)
                                    gcCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTg                                                                                                          <  2:1505191/150‑1 (MQ=255)
                                       aTCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCa                                                                                                       >  2:206156/1‑150 (MQ=255)
                                            tGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGcc                                                                                                  >  2:622414/1‑150 (MQ=255)
                                                  cccGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCATGTCGGTCAGCCTGTACGGCCTCAGTGGCATTGccactgag                                                                                            <  2:946668/150‑7 (MQ=255)
                                                   ccGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGATGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCagcagc                                                                                           >  2:534932/1‑150 (MQ=255)
                                                    cgcgCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgcc                                                                                          >  2:174073/1‑150 (MQ=255)
                                                          cGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgccgcgccg                                                                                    <  2:1505951/150‑1 (MQ=255)
                                                                  gCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGt                                                                                                                              <  1:843495/100‑1 (MQ=255)
                                                                  gCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGt                                                                                                                              >  2:843495/1‑100 (MQ=255)
                                                                    cAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGAGCGCAGCACCAGGGCGGTCAGCCTGTAAAGCCTCAGTAACCATGCCCAGCAACCCCGCCGCCTccccctg                                                                          >  2:732979/1‑148 (MQ=255)
                                                                      gCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGt                                                                        <  1:890027/150‑1 (MQ=255)
                                                                       cTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTg                                                                       <  1:737474/150‑1 (MQ=255)
                                                                                 ggggCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGCCGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAgg                                                             <  1:588293/150‑1 (MQ=255)
                                                                                   ggCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGAt                                                           <  2:1640203/150‑1 (MQ=255)
                                                                                      tGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGAtgct                                                        <  2:478903/150‑1 (MQ=255)
                                                                                         caccaGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAAGGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGcc                                                     >  2:219817/1‑150 (MQ=255)
                                                                                          accaGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgccgcgcc                                                                                     >  2:1638515/1‑117 (MQ=255)
                                                                                          accaGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgccgcgcc                                                                                     <  1:1638515/117‑1 (MQ=255)
                                                                                           ccaGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCAc                                                   >  2:123146/1‑150 (MQ=255)
                                                                                             aGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTg                                                 <  1:206156/150‑1 (MQ=255)
                                                                                                   ccGGCACCTTTGAGTACACAGATGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAgg                                           <  1:534932/150‑1 (MQ=255)
                                                                                                   ccGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAgg                                           <  2:1466486/150‑1 (MQ=255)
                                                                                                   ccGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGAAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAgg                                           <  2:271143/150‑1 (MQ=255)
                                                                                                    cGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCcagca                                                                                             >  2:471651/1‑99 (MQ=255)
                                                                                                    cGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCcagca                                                                                             <  1:471651/99‑1 (MQ=255)
                                                                                                    cGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgccgcgc                                                                                      >  1:728130/1‑106 (MQ=255)
                                                                                                    cGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAgccgcgc                                                                                      <  2:728130/106‑1 (MQ=255)
                                                                                                    cGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAggg                                          >  2:1049867/1‑150 (MQ=255)
                                                                                                     ggCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCca                                                                                                >  1:810227/1‑95 (MQ=255)
                                                                                                     ggCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCca                                                                                                <  2:810227/95‑1 (MQ=255)
                                                                                                        aCCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTaa                                                               >  2:1338485/1‑125 (MQ=255)
                                                                                                        aCCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTaa                                                               <  1:1338485/125‑1 (MQ=255)
                                                                                                        aCCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGAtgc                                                         >  1:1094188/1‑131 (MQ=255)
                                                                                                        aCCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGAtgc                                                         <  2:1094188/131‑1 (MQ=255)
                                                                                                          cTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCg                                    <  1:870017/150‑1 (MQ=255)
                                                                                                            ttGAGTACACAGACGCTGTTATACCGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCggg                                  <  1:520552/150‑1 (MQ=255)
                                                                                                            ttGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGc                                     <  2:531391/147‑1 (MQ=255)
                                                                                                            ttGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGc                                     >  1:531391/1‑147 (MQ=255)
                                                                                                            ttGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCggg                                  <  1:123146/150‑1 (MQ=255)
                                                                                                              gAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGt                                >  1:1323823/1‑150 (MQ=255)
                                                                                                                  acacAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTcgc                            <  2:687818/150‑1 (MQ=255)
                                                                                                                   cacaGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTcgcg                           >  2:438040/1‑150 (MQ=255)
                                                                                                                       gACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTCGCGCCAg                       <  2:955792/150‑1 (MQ=255)
                                                                                                                                ataACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTCGCGCCAGCAAATTCaa              <  1:174073/150‑1 (MQ=255)
                                                                                                                                 taACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTgg                                                                                                         >  1:773123/1‑58 (MQ=255)
                                                                                                                                 taACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTgg                                                                                                         <  2:773123/58‑1 (MQ=255)
                                                                                                                                    ctgcgcgcAAGCTGGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCgggg                                 >  2:972629/3‑127 (MQ=255)
                                                                                                                                    ctgcgcgcAAGCTGGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCgggg                                 <  1:972629/125‑1 (MQ=255)
                                                                                                                                     tgcgcgcAAGCTGGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCa                                             <  1:641743/113‑1 (MQ=255)
                                                                                                                                     tgcgcgcAAGCTGGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCa                                             >  2:641743/2‑114 (MQ=255)
                                                                                                                                            cAGCTTGCGCGCAGCACCCGGGCGGTAAGCCTGTACGGCCTGAGTGGCAATGCCCAGCAGCCGCGCTGCCTCCCCAGGCTGCGCAGTACGGATGCTGCAACTGGGAAGGGCCAGCGGGATTTTCGCCATGAAATTCAAggcatcggcatc  <  2:431047/150‑1 (MQ=255)
                                                                                                                                            cAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTCGCGCCAGCAAATTCAAggcatcggcatc  <  1:438040/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGCACGCCAGTCAGCACATCGCCCAGGCCGGCGCCGGCCATCGCTGGGTGCCCGCGCTCGCACAACGCCAGCTGCCCGCTCGGGGCTGCCACCAGGGTGCCGGCACCTTTGAGTACACAGACGCTGTTATAACGGCGCGCCAGCTTGCGCGCAGCACCAGGGCGGTCAGCCTGTACGGCCTCAGTGGCAATGCCCAGCAGCCGCGCCGCCTCCCCCGGGTGCGGAGTAAGGATGCTGCCACTGGGCAGGGCCAGCGGGGTTCGCGCCAGCAAATTCAAGGCATCGGCATC  >  NC_002947/5569611‑5569900

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: